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Introduction 



Charmonium, the bound system of charmed quark and antiquark, is studied experimen- 
tally since 1974, when the J/'4> meson was discovered by two groups at SLAC and BNL. 
Since that time eight charmonium states below the mass of open charm {DD) threshold 
were observed. The theoretical description of charmonium started with the non-relativistic 
potential model, which was motivated by the large mass of the c quark. Improved by incor- 
poration of spin-dependent terms and relativistic corrections it provides a good agreement 
with the mass spectra of charmonium below the DD threshold. The recent progress in 
approaches based on Lattice QCD and Effective Field Theory (NRQCD) also allowed a 
good theoretical description of the charmonium system. However, there are still many open 
question in this area, which can only be answered with high precision and high luminosity 
measurements. 

The PANDA experiment, which is a part of the future Facility for Antiproton and 
Ion Research (FAIR) at Darmstadt, will provide possibility for the precise spectroscopy 
of charmonium produced in antiproton-proton annihilations. The cooled antiproton beam 
circulated in the High Energy Storage Ring (HESR) with a momentum between 1.5 GeV/c 
and 15 GeV/c will provide a peak luminosity up to 2 • 10^^ cm~^s~^. In contrast to the 
e^e~ annihilations which excite only the states with quantum numbers of the photon, 
J^^" = 1 , in the pp collisions the charmonium states with all the quantum numbers can 
be formed directly. Therefore their masses and widths can be measured with high precision 
determined by the parameters of the beam and not limited by the detector resolution. 

The subject of this thesis is the study of electromagnetic transitions in charmonium 
with PANDA. The possible registration of the he state in charmonium, observed recently 
by the E760 and E835 Fermilab experiments, is used for demonstration of the physics 
performance of the PANDA detector. The measurement of the angular distribution of 7- 
rays from radiative transitions can be used for the verification of its J^*" quantum numbers. 
The electromagnetic calorimeter of the PANDA detector is the crucial component for this 
studies from the point of view of signal registration and background suppression. 
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CONTENTS 



In Chapter [T] an overview of the charmonium physics is given, which includes the 
status of experimental studies as well as development of the theoretical description of 
charmonium. Chapter |2] gives an overview of the current design status of the PANDA 
detector. The electromagnetic calorimeter (EMC) is described in more detail because of 
its interest for the presented work. 

In Chapter [3] the results of measurements of energy resolution of PhWO^ scintillators 
with the Avalanche Photodiode (APD) readout are presented. Low energy protons and 
7-rays produced in the radiative capture reaction 7)^^(7 were used for the response 

studies. 

Chapter |4] presents the results of Monte Carlo simulations which permit to evaluate the 
performance of the PANDA detector in clarifying several questions related to physics of 
the he charmonium state. Chapter 5 contains conclusions. 



Chapter 1 



Physics motivation 



1.1 Physics of charmonium 
1.1.1 Quark model and QCD 

The Standard Model of particle physics is f/(l) ® SU{2) ® SU{3) field theory and according 
to it all the matter is composed from the 6 quarks and 6 leptons divided into 3 generations 
and gauge bosons corresponding to 3 groups of symmetry of the Standard Model. One 
of the main components of the Standard Model is Quantum Chromodynamics or QCD. 
QCD describes interaction between quarks mediated by massless gauge bosons - gluons. 
QCD is responsible for 99% of the mass of proton and as a result of the baryonic matter 
of Universe. 

The notion of quark was introduced before establishing of QCD in quark model by 
Murray Gell-Mann in 1964 |6T], who proposed together with George Zweig the scheme of 
classification of mesons and baryons known at that time. The lightest mesons and baryons 
were classified according to the multiplets of SU (3) fiavour group in the following way. In 
case of the three fiavours of quarks they lie in the fundamental representation, 3 (called 
the triplet) of fiavour SU(3). The antiquarks lie in the complex conjugate representation 
3*. Mesons consisted from a quark and an anti-quark and nine states made out of quark- 
antiquark pairs can be decomposed into a trivial representation (singlet) and an adjoint 
representation (octet): 

3®3* = 8©1, (1.1) 



which are filled with the 9 lightest spin mesons (Fig. 1.1). 

In case of baryons, the basis quark states can be considered as the six states of three 
fiavours and two spins per fiavour (the so called approximate spin-fiavour SU(6) symmetry). 



7 



8 



1.1 Physics of charmonium 



TT- 




-1 _ , -1/2 

— ua 



Figure 1.1: Octet of the hghtest spin mesons, located on the plane strangeness (S) - isospin 
Is projection. Q is charge axis. 



The spin-flavour part of a barion wavefunction is decomposed in the following way: 

6 ® 6 ® 6 = 565 © 70m © 70m © 20a, (1-2) 

where the 56 symmetric states are decomposed under flavour SU(3) into: 

56 = lot ©8^; (1.3) 

i.e. a decouplet of spin-3/2 and an octet of spin-1/2 baryons. The octet of lightest spin-1/2 
baryons in presented in Fig. [1.2 

The quantum number of color was proposed for an explanation of the baryon 
properties [HB]- It has spin 5 = 3/2 and is composed of the three u quarks with parallel 



spins and relative angular momentum. In order to have an antisymmetric total wave- 
function to conform with Pauli principle there should be a hidden quantum number, which 
was called color. It can have three values: red, green and blue and the combination of 
all three colors gives white, i.e. a color neutral object. This corresponds to the case of 
baryons, where all three quarks have different color quantum numbers. In case of mesons 
each of them is built from the quark- ant iquark pair in which a quark carries color and an 
anti-quark carries the corresponding anti-color. As a result all the observable particles are 
color-neutral. 

The existence of three colors came after the measurements of the ratio of electron- 
positron annihilation into hadrons and muon pairs: 

a{e~^e~ — )■ hadrons) 



R 
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Figure 1.2: Octet of the lightest spin 1/2 baryons, located on the plane strangeness (S) - isospin 
I3 projection. Q is charge axis. 



In the lowest order it is equal to N^.qf, where is the number of colors and the sum 
is taken over quarks with 2mg < E, E - energy of the process. In the energy region below 
the 6-quark production, the sum is equal ^ for = 3, which is in good agreement with 
experimental data. 

The additional step between discovery of quarks, color quantum number and the QCD 
is establishing of dynamical role of color octet gluons, which was proposed by Fritzsch and 
Gell-Mann [S^. Although the word Lagrangian was not mentioned in their work, they 
made complete description of two terms that make up the QCD Lagrangian. 

QCD is the relativistic quantum field theory of strong interactions and describes inter- 
action of quarks and gluons according to the QCD Lagrangian: 

^QCD = -\g>:GI^ + E^t^^'^M - ^fMf^ (1-5) 

/ 

where 

Gr = d'Al - d^A^, + gsf^A^.A: (1.6) 
is the gluon field strength tensor, and 

D^^ = Q^^_ i^A^X" (1.7) 

is the gauge covariant derivative, gs is the gauge constant, / denotes the quark flavour, 7^^ 
are the Dirac 7- matrices, A'^ correspond to the gluon field with a running from 1 to 8, A" 
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correspond to 3 x 3 matrices - generators of SU(3) group, f^'^ are structure constants of 
SU (3) group and fi, v are Lorentz indexes running from 1 to 4. 

QCD has some similarities to QED with the essential differences coming from it being 
non-abelian theory (SU(3) group, which is a base group of the QCD symmetry is non- 
abelian in contrast to U(l) symmetry of QED). This property made implications to the 
basic interaction vertex. In QED the vertex connects only two fermions (electrons) with a 
gauge boson (photon), whereas in QCD in addition to a two quark-gluon vertex there are 
tree- and four-gluon interaction vertices, which describe interaction between gluons. 

The gauge constant is related to an effective coupling constant of the strong interaction 

by: 



The coupling constant depends on the energy scale /x and in the lowest order can be 
expressed as: 

In 

"^^^^ = [ll-(2/3)n,]ln(/./A)' ^^'^^ 

where n/ - number of quark flavours with the mass less than yU and A - QCD scale parameter. 
At the mass of Z° boson a,(Z°) = 0.1172 ± 0.0020. 

As can be seen from the above expression (Eq. 1.9) and from the measured values 
(Fig. 1.3), Us decreases for higher energies or small distances and approaches zero at /i — )■ 
oo, a property called the asymptotic freedom. 

Although the Feynman rules can be derived from the QCD Lagrangian (Eq. 1.5), such a 
technique can be used for the calculation of strong processes only at large momenta, where 
coupling constant becomes very small. At small momenta, and correspondingly large 
coupling constant the perturbative methods do not work, in that case different effective 
theories or numerical solutions such as Lattice QCD are used to calculate strong processes. 



Their application for description of the charmonium is discussed in Sec. 1.2 



1.1.2 History of charmonium physics 

The review of charmonium physics outlined below is based on the review articles by Bianco 
et al.|3^, which covers theory of charmonium, Bettoni and Calabrese [29] which describes 
the current status of experiments, QWG report [34J which covers both theoretical and 
experimental aspects of charmonium physics, reviews given in several PhD theses (e.g. 
[50] , |67j ) on charmonium physics and on original papers which are referred to in the text. 

The first proposal for existence of the fourth quark flavour appeared in 1964 in the 
work of Bjorken and Glashow [32]. They proposed the quark flavour which they termed 



1 Physics motivation 



11 



T \ — I I I I I I I 




1 10 10 

fiGeV 

Figure 1.3: Summary of the measured values of a<j(/i) [1]. 

"charm" to combine quarks in doublets by analogy with lepton doublets (e~,z/e), (/i",!/^) 
and as a consequence extension of the SU(3) flavour group to SU(4). The argument for 
such an extension comes from the requirement of renormalizability of the theory. The 
triangle diagram, which is represented by fermion loop with three external spin-one lines, 
creates the so-called Adler-Bell-Jackiw anomaly. The non-conservation of an axial current 
produced by this diagram leads to infinities to higher orders. The only way to get rid of 
them is to cancel them among the different fermion loops. In the Standard Model this 
requires that electric charges of all the fermions - quarks and leptons should add up to 
zero. For electron, muon, up, down and strange quarks their charges add up to -2. This 
means that the fourth quark with the three colors is needed with the charge +| to make 
the sum equal 0. 

The experimental evidence for the existence of the fourth quark came from the Cabibbo 
theory [38] • He proposed a model which explained suppression of the strangeness-changing 
semi-leptonic weak decay {AS = 1) with respect to strangeness conserving decay {AS = 0). 
According to his idea in weak interactions the flavour eigenstate d does not participate but 
rather a mixture of d and s quarks, dc = dcos{6c) + ssin{6c). The angle 6c is called Cabibbo 
angle with an experimental value of 0.25. However, this model predicted the existence of a 
neutral current with AS* = 1, which leads, in particular, to the contribution to Kl — )■ fi~^fi~ 



12 



1.1 Physics of charmonium 



several orders of magnitude larger than experimentally observed. 

The solution to this problem was proposed by Glashow, Iliopoulos and Maiani 
in 1970. They proposed the existence of the fourth quark with a charge +2/3 and a 
flavour called charm which participates in weak interaction in a doublet with the state 
Sc = scos{9c) — dsm{9c) orthogonal to the dc state; i.e. the two quark doublets which 
participate in weak interaction are: 



u 

dc = dcosOc + ssinOcJ \Sc = scosOc — dsinOc 



(1.10) 



Due to dc and Sc orthogonality the strangeness-changing (AS* = 1) term in the neutral 
current, 

J^+^ = . . . + {scdc + ~s~cdc) sin Oc cos 9c , (l-H) 

vanishes. 

But in spite of these arguments, the existence of the fourth quark and validity of the 
quark model itself became widely accepted only with the discovery of the J/ip meson in 
1974. The J/ip was discovered almost simultaneously by two groups from the USA, neither 
of which was searching for charm. 

The group of Richter at SLAC's SPEAR e^e~ collider, reported the resonance called ip 
which was observed as a peak in cross-section of the reaction e+e~ — e~ , fi'^ fi~ , hadrons 
[21] . The reported mass of the resonance was 3.105 ± 0.003 GeV and its width F < 1.3 



MeV. The corresponding discovery plot is presented in Fig. 1.4 (a). The three plots top 
to bottom show enhancements in the cross-section for e^e~,fi'^fi~ and hadron final states, 
respectively. 

In the Brookhaven National Laboratory (BNL) the group led by Ting reported a particle 
called J in the e~^e~ invariant mass in the reaction 

p + Be^e^e' +X (1.12) 

in the fixed target experiment at 28 GeV [20j. The reported mass of the enhancement is 



M(e+e~) = 3.1 GeV. The plot is presented in Fig. 1.4 (b). 

Immediately after announcing the discovery of new states by these two groups the 
ADONE e^e~ collider at Frascati, which was designed for a maximum center of mass 
energy of 3.0 GeV, boosted current in their magnets over the design limits and confirmed 
existence of the 3.1 GeV resonance [22]. The papers announcing new resonance discovery 
appeared in the December 1974 issue of the Physical Review Letters and the particle has 
been henceforth called J/ip. Another confirmation of J/ip came from DESY in the wide 
angle elastic e~^e~ scattering experiment [35] . 
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Figure 1.4: J/tp discovery plots: by the SLAC group in e~'~e collisions [21] (a) and by the BNL 
group in p + Be — )■ e+e" + X reaction [20] (b). 
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In all the quoted experiments the observed width of the peak was dominated by ex- 
perimental resolution. The proper width of the J/ip state could be determined from the 
measurements of the total integrated cross-section and the leptonic branching ratio, both 
of which were measured experimentally. From the Breit-Wigner formula the J/ip width has 
been estimated at the level of one hundred keV, which is much smaller than the MeV width 
dominated by beam resolution. Most resonances in strong interaction known at that time 
were much wider, having widths of the order of MeV and larger. The narrow width of the 
discovered J/ip particle led to its interpretation as a bound system of the fourth "charm" 
quark- ant iquark pair cc. Such an interpretation was proposed by Appelquist, and Politzer 
[8] and De Rujula and Glashow [53]. 

The small width of J/ip can be interpreted in terms of the OZI (Okubo-Zweig-Iizuka) 
rule [HE], |121j . [77] . The OZI rule states that strong processes, which can be described by 
Feynman diagrams containing only disconnected quark lines between the initial and final 
states are strongly suppressed in comparison with diagrams with connected quark lines. 
Two examples related to the processes with a charmed quark-antiquark pair are presented 



in Fig. 1.5 The left diagram corresponds to the OZI-allowed decay of the charmonium 
state. It is assumed that the mass exceeds the DD threshold, which corresponds for exam- 
ple to the Tp" decay (see Fig. 1.7). The right diagram corresponds to the OZI-suppressed 
J/ip decay into light hadrons. 

The OZI rule can be explained qualitatively in the following way. Processes with 
disconnected quark lines in the initial and final states are connected via gluons. All the 
intermediate gluons must combine in a way to preserve all the strong interaction quantum 
numbers. Since a gluon carries a color quantum number and mesons are colorless there 
should be at least 2 intermediate gluons, connecting the initial and final states. On the 
other hand vector mesons, such as J/ip, have charge parity quantum number C = —1, which 
requires an odd number of gluons. As a result, at least 3 intermediate gluons participate in 
J/ip strong decay. Since the decay mesons are massive, the gluons must be very energetic 
("hard") and the coupling constant is small in comparison with the OZI-allowed decay 
where gluons can be "soft" and the coupling constant is large. 

The first radial excitation of charmonium, ip', was found at SLAG [2] soon after the J/ip 
observation. This state also appeared to be very narrow with the observed width consistent 
with beam resolution 1 MeV). Because the mass of ip' is elevated (3686 MeV), but still 
below the DD threshold, the narrowness of the resonance was also interpreted as a result 
of the OZI rule. 

Other charmonium states {xco, Xa, Xc2, Vc) were measured at SLAG with the series of 
improved detectors Mark I, Mark II, Mark III and Grystal Ball. The Xcj and rjc states 
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Figure 1.5: Feynman diagrams demonstrating OZI rule. The left diagram demonstrates the 
OZI-allowed decay of a charmonium state into a D^D^ pair, right diagram corresponds to an 
OZI-suppressed decay of a charmonium state into light hadrons. 



observed in the radiative decay of ip' in the Crystal Ball experiment |102] are presented in 



Fig. 1.6 



The radial excitation of the charmonium ground state, rj'^, was firmly established in 
2002 by the Belle experiment j37] as a peak in the invariant mass distribution of KsK~tt^ 
with the mass 3654 MeV in a sample of exclusive B — )■ KKsK~tt^ decays collected with 
the Belle detector at the KEKB asymmetric e^e~ collider. The 1]'^ was confirmed soon 
afterwards by the CLEO [TB] and Babar [12] experiments. 

The he is the last among the observed discrete charmonium states. This is a major topic 



of the present thesis, therefore a separate chapter (Sec. 1.3.4) is devoted to the history of 
its discovery. 



Fig. 1.7 shows the spectrum of the charmonium states below the DD threshold discov- 
ered so far. The states of charmonium are classified according to the standard spectroscopic 
notation n^^^^Lj, where n numerates the radial excitations and L denotes the orbital an- 
gular momenta S, P, D, ... for L = 0, 1, 2, ... respectively. Each constituent quark has 
spin ^ therefore the total spin quantum number can be either S=0 or 1. J is the quan- 
tum number of the total momentum, a vector sum of the spin and orbital momenta. An 
alternative notation used to denote the charmonium states is J^*", where P is the space 
parity and C - is the charge conjugation parity. For a quark-antiquark system these can 
be expressed as: 

P=(-l)^+\ (1.13) 

and 

^=(-1)^+^, (1.14) 



respectively. 
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According to Eqs. 1.13 and 1.14 the states with natural spin-parity P = (—1)'^ should 
have 5* = 1 and CP = +1, respectively. The states with natural spin-parity and CP = —1 
(0^~, 1^^, 2+~, etc) are forbidden in the qq' model as well as the state J^'" = . 
However, exotic quantum numbers can be realized for exotic mesons (glueballs, hybrids, 
diquark-antidiquarks, etc). 

All the known charmonium states below the DD threshold with their masses and widths 
are summarized in Table 11.11 

The first two discovered charmonium states had J^'" = 1 , i.e. the quantum num- 
bers of the photon. Only these states can be directly formed in e~^e~ annihilations. The 



Feynman diagram for the corresponding process is shown in Fig. 1.8 All the other states 
accessed in e~^e~ annihilations are reached in the decay of such vector states, in particu- 
lar via radiative transitions. For example, the energy balance in the cascade transitions 
V^' Xco,i,2 + 7i — ^ J/i' + 7i + 72, via the putative states Xcj permitted to assign them 
to charmonium in the Crystal Ball experiment |102j (see Fig. 1.6). Other charmonium 
states, besides J^^ = 1 , can in principle be obtained in e^e~ annihilations with two or 
three intermediate photons. However, a process with two intermediates photons contains 
additional factor in probability and correspondingly the event rate a factor 10~^ smaller 
than in processes with one intermediate photon. The technique using e~^e~ collisions was 
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Notation 




jPC 


Mass, MeV 


Width, MeV 


Vc 


I'So 


0-+ 


2980.3 ± 1.2 


26.7 


J/^ 




1 — 


3096.916 ±0.011 


0.0932 


XcO 


l^Po 


0++ 


3414.75 ±0.31 


10.2 


Xcl 




1++ 


3510.66 ±0.07 


0.89 


Xc2 




2++ 


3556.20 ± 0.09 


2.03 


he 




1+- 


3525.93 ±0.27 


< 1 




2'So 


0-+ 


3637 ± 4 


14 


^' 




1 — 


3686.09 ± 0.04 


0.317 



Table 1.1: Charmonium states below the DD threshold. 




Figure 1.8: Formation of the vector (1 ) charmonium in e'^e annihilations, decaying into 
lepton pairs. 



the only means to study charmonium until the middle of the 80-ties. 

An alternative mechanism of charmonium production in proton-antiproton annihila- 
tions for charmonium spectroscopy was demonstrated for the first time in the R704 exper- 
iment [23] at CERN. This technique made possible the formation of states with all J^*" 
values via intermediary of an appropriate number of gluons or qq pairs. An example of 



the process with two or three gluon annihilation is shown in Fig. 1.9 and corresponds to 
formation of the states with C = ±1 an C = — 1, respectively, due to negative internal 
charge parity of a gluon. This technique became successfully applicable thanks to the de- 
velopment of stochastic cooling of the antiproton beam. In R704 for the first time the '", 
1"*"+ and 2++ charmonium states were formed directly. The technique developed at CERN 
was later successfully used at Fermilab in the two subsequent experiments E760 and E835 
|60j . A significant improvement has been achieved in resolution of non — 1 charmonium 
states. The main difficulty in application of the pp technique in charmonium spectroscopy 



1 Physics motivation 



19 



g 



sQOQOQQQQ / 



ill 



OQQOQOQQQ 



n P 



Figure 1.9: Formation of charmonium in 2- and 3-gluon pp annihilation. 



is dealing with huge hadronic background. The total pp cross-section in the energy range 
of charmonium is nearly 60 mb, whereas the highest peak cross-section for an individual 
charmonium state production (J/ip) is about 270 nb [T3] . 

It is worthwhile to mention here the states above the DD threshold at 3.73 GeV. 
This region is rather poorly known. The first experiments with e~^e~ collisions reported 
a number of wide vector states, which were assigned to ■ip{2S), ipi'^S), ip{4:S) [58] but not 
all of them were confirmed in a much more precise BES experiment. In addition to vector 
states this region should contain first radial excitations of the singlet and triplet vector 
states as well as the narrow D-states. Only one of them, \Ef(3770), which is considered to 
be largely ^-Di, has been observed so far |106] . 

Moreover, in the last years in analyses of high statistics data samples by CLEO, Babar 
and Belle experiments a number of new states labeled X, Y and Z were discovered with 
masses in the range of charmonium [65]. Some of them can be interpreted as charmonium 
states, for example X(3872) has a possible interpretation as 1^D2 or l^-Ds with alternative 
interpretations as charmonium hybrid or a D^D^* molecule. The other two states Z(3931) 
and X(3940) have the most probable interpretation as Xc2(2-P) and ric{3S), respectively, 
while the nature of Y(3940), Y(4260), Y(4320), Y(4660) and Z(4430) is an open question. 
To conclude, the situation above DD still needs further work to be fully understood. 



1.2 Theoretical approaches to charmonium 

In description of charmonium spectra two different approaches are used: phenomenological 
and QCD-oriented. The phenomenological approach tries to model the main features of 
QCD which are relevant to describe the spectrum of charmonium; it is represented mainly 
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by potential models. The QCD-oriented approaches use QCD-based calculations, whether 
directly (lattice QCD) or within different approximations (effective theories). 

1.2.1 Phenomenological approach. The potential model 

QCD motivated potential models have been used for describing charmonium spectrum since 
the discovery of charmonium states. However, their relation to QCD is rather loose and 
the model quantities like quark mass, potential parameters etc. cannot be reliably related 
to the Standard Model parameters. Nevertheless, the application of potential models 
to the physics of charm was rather beneficial and helped a lot in interpretation of the 
experimental results. Non-relativistic potential models are based on the assumption that 
the quark-antiquark interaction can be approximated with the aid of a static potential. 
Such an assumption is based on the expectation that for sufficiently heavy quarks the 
characteristic time scale associated with the relative motion of the constituent quarks is 
much larger than that associated with the gluonic or the sea quark degrees of freedom. In 
this case the adiabatic approximation applies and the effect of gluon and sea quarks can 
be represented by an average instantaneous interaction potential between the heavy quark 
sources. Despite some doubts on the existence of a static potential between two quarks, 
calculations based on this assumption were successful in predicting the masses of bound 
states of charmonium. The non-relativistic treatment of charmonium is possible due to 
relatively high charmed quark mass, which leads to f^/c^ ~ 0.25 in the bound states. 
The non-relativistic potential model is based on the Schrodinger equation: 



where T is the kinetic energy term and V{r) is the potential energy term. The potential 
should reproduce two main properties of QCD: the asymptotic freedom at small distances 
and the confinement at large distances. The quark-gluon interaction in QCD is similar 
to the electron-photon interaction in QED with the Born term for qq interaction of the 
Coulomb form at small distances. This was first recognized in 1975 by Appelquist and 
Politzer [S]. The difference from QED is that the gluon self-coupling results in a slow 
decrease of the strong coupling constant at small distances where quark-antiquark potential 
due to one gluon exchange has the form: 




V{r) 




(1.16) 



3 r 



with the factor 4/3 arising from SU(3) color factors. 
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Appelquist and Politzer proposed the term "orthocharmonium" for the J/tp state in 
analogy with ^Si orthopositronium and predicted the existence of ^5*0 paracharmonium. 
Later they were able to predict [IDj the complete spectrum of bound states of charmonium 
based on the potential intermediate between Coulomb and harmonic oscillator in its radial 
dependence. 

However the most popular among the potential models became the Cornell potential 
proposed by Eichten et al. in 1975 [55], which combines one-gluon exchange Coulomb term 
and the additional term proportional to r, known as the confinement term. At distances of 
roughly 1 fm, which correspond to energy scale of Aqcd the as increases significantly and 
one-gluon-exchange is no longer a good approximation. Qualitatively the chromoelectric 
lines of force can be imagined to bunch together into a flux tube which leads to a distance 
independent force or the potential V{r) = ar. Such a term satisfies confinement, i.e. the 
experimental fact that colored free quarks are never experimentally observed. The Cornell 
potential has the form: 

V{r) = — + -, (1.17) 

r 

where the constant 1/a^ is of order ~ 1 GeV/fm. 

Soon afterwards many alternative formulations of the quark- ant iquark potential ap- 
peared. Quigg and Rosner proposed the logarithmic potential |105j in 1979: 

V(r) = Cln(r/ro), (1.18) 

which was motivated by the experimental observation that radial excitations of the S 
ground state are nearly equal for charmonium (589 MeV) and bottonium (563 MeV). The 
independence on the quark mass can be satisfied with the logarithmic potential. 

Richardson proposed in 1978 |11U] to construct the potential in momentum space in 
order to take into account the momentum dependence of a^. The proposed form of the 
potential was: 

V{q') = -l^^^^ll' X ln(l + qVA'r\ (1.19) 

where nj is the number of fermion flavors with mass below q. An advantage of this 
potential is that only one parameter Aqcd is engaged. It reproduces 1/r dependence at 
small distances and yields linear behavior in r at large distances. 

Additional potentials were proposed by BuchmuUer and Tye fST] , Martin [HH] and Cel- 
master and Georgi [51]. It should be stressed, however, that practically all the potentials 
are identical in the region between 0.2 fm and 1 fm, which is the physical radius of the 
charmonium bound states. An additional remark here in place is that calculations of a 
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static spin-independent quark-antiquark potential from Lattice QCD are fitted perfectly 
with the Cornell potential (Eq. 1.17) with parameters k = 0.322 and 1/a^ = 0.54GeV//m 



However the models with spin-independent potentials were not accurate enough to 
predict mass differences between states in the same orbital angular momentum or spin 
multiplets. Additional spin and orbital momentum dependent terms were introduced into 
the Hamiltonian to describe charmonium spectra in more details. 

The development of spin-dependent potential started with the work of Eichten et al. 
[55], Pumplin et al. |104j and Schnitzer |112j . A systematic investigation of spin-dependent 
forces in heavy quark systems was performed in 1979 by Eichten and Feinberg [51]. 

The spin-dependent part of the potential is separated into three terms: 

• Spin-orbit term which splits the states with the same orbital momentum depending 
on the < L ■ S > value (fine structure): 

=< ^ ■ 5 > (3^ - ^)/(2mH (1-20) 

where Vs and Vy are the scalar and vector components of the non-relativistic poten- 
tial, respectively. The expectation value of the L ■ S operator is equal to: 

<L-S>= ^[J(J + 1) - L{L + 1) - SiS + 1)]. (1.21) 

• Spin-spin term describes the interaction between spins of the quarks. It is responsi- 
ble for splitting between the spin singlet and spin triplet states (hyperfine structure): 



2< S,-S2> ^,, 
?>ml 



"^ss = 0I2 VVv-(r), (1.22) 

where the expectation value: 

<Sr-S2>=\[S{S+l)-^-]. (1.23) 

Tensor term contains the tensor effect of the vector potential: 



< Ti2> ,ldVv _ (PVv 
2m^ r dr dr 

where 



= (1.24) 



Tv2 = \\i(S-f)(S-f)-S\ (1.25) 

with the expectation value 

-<L-S>^ -1/2 <L-S> +1/3 < X > 

(2L-f3)(2L- 1) ^ ' 
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The tensor and spin-orbit interactions are responsible for mass splitting of the Xcj{^^Pj), 
J = 0, 1, 2 states. The spin-spin interaction, which splits the vector and pseudoscalar states, 



is responsible for mass differences between J/^p and rjc and between ijj' and rj'^ (see Fig. 1.7 



and Table 1.1). 

Additional important point is that spin-dependent terms in the Hamiltonian are more 
singular than and therefore are illegal operators in the Schrodinger equation. It was 
shown in [M] that the relativistic potential V{p, r) differs from its non-relativistic limit 
in that the coordinate r (which in the non-relativistic limit is the relative distance ri2 = 
ri — becomes smeared over the distance of the order of inverse quark mass. The smearing 
of the potential has the consequence of taming all the singularities in the spin-dependent 
potential, so that the contributing terms become legal operators. The usually implemented 
form of the smearing function is: 



with 



2 2 



pij[,r - r) 



^ -4(r'-r)2 
^3/2 e 



1 1 




4" 




2mimj 


2^2 


{rrii + mj)2 




rrii + rrij 



1.27) 



:i.28) 



where cxo = 1.8 GeV and s = 1.55. 



1.2.2 QCD sum rules 

One of the first calculational techniques successfully applied for mass predictions of the 
charmonium states is QCD sum rule. In general sum rule technique is applicable when 
sums or integrals over observables can be related to normalisation conditions reflecting 
unitarity etc. or to a quantity that can be calculated in the underlying theory. QCD sum 
rule introduced in 1979 by Shifman, Vainstein and Zakharov |115j allows to express low 
energy hadronic quantities through basic QCD parameters. QCD vacuum is typified by 
large fluctuating fields whose strength is characterised by the quark and gluon condensates 
< 0|gg|0 >, < 0\G[^uG^°'^'^''\0 >. These quantities parameterize non-perturbative dynamics 
and vanish in the perturbation theory. When a pair of quarks is injected into the vacuum, 
the subsequent formation of a hadron can be calculated via dispersion relations. The 
heavy quark-antiquark bound states, charmonium and bottonium, appear to be sensitive 
only to the gluon condensate. The success of QCD sum rule technique proved in the 
correct prediction of the mass of charmonium ground state rjc in 1978 |114] and the IP 
charmonium states in 1981 |107] . However no other charmonium states were calculated 
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with this technique due to large power corrections in the expansion of vacuum polarization 
operator. 

1.2.3 pQCD 

The perturbative calculations in QCD (pQCD) are valid for small running coupling con- 
stant, i.e. as it follows from asymptotic freedom in QCD, at large momenta. This technique 
has been applied for charmonium annihilation calculations under the assumption that an- 
nihilation is the process occurring at short distances. Also this technique was applied for 
the calculation of strong radiative corrections to the electromagnetic, radiative (decay to 
77) and hadronic decays of the S and P states in quarkonia by Barbieri et al. How- 
ever the pQCD predictions met with several problems. The first one is that in contrast 
to positronium in charmonium m{qq) 7^ 2m(g), as a result there is an ambiguity whether 
to use m{q) or m{qq)/2 in calculations, which affects the calculated hadronic decay width. 
The second problem is in that the coupling constant is large enough, as ~ 0.3, for the 
charmed quarks so that the lowest order gluon radiative corrections can reach up to 100 
%. However radiative corrections calculated with pQCD are in some cases in good agree- 
ment with experiment, for example the radiative width of rjc is calculated with precision 
of 10% [25]. 

1.2.4 Lattice QCD 

The lattice QCD technique was proposed by Wilson in 1974 |120] . The main intention 
behind lattice QCD was to define an entirely non-perturbative regularisation scheme for 
QCD, based on the principle of local gauge invariance. The four-dimensional space-time 
continuum is replaced by a discrete lattice with the spacing a between lattice sites. Dis- 
tances ~ a or smaller cannot be treated in this way, in the other way stated, the finite 
spacing introduces an ultra-violet cut-off ~ vr/a on momenta in the lattice QCD. In the 
more recent calculations asymmetric lattices have been used, in which lattice spacing for 
time, at, is chosen smaller than for space, a^, with the ratio as/ at ~ 3. The typical lattice 
spacing is ~ 0.07 fm and the overall size of the lattice is 1-4 fm. 

Observables are calculated taking their expectation values in the path integral approach: 
calculated quantities are averaged over all the possible configurations of gauge fields on the 
lattice and weighted with the respective exponent of the action. From the point of view of 
numerical calculations it is convenient to work in Euclidean space-time in which the path 
integral measure can be defined. Most results that are obtained in Euclidean space can be 
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related to the space-like region of the Minkowski world and can, in principle, be analyt- 
ically continued into the time-like region of interest. However, with the results obtained 
on a discrete set of points such a continuation is not straightforward. Fortunately, the 
mass spectrum remains unaffected by the rotation to imaginary time, as long as reflection 
positivity holds, which is the case for most lattice actions used for the charmonium mass 
calculations. 

Past lattice results were often obtained in the "quenched" approximation, i.e. no quark- 
antiquark pairs were allowed to be excited from the QCD vacuum. Recently, the progress 
in computers and computational techniques allowed to do unquenched calculations with 
uu, dd and ss pairs excited from the vacuum. 

1.2.5 Effective theories. NRQCD 

From the point of view of QCD the description of hadrons containing two heavy quarks is a 
rather challenging problem and a proper relativistic treatment of the bound state based on 
the Bethe-Salpeter equation has proved difficult. A non-relativistic treatment of the heavy 
quarkonium dynamics, which is suggested by the large mass of a heavy quark has clear 
advantages. The velocity of quarks in the bound state provides a small parameter in which 
the dynamical scales can be ordered and the QCD amplitudes dynamically expanded. A 
non-relativistic bound state is characterised by at least three scales: the scale of the mass 
m (hard) , the scale of the momentum transfer p ~ mv (soft) and the scale of kinetic energy 
of the quark and anti-quark in the center of mass frame E ~ p^/m ~ mv^ (ultrasoft). In 
a non-relativistic system {v <^ 1) the above scales are hierarchically ordered: m ^ mv ^ 
mv'^. The wide span of energy scales involved makes the lattice QCD calculations extremely 
challenging because of the conflicting requirements that the space-time grid should be large 
in comparison with the largest length of the problem (l/mv^) and the lattice spacing be 
smaller compared with the smallest length 1/m. The hierarchy of scales has the advantage 
that it can be used to substitute QCD with a simpler but equivalent Effective Field Theory 
(EFT). An EFT is the quantum field theory with the following properties: a) it contains 
the relevant degrees of freedom to describe phenomena that occur in a certain limited 
range of energies and momenta and b) it contains an intrinsic scale that sets the limit of 
applicability of the EFT. 

The prototype of an EFT for heavy quarks is the Heavy Quark Effective Theory 
(HQET), which is an EFT of QCD suitable to describe systems with only one heavy 
quark. These systems are characterised by two energy scales: m and Aqcd- HQET is 
obtained by integrating out the scale m and use of an expansion in powers of Aqcd/tti. 
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Following the above discussion, one may conclude that the bound states made of two heavy 
quarks are characterised by more than two scales. Integrating out only the scale m leads 
to an EFT called Nonrelativistic QCD (NRQCD) that still contains the lower scales as 
dynamical degrees of freedom. Starting from NRQCD two approaches can be followed for 
spectrum computation: direct lattice calculations or further integration of the soft scale 
to arrive at an EFT in which only the ultrasoft degrees of freedom remain dynamical, an 
approach called pNRQCD. The lattice calculations with NRQCD and their application in 
charmonium physics are reviewed in |24j . 

1.3 he meson 

The present chapter is devoted to the he singlet state of charmonium (l^Pi). It is the 
most elusive charmonium state with mass below the DD threshold, since its experimental 
observation met with many difficulties. Its quantum numbers make this state not directly 
accessible in e~^e~ annihilations. Also, it cannot be reached in the radiative El transition 
from since the C-parity conservation forbids the transition from 1 to 1^~. Its obser- 
vation at the B-factories in B decay via an intermediate ric{2S) state, which can decay via 
an El transition to h^., is suppressed because of a large hadronic width of the ric{2S). 

In this chapter the main arguments behind the interest in precise measurements of the 
mass and width of he are presented, the theoretical predictions of mass and decay width 
are summarized, and previous attempts of the experimental observation of he are reviewed. 
The results of a Monte- Carlo study of the PANDA detector performance in the foreseen 
he measurements are presented in Chapter |4| 

1.3.1 Spin dependence of the confinement potential 

The main argument besides the interest in precision mass measurement of he is determina- 
tion of the mass difference (hyperfine splitting) between the singlet (^-Pi) and the centroid 
of the triplet P-states (^-Pj). In the framework of potential models the hyperfine splitting 
arises from the spin-spin term of the charmonium potential: 



A 



HF — 



-Vss = - M, 



(1.29) 



where M^^^^ is the mass centroid of the triplet ^Pj states: 



E(2J + 1)M, 
E(2^+l) 



9 



X2 



(1.30) 



1 Physics motivation 



27 



Since the spin-spin potential is a contact one {Vss ~ V VV(r), see Sec. 1.2.1) its 
expectation value is finite only for the wavefunction finite at the origin. Therefore, it 
gives rise to the hyperfine splitting between the triplet and singlet states only in L = 0, 
which is reflected, for example, in the mass difference between J/ip and rjc. The hyperfine 
splitting for the P-wave states should be zero and the higher order corrections should 
provide no more than a few MeV deviation from this result [9], [66]. This is a consequence 
of the assumption that the long-range confinement potential is a pure Lorentz scalar. The 
magnitude of hyperfine splitting larger than a few MeV would be an indication of a vector 
component of the confinement potential, but the experimental observations of he performed 
so far indicate that the latter should not be significant. 



1.3.2 Mass predictions 

Mass predictions for the singlet ^Pi state in charmonium might serve as a validation test 



of the theoretical approaches mentioned above in Sec. |1.2[ Most of the he mass predictions 
are based on the potential model calculations including variation of the strong coupling 
constant and the charmed quark mass rric- The most commonly used form is that 
of the Cornell potential [SS] and the potential proposed by Richardson [llOj . The linear 
confinement term is usually assumed to be scalar but some authors consider an admixture 
of the vector confinement. In some cases relativistic corrections at the level of f^/c^ are 



taken into account. Table 1.2 summarizes all the predictions of mass expressed as the 
hyperfine splitting Ahp = M^^^^ — M^^. It is based on the summaries presented in 
and [HQj. The value {M^^^g)exp = 3525.3 ±0.1 MeV is known fJ] from experiment. 



1.3.3 Decay modes of he 

Theoretical calculations predict he as a narrow state with < 1 MeV total width. The 
possible decay modes of he are electromagnetic or hadronic transitions to lower charmonium 
states or a decay into light hadrons. All predictions for the partial and total widths for the 



he decay are summarized in Table |1.3[ Below some details concerning width calculations 
are discussed. 

The radiative transitions in charmonium occur, as in atoms, from the higher to lower 
mass states with the emission of a photon. The electromagnetic transitions are divided into 
electric dipole transitions (El) and magnetic dipole transitions (Ml). The El transitions 
obey the AL = ±1, AS = selection rule, for the Ml transitions AL = 0, AS* = ±1. For 
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Author 


Year 


Ahf, MeV 


Approach 


Eichten [SI] 


1979 





Potential model with Cornell potential and 
scalar confinement term. 


Ono |l()()j 


1982 


+9.8 


Potential model with smeared hyperfine in- 
teraction. 


McClary [M] 


1983 


-5.2 


Potential model with smeared hyperfine in- 
teraction and relativistic corrections. 


Moxhay [93] 


1983 





Potential model with long-range longitudinal 
color electric field. 


Grotch [69] 


1984 


+ (6-17) 


Potential model with Lorentz scalar and vec- 
tor confinement terms. Predictions depend 
on the fraction of scalar versus vector contri- 
bution. 


Godfrey [61] 


1985 


+8 


Potential model with smeared short-range 
hyperfine interaction. 


Pantaleone jlUlj 


1986 


-3.6 


pQCD calculation (a^ = 0.24). 


Gupta [70] 


1986 


-2 


Potential model with 1-loop QCD correc- 
tions. 


igi m 


1987 


+24.1 


Potential model with short-distance 2-loop 
QCD calculation. 


Pantaleone 


1988 


-1.4 


pQCD calculation (a^ = 0.33). 


Lichtenberg [85] 


1992 


+4 


Potential model with pQCD corrections at 
short distances. 


Halzen [72] 


1992 


-0.7 ± 0.2 


pQCD calculation (a^ = 0.28). 


Bali [21] 


1997 


-1.5 ± 2.5 


Lattice QCD 


Manke [HI] 


2000 


+ (1.7-4.0) 


Lattice QCD 


Okamoto [97] 


2002 


-1.5 ± 2.6 


Lattice QCD 



Table 1.2: Summary of the he mass predictions expressed in terms of the predicted hyperfine 
spHtting Ahf = M^^„^ - Mh^. 
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he the El transition to rjc is possible with the width given in non-relativistic hmit by: 

TiK ^Vc + 7) = ^e^ae\E,f\^, (1.31) 

where k = {Mf — M|)/2Mj is the photon momentum, \Eif \ is the transition dipole matrix 
element: 

I'OO 

\E,f\ = / drr^Ri{r)Rf{r) (1.32) 

where Ri{r) [-R/(r)] are the radial parts of the charmonium initial [final] wave functions, 
which are model dependent. Predictions for El decay width have been made by several au- 
thors, most of the predictions relate the partial width for he decay to the already measured 
partial widths for other charmonium states. For example, Bodwin et al. [33] predict the 
he decay width by scaling with photon energy the radiative decay of Xi,2,3 to J/ip. Renard 
|1U8] also predicts a few keV for the Ml transitions to the triplet P states in addition to the 
El width. Barnes et al. ^27j evaluated the radiative width using two models, the relativised 
Godfrey-Isgur model and a nonrelativistic potential model for the wave functions. 

The hadronic transitions of he to J/ip are allowed energetically with the emission of 
one or two 7r°. The J/ipirir width is predicted to be larger than J/ipn^. The predictions 
for J/tpTCTf and J/ipTc^ by Kuang et al. [82] were made in the framework of QCD multipole 
expansion. In the work of Chen and Yi [H] the same approach is used, taking into account 
also the 2-loop corrections. 

The partial width of he decay into light hadrons is actually the width of decay into 
three gluons (ggg) or two gluons plus a gamma (ggj), where gluons hadronize. In most 
works scaling of the partial width from the known decay is used, as for example Bodwin 
et al. [33] use the Xi partial width into hadrons for their prediction. It is important to 
mention here that he decay into light hadrons is OZI suppressed, which is typical for a 
charmonium state (see Sec. 1.1.2). 

The total width prediction is the sum of partial widths. Some authors quote the sum 
of their own predictions adding the results of other authors for the missing channels. 



1.3.4 Experimental observations 

The first attempts for observation of the he date back to 1985 and the R704 experiment 
at the CERN Intersecting Storage Rings (ISR), where the proton-antiproton annihilation 
has been used to produce the charmonium states. The expected high hadronic background 
was an argument to chose the J/ip — >■ e'^e~ decay as a tag for charmonium formation in the 
R704 experiment. To detect high mass e'^e~ pairs R704 was constructed as a nonmagnetic 
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Author 


Year 




T{J/tlJ7r^) 


r(J/^/^7rV) 


T{hadrons) 


T {total) 






(keV) 


(keV) 


(keV) 


(keV) 


(keV) 


Renard [T08j 




1976 


240 






370 


500-1000 


Galkin [59] 


1990 


559 










Novikov [95] 


1977 


975 






60-350 




McClary [HH] 


1983 


483 










Chao [12] 


1992 


385 










Bodwin [33] 


1992 


450 






530 


980 


Casalbuoni 


10] 


1992 


450 










Ko [HI] 


1992 


400 


1.6 








Gupta [H] 


1993 


341.8 










Barnes [27] 


2005 


352-498 










Kuang [82] 


1988 




2 


4.12 


53.7 


394 


Chemtob [13] 


1989 




0.0061 


52.6 






Chen [H] 


1992 




0.29-0.58 


4.1-7.1 


19-51 


360 



Table 1.3: Partial and total width predictions for the he decay. 



spectrometer with two arms positioned symmetrically relative to the beam axis to detect 
lepton pairs. A 2.3 a signal was claimed with the observation of 5 — > J/ip + X events 
[23] near the Xc center of gravity. This signal was interpreted as the h^ — J/ipT^^ decay 
with the he mass 3525.4 ± 0.8 ± QAMeV/c^. 

After shutdown of the ISR, the E760 experiment was proposed at the antiproton source 
complex of Fermilab to continue the study of h^ in proton-antiproton annihilations. An 
antiproton beam (up to 4.0 ■ 10^^ accumulated p) intersecting with an internal hydrogen 
jet target produced instantaneous luminosity in the range (3^9)- 10^° cmr'^s~^ . A system 
of stochastic cooling compensated the effects of scattering and energy loss in the target to 
keep the beam energy spread at the level Ap/p < 2.5 • 10^^. The search for he has been 
performed in the vicinity of the Xc center-of-gravity in small energy steps (< 500 keV). 
The E760 experiment was a non-magnetic spectrometer optimised for the detection of the 
electromagnetic final states with an electromagnetic calorimeter, composed of lead-glass 
scintillators, covering polar angles from 2° to 70° and having full azimuthal coverage. The 
following decay modes of he have been searched for: 

p-hp ^ /ic ^70 7 -> (77)7, (1-33) 
p + p^he^ J/i)Tx^ ^ {e'^e-)TX^, (1.34) 
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3522 3523 3524 3525 3525 3527 3523 3529 3530 



Figure 1.10: Evidence for the he — )• J/i/jtt^ observation in the E760 experiment 



c.o.g. 



indicates the position of the center of gravity (Eq. 1.30) of the Xcj states. 



p + p ^ he ^ J/ip2'K ^ {e^e )27r. 



(1.35) 



According to theoretical predictions (see Table 1.3) the first decay mode should be the 



dominant one, but due to low branching ratio for the transition rjc — )■ 77, the event rates 
for all three decay modes were assumed comparable. The evidence for the resonance in 
the p + p — 7- J/ipTT^ channel has been demonstrated with the 59 observed events 



The observation plot is presented in Fig. 1.10 The mass of the resonance was estimated at 
Mr = 3526.2±0.15±0.2 MeV/c^. An upper limit on the resonance width T < 1.1 MeV was 
set at the 90 % confidence level. No events were found to fit the reaction p+p — )■ J/ipn^Ti^ 
or p + p — )■ J/ip TT^n^. 

The E835 experiment, which was an upgrade of E760, continued the study of he in two 
runs with the collected luminosities of 141.4 pb~^ in Run-I (1997) and 113.2 pb^^ in Run-II 
(2000). The most important modification of the detector with respect to E760 was an 
upgrade of electronics of the calorimeter, which included pulse shaping and installation of 
a time-to-digital converter for each calorimeter channel to reduce pileups and replacement 
of the forward calorimeter with lead-glass modules. The two reactions p + p ^ he ^ rjcj 
and p + p^he—>J/ipn^ were investigated. 

With the 6 times larger collected statistics in comparison with E760 no evidence of 
the he — )■ J/tpTc^ signal was observed. In conclusion the resonance reported by E760 was 
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Figure 1.11: Observation of the he — )• ?7c7 in the E835 experiment [6]. The inset shows the he 
search region in an expanded scale. 



not confirmed. However the evidence of the signal in the he — )■ rje'f mode was observed 



with rje — )■ 77 (1.33) with the statistical significance (P- value) 0.001 [6]. The observation 



plot is presented in Fig. 1.11 The reported mass is M/j = 3525.8 ± 0.2 ± 0.2 MeV and 
the resonance width F < 1 MeV. The signal was observed in both runs with comparable 
strength. 

The most recent observation of he was made by the CLEO experiment |111] at the 
Cornell Electron Storage Ring in the reaction: 

e+e" ^ ^(2^) ^ Tc^he, he Vd, vr° ^ 77. (1.36) 

The study has been performed for both CLEO III and CLEO-c detector configurations, 
which are characterised by 93 % solid angle coverage for both charged and neutral product 
particles. 

Two methods of the analysis have been used: exclusive analysis with the reconstruction 
of rje decay which has higher signal purity and inclusive analysis with higher signal yield. 
Both methods identify he as an enhancement in the spectrum of neutral pions from the 
decay ip{2S) — >■ Tr^he- For the exclusive analysis rje was reconstructed in the following chan- 
nels: KgK^TT^, KIK^tt^, K~^K~7i~^'n'~ , Tr+vr'Tr+Tr", K^K~tt^ and tt^tt't]. The reported 
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Figure 1.12: PDG [4j plot for the he mass averaging. 



[Ill] mass is M{hc) = 3524.4 ± 0.6 ± OAMeV and the significance of the observation is 
greater than 5a. 

The current PDG mass of M{hc) = 3525.93 ± 0.27 MeV 0, is based on the four 
observations mentioned above. The plot demonstrating the mass averaging from PDG [1] 
Fig. [1.12] shows that there are some discrepancies between the observations, which motivate 
improved experimental search with higher statistics and better resolution. 



1.4 Angular distributions of photons in radiative tran- 
sitions in charmonium 

Formation of charmonium as a resonance in the pp annihilation provides an experimen- 
tal tool to determine the electromagnetic radiative multipole structure for the transition 
between the resonance and the selected final state. The study of a radiative angular distri- 
bution provides information on the relative annihilation amplitudes through the possible 
helicity states. In those cases, in which only one multipole and one helicity state is allowed 
by the conservation laws, the angular distribution can be used to verify J^*^ characteristics 
of the newly discovered states, such as he. The two independent formalisms can be applied 
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for the description of radiative angular distributions, helicity amplitudes and transition 
multipoles. 

1.4.1 Helicity formalism 

Here a short summary of the helicity formalism is provided, based on |90] . The formalism of 
helicity amplitudes was proposed by Jacob and Wick in 1959 [7H] to describe the collisions 
of particles with spin. In this formalism states are labeled with the projection of the total 
angular momentum along the direction of particle momentum, which avoids relativistic 
complications resulting from division of the total angular momentum into spin and orbital 
parts. Additional useful property of the helicity formalism is that for the two-particle 
states, in the center-of-mass frame, the orbital angular momentum is oriented perpendicular 
to the direction of relative motion of the two particles. With the helicity defined as: 

n = (1.37) 

the component of the spin along the direction of relative motion will be the component 
of the total angular momentum along this direction, since L + S and L ■ f = 0. 

Moreover, the helicity quantum number is invariant under the ordinary rotation, however 
if the momentum vector is rotated by 180° the helicity quantum number changes sign. 

1.4.2 Prediction for the angular distribution of 'j in pp ^ rjc'j 

The angular distribution of 71 in the pp ^ he ^ liVc lilli based on the helicity 
formalism, is provided in [99j|. The joint angular distribution for this process can be 
expressed as: 

A ^l,^J.'=±l 



with the angles and helicities as indicated in Fig. 1.13 Here the angles 6', (j)' specify the 
Vc 77 decay in the r]^ rest frame. W{6; 9',(f)') does not depend on the angles 9' and 0' 
since the pseudoscalar state rjc decays into 77 isotropically. 
A normalised joint angular distribution is defined as: 

Bl + 2BI 

A constant R is introduced defining the fractional contribution of the helicity-one am- 
phtude into the initial production process: 

2R? , , 



#(g;g^0O = g ;r, o'^2 • (1-39) 
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Figure 1.13: Angles and hclicities involved in the process pp ^ he liVc 7i77; 7i is the 
photon emitted in the he decay. 



For each transition the joint distribution W{9; 9' , (f)') can be expressed in terms of the 
observables {Ki} and the elementary trigonometric functions. {Ki} can be expressed in 
terms of the helicity amplitudes Ai, which in turn can be written in terms of the multipole 
transition amplitudes a^. The general relation between the two sets of amplitudes is: 

= E «'^[|^]'^'(^' 1; 1' ^ - Mf, (1-41) 

where <; | > is the Clebsch-Gordan coefficient. 

The he — > rjc'fi decay proceeds via the El amphtude only, so that ai — Aq — 1. The 
angular distribution is: 

—w{e;e',<f>') = {K, + K2cos^e), (1.42) 



with 



and 



so that 



with 



K, = l- ^R, (1.43) 
K2^^R-1, (1.44) 
W(e) ^W{7r/2)(l + acos^e), (1.45) 



3^-2 , , , 

a = (1.46) 
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Because of C-parity conservation the Bi hehcity state does not enter into he production 
and R — 0. Therefore, the angular distribution of 7-rays in the he decay has particularly 
simple form: 

W{e) ^W{7r/2)sm'^e. (1.47) 



Chapter 2 

The PANDA experiment 



The PANDA (Pbar ANnihilations at DArmstadt) is a hadron physics experiment which is 
planned as part of the future Facihty for Antiproton and Ion Research (FAIR) at Darm- 
stadt. The cooled antiproton beam with a momentum between 1.5 GeV/c and 15 GeV/c 
will be provided for the PANDA experiment in High Energy Storage Ring (HESR). The 
antiprotons will interact with an internal hydrogen target to reach a peak luminosity up 
to 2 ■ 10^2 cm-^s-^ 



2.1 The PANDA physics program 

Experimental program of the PANDA is concentrated on the studies of strong interaction 
and hadron structure performed with the aid of an antiproton beam. The major topics of 
the PANDA physics program are summarized below. 



2.1.1 Charmonium spectroscopy 

Precise measurements of the properties (energies, total and partial widths) of the charmo- 
nium states below the DD threshold have fundamental importance for understanding of 
QCD. The theoretical predictions for the charmonium spectrum have been performed in 



several theoretical approaches outlined in Sec. L2 starting with the simple non-relativistic 
potential model and ending with the latest calculations using the effective field theoretical 
approach (NRQCD). Most of them can be verified with high precision with the experi- 
mental results provided by PANDA. Moreover, the recent discoveries of new states (X, Y, 
Z) at the B factories renewed interest in physics of charmonium and stimulated intense 
experimental and theoretical activities. All charmonium states can be formed directly in 
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pp annihilation, in particular those with quantum numbers other than J^^ = 1 and 
their masses determined with precision of several 100 keV, which is 10 to 100 times better 
than achieved in the e+e~ experiments. 



2.1.2 Search for gluonic excitations (hybrids and glueballs) 

Glueballs and hybrids are the bound states predicted by QCD in which the gluonic degree 
of freedom explicitly manifests itself. Glueballs are the states of pure glue. Prediction of 
their mass spectrum was one of the early achievements of the lattice QCD. There is about 
15 glueballs in the mass range accessible with the PANDA experiment [92]. Glueballs 
with exotic quantum numbers are called oddballs; they cannot mix with the "normal" 
mesons. As a consequence, they are predicted to be rather narrow and easy to identify 
experimentally. The hybrids consist of a qq pair and an excited glue. The additional 
gluonic degree of freedom manifests itself as a contribution to quantum numbers. In the 
simplest scenario this corresponds to adding the quantum numbers of a gluon to a simple 
qq pair. Charmonium hybrids are expected to exist in the 3-5 GeV mass region. The lowest 
lying charmonium hybrid is expected with J^^ = and the mass of 4.34 GeV [25] . 



2.1.3 Study of hadrons in nuclear matter 

The change of hadron masses inside nuclear medium has been proposed as the manifestation 
of partial restoration of the chiral symmetry breaking predicted by QCD due to finite 
density, and thus to be an indicator of changes in quark condensates. In particular, the 
attractive mass shifts reflecting reduced quark pair condensate at finite density have been 
predicted for vector mesons. Numerous experiments have been devoted to deduce in- 
medium mass shifts or nuclear potentials for hadrons in the light quark sector {u,d,s), 
both for the vector and for the pseudoscalar mesons. The energy range of HESR and 
detection capability of PANDA will permit to extend this type of studies into the sector 
of charmed hadrons for both open and hidden charm. 

Another open question where PANDA can contribute is the J/ipN dissociation cross 
section on a series of nuclear targets. Knowledge of this parameter is important for inter- 
pretation of the J/ip suppression effect as the signature of quark-gluon plasma formation 
in ultra-relativistic heavy ion collisions. 
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2.1.4 Open charm spectroscopy 

The D meson, consisting of the heavy and hght constituents, is a very interesting object 
since it combines the aspect of the heavy quark as a static color source on the one hand, 
and the aspect of chiral symmetry breaking and restoration due to the presence of the 
hght quark on the other hand. The phenomenological quark model was able to describe 
the excitation spectra of heavy-light systems for the lowest states. However the series of 
new observations in the open charm spectrum starting with the unexpected discovery of a 
narrow Ds(2317) state by BaBar [18] attracted much interest since the new states do not 
fit well into the quark model predictions. Suggestions appeared that the scalar and the 
axial-vector Dg mesons have an exotic non-cs structure, with possible tetraquark |118j or 
DK molecule interpretation [21] ■ Precise measurements of the widths and decay modes 
will allow to distinguish between different interpretations. 



2.1.5 Hypernuclear physics 

Hypernuclei are the systems in which one or more protons or neutrons are replaced by the 
hyperons. Although single and double A-hypernuclei were discovered decades ago, only 6 
double A-hypernuclei have been discovered up to now. The p beam at FAIR will allow 
efficient production of hypenuclei with more than one strange barion. These studies are 
expected to provide valuable information on hyperon-hyperon interaction. At PANDA not 
only S~ but also Q~ atoms can be studied for the first time supplying unique information 
on hyperon-nucleus interaction. 



2.1.6 Electromagnetic Processes 

The experimental setup of PANDA will offer the opportunity to study also a certain class 
of hard exclusive processes such as Wide Angle Compton Scattering (WACS) and determi- 
nation of electromagnetic formfactors in the time-like region from the process pp — > e+e~ 
and/or . These processes will be accompanied with a significant hadronic background, 
therefore the particle identification capability of PANDA will be essential to execute them 
successfully to complement the excellent performance of the Electromagnetic Calorimeter 
used to detect lepton pairs. 
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2.2 FAIR and HESR storage ring 



The future FAIR facility will be built profiting from the existing GSI facility. Its layout is 
presented in Fig. 2A FAIR will provide good quality antiproton beam whose production 
and delivery chain is briefly described below. The proton linac will inject 70 mA of protons 
at 70 MeV into the existing heavy ion synchrotron SIS18. The SIS18 will be used to boost 
the proton energy to 2 GeV for injection into the superconducting synchrotron SISIOO. The 
SISlOO in turn is employed to accelerate protons to the envisaged antiproton production 
energy of 29 GeV. Up to 4 ■ 10^'^ protons per SISlOO-cycle can be accelerated and before 
fast ejection to the antiproton production target the proton beam circulating in SISIOO has 
to be compressed to a single bunch of approximately 7.5 m length ( 25 ns duration). With 
a bunch containing 4 ■ 10^^ primary protons at 29 GeV approximately 2 ■ 10^ antiprotons 
are expected from the antiproton production target within the phase space acceptance of 
the following separator and the Collector Ring (CR). It has been demonstrated at the 
antiproton-facilities at CERN and FNAL, that collection and pre-cooling of antiproton 
beams on the one hand and accumulation on the other hand are optimized by means of 
two differently designed rings serving different purposes. In the CR antiprotons with a 
mean momentum of 3.8 GeV/c are stochastically cooled and rebunched to be transferred 
to RESR. In RESR accumulation of antiprotons takes place. There is stochastic cooling 
applied during the whole accumulation process so that up to 7 ■ 10^° antiprotons per hour 
should be obtained. Finally the antiproton beam is rebunched into one single bunch and 
transferred to the HESR. 



The synchrotron HESR (Fig. 2.2) accelerates or decelerates antiprotons to the desired 
momentum and serves as a storage ring for the internal target experiment PANDA. Stochas- 
tic and electron cooling of the antiproton beam are foreseen to satisfy requirements of the 
planned experiments. HESR is designed as a racetrack shaped ring with a circumference 
of 574 m, including the 132 m long straight sections. Within one of these sections the 
PANDA experiment will be installed. The opposite side is nearly completely used by the 
electron cooler. At both ends of the electron cooler compensation solenoids are foreseen. 
The equipment for stochastic cooling is installed in both straight sections. Two operation 
modes of the HESR are foreseen: high-luminosity and high-resolution with the parameters 



presented in the Table 2.1 
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Figure 2.1: Layout of the future FAIR facility at GSI. 




Figure 2.2: Schematic view of the HESR with the indicated position of antiproton injection from 
RESR, coohng devices and experimental installations. Direct injection of protons from SIS18, 
which is useful for tuning of PANDA is foreseen in the later stage of construction of the system. 
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High resolution (HR) 


Luminosity of 2 ■ 10^^cm~^s~^ for 10^" p, 
RMS momentum spread (Jp/p < 4 ■ 10~^, 
1.5 to 9 GeV/c, electron cooling up to 9 GeV/c 


High luminosity (HL) 


Luminosity of 2 ■ 10^^cm~^s~^ for 10^^ p 

RMS momentum spread <7p/p ~ 10"'^, 

1.5 to 15 GeV/c, stochastic cooling above 1.5 GeV/c 



Table 2.1: Parameters of the HESR operation modes. 



2.3 PANDA detector overview 

The main design requirements for the PANDA detector, which are determined by the 
physics program, are: An acceptance, high resolution for tracking and calorimetry, particle 
identification and high rate capability. The detector is divided into two parts, the Target 
Spectrometer (TS) with the superconducting solenoid magnet and the Forward Spectrom- 
eter (FS) based on the dipole magnet. The main components of the PANDA detector 



(Fig. 2.3) are briefly described below according to the PANDA Technical Progress Report 



[50] and the PANDA Physics Performance Report [52]. 
2.3.1 Target Spectrometer 

The target spectrometer, surrounding the interaction point, is split into the barrel part 
which covers polar angles above 22° and the two end-caps, forward and backward. The 
forward end-cap is covering angles down to 5° in the vertical and 10° in the horizontal plane. 
The main design requirement to the target spectrometer is its compactness to minimize 
the cost of the magnet; this requirement imposes limitations on the size of the inner layers 



of the detector. The target spectrometer (Fig. 2.4) will contain the micro vertex detector 
(MVD), the tracking device with two possible alternatives presently considered, straw tube 
tracker (STT) or time projection chamber (TPC), forward GEM detectors, the PID de- 
tectors (Barrel DIRC, Forward Endcap DIRC and Barrel Time-of-Flight), electromagnetic 
calorimeter (EMC) and muon detectors. All the subsystems of the target spectrometer are 



described below except for EMC, which will be covered in a separate section (Sec. 2.4). 
2.3.1.1 Target 



The requirements to the target are dictated by two factors, the limited space available due 
to the desired minimal distance from the interaction point to the vertex tracker and the 
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Figure 2.3: Artistic view of the PANDA detector. 




IS 



Figure 2.4: Section of the target spectrometer showing the elements hsted in Sec. 



2.3.1 
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need to provide about 4 • 10^^ hydrogen atoms per cm^ in order to achieve the designed 
luminosity of 2 • 10^^ cm~'^s~^ with 10^^ antiprotons stored in HESR. To satisfy these 
requirements two complementary techniques have been proposed: the cluster-jet target 
and the pellet target. 

The cluster-jet target is a homogeneous, variable density stream of hydrogen clusters 
formed by an expansion of pressurized cold hydrogen gas into vacuum through a Laval-type 
nozzle. This leads to condensation of hydrogen molecules into a narrow jet of hydrogen 
clusters with the size ranging from 10'^ to 10^ molecules. However, the size of the interaction 
zone determined by the lateral spread of clusters is estimated at 15 mm at present, 
which exceeds the required. The desired density has not yet been reached with the 8 ■ 10^^ 
atoms/cm^ achieved with the target prototype, therefore R&D work of the participating 
teams is continuing to reach the goal. 

The pellet target provides a stream of frozen hydrogen micro-spheres, called pellets, 
which traverse the interaction point at the velocity of about 60 m/s with the typical pellet 
size of 25 - 40 fim. Such target type is used presently at the WASA detector at the COSY 
machine in FZ-Jiilich. It provides the density of 5 ■ 10^^ hydrogen atoms per cm^ required 
by PANDA, however the pellet stream is far from being uniform with the average distance 
between pellets of 2 mm. Anyway, the pellet target is an alternative of inferior quality 
since it provides reaction products in the pulse regime. Every time a pellet crosses the 
circulating antiproton beam there is a pulse of products. 

2.3.1.2 Solenoid Magnet 

The magnetic field of 2 T will be generated by a superconducting solenoid coil with the 
inner radius of 90 cm and the length of 2.8 m. The field homogeneity, which is a critical 
parameter for the tracking devices, is foreseen to be better than 2 % over the volume of 
the vertex detector and the central tracker. 

2.3.1.3 Microvertex Detector 

The purpose of the microvertex detector (MVD) within PANDA is to detect secondary 
vertices from the D meson and hyperon decay and to improve the transverse momentum 
resolution of tracking. MVD will consist of the barrel part and eight disks perpendicular 
to the beam axis in the forward direction. The barrel part will consist of four layers, with 
the two inner layers of radiation hard silicon pixel detectors and the two outer layers of 
double sided silicon strip detectors. An inner radius of the barrel MVD is 2.5 cm and an 
outer radius is 13 cm. The first two of the forward disks of MVD are made entirely of pixel 
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detectors, the following four are a combination of strip detectors on the outer radius and 
pixel detectors closer to the beam pipe. The last two disks, made entirely of silicon strip 
detectors, are placed further downstream to achieve a better acceptance for the hyperon 
cascades. The readout of pixel detectors via bump-bonded wafers with ASICs is foreseen 
as a default solution, while the pixel readout chip based on the 0.13 /xm CMOS technology 
is under development for PANDA. 

2.3.1.4 Central Tracker 

The central tracker for PANDA will consist of two subsystems. For the barrel part around 
MVD two possible options are considered at present, Straw Tube Tracker (STT) or Time 
Projection Chamber (TPC). The second subsystem consists of three sets of GEM trackers 
covering forward angles. The main requirements addressed to the tracking devices are a 
good momentum resolution at the percent level and the ability to handle high count rates, 
which correspond to the designed luminosity of 2 • 10^^ cm~^s""^. 

The STT consists of the aluminised mylar tubes, called straws, 10 mm in diameter. 
The straws will be filled with a gas mixture of Argon and CO2 used as a quencher at 1 bar 
overpressure. They are arranged in 24 planar layers mounted in a hexagonal shape with 
4200 straws in total. The resolution in x and y coordinates of about 150 /im is expected. 
The 8 central layers of the straws are tilted to achieve the acceptable resolution of 3 mm 
also in z. A material budget is estimated at 1.3% of the radiation length. 

The main advantages of TPC as an alternative to STT are low material budget and 
an additional particle identification capability through energy loss measurements. The 
TPC consists of two large gas-filled half-cyhnders enclosing the target and beam pipe and 
surrounding the MVD. An electric field along the cylinder axis separates positive gas ions 
from electrons created by ionizing particles traversing the gas volume. The electrons drift 
with constant velocity towards the anode at the upstream end face and create avalanches 
detected at the pad readout plane yielding information on two coordinates. The third coor- 
dinate of the track comes from the measurement of the drift time of each primary electron 
cluster. Traditionally in TPCs an amplification with multi-wire proportional chambers 
gated with an external trigger is used to avoid a continuous backflow of positive ions into 
the drift volume. The high counting rate expected with PANDA does not allow for ap- 
plication of the external trigger and a novel readout scheme was proposed based on GEM 
foils as amplification stages. Additional challenge for the TPC operation comes from the 
slow drift, therefore the necessity to deal with a large number of tracks coming from dif- 
ferent interaction events at the same time. For the assignment of a track to the specific 
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interaction the tracklet reconstruction has to take place close to the readout electronics 
with the following matching by time correlations with other detectors. 

To perform the tracking at angles below 22° in the forward direction, the gaseous 
micropattern detectors based on GEM foils as amplification stages are proposed. These 
detectors have rate capabilities three orders of magnitude higher than drift chambers and 
therefore can handle the particle flux of 3-10^cm~^s~^ expected in the vicinity of the 5 cm 
diameter beam pipe. Current design of PANDA foresees three double planes of GEM 
detectors placed at 1.1, 1.4 and 1.9m downstream of the target. The readout is foreseen 
with the same front-end chips as are planned to be used for the silicon microstrips. 

2.3.1.5 Pcirticle Identification 

Several systems are foreseen for the charged particle identification over a wide range of 
energies and momenta to fulfil the physics objectives of PANDA. The main part of the 
momentum spectrum above 1 GeV/c will be covered by Cherenkov detectors. For the low 
momenta the energy loss measurements can be performed with tracking detectors. In 
addition a time-of-fiight barrel can be used for the identification of slow particles. 

The Cherenkov detectors are based on the principle that a charged particle propagating 
in the medium with an index of refraction n, with velocity f^c > 1/n, emits a cone of light 
with an opening angle Qc — arccos(l/n^) to the direction of its propagation. The velocity 
information determined from Qc, in combination with the momentum information from 
the tracking detectors determines the mass of the particle. For the polar angles between 
22° and 140° the barrel DIRC (Detection of Internally Refiected Cherenkov light) detector 
based on fused silica (synthetic quartz) with a refractive index of 1.47 will be used. The 
1.7 cm thick perfectly polished quartz slabs will surround the beam line at a radial distance 
of 45 - 54 cm. The hght emerging from each slab will be focused with lenses onto micro- 
channel plate photomultipher tubes (MCP PMTs) which are insensitive to magnetic field. 

The polar angles between 5° and 22° will be covered with the Forward Endcap disk 
DIRC based on a similar concept. The fused silica radiator will be arranged in shape of 
a disk with the focussing done by mirroring quartz elements reflecting onto MCP PMTs. 
The disk will be 2 cm thick and will have an outer radius of 110 cm. 

For slow particles at large polar angles particle identification will be provided by a 
time-of-flight detector. In the target spectrometer the flight path is only of the order of 50 
- 100 cm. Therefore the detector must have a very good time resolution between 50 and 
100 ps. As detector candidates scintillator bars and strips or pads of multi-gap resistive 
plate chambers (RPC) are considered. 
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2.3.1.6 Muon Detectors 

The range tracking system will be implemented in the yoke of the solenoid magnet for the 
proper separation of primary muons from pions and decay muons. A fine segmentation 
of the solenoid yoke, playing the role of an absorber with aluminium drift tubes used 



as interleaved tracking detectors is foreseen (see Fig. 2.4). In the barrel region the yoke 
is segmented in such way that a first layer of 6 cm iron is followed by 12 layers of 3 cm 
thickness with 3 cm wide gaps for the detectors. In the forward end-cap more material is 
needed, therefore, six detection layers will be placed around five iron layers each 6 cm thick 
within the door, and a removable muon filter with the additional five layers of 6 cm iron is 
installed in the space between the solenoid and the dipole. 



2.3.2 Forward Spectrometer 



The Forward Spectrometer (FS) (see Fig. 2.3), based on a wide-aperture dipole magnet, 
consists of the following subsystems: forward tracker (drift chambers or straw tubes), 
particle identification, represented by the RICH detector and time-of-fiight wall, forward 
electromagnetic calorimeter and forward muon detector. 



2.3.2.1 Dipole Magnet 

The purpose of the dipole magnet with a maximum bending power of 2 Tm is momentum 
analysis of charged particles emitted in the forward direction. In the current design, the 
magnet yoke will occupy about 2.5 m in the beam direction starting at 3.5 m downstream 
of the target. Thus, it covers the entire angular acceptance of the forward spectrometer of 
±10° in the horizontal and ±5° and in the vertical direction, respectively. 



2.3.2.2 Forward Trackers 

The forward tracker, consisting of the set of six wire chambers (two located in front, two 
within and two behind the dipole magnet), will define the trajectories of charged particles 
deflected with the dipole magnet. The drift chambers with small-size cells and straw tubes 
are considered as possible alternatives for the forward tracker. Each set will contain three 
pairs of detection planes, one pair with vertical wires and two pairs with wires inclined by 
+10° and -10° to reconstruct tracks in each chamber separately. The expected momentum 
resolution of the system for 3GeV/c protons is 6p/p = 0.2%. 
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2.3.2.3 Forward Particle Identification 

The RICH detector is proposed for tlie t^/K and -ft'/p separation at tlie liigliest momenta 
with a dual radiator design similar to the one used at HERMES. Using two radiators, 
silica aerogel and C4F10 gas, makes it possible to achieve the vr/fC/p separation in a broad 
momentum range 2-15 GeV/c. A lightweight mirror focusing the Cherenkov light on an 
array of photo tubes placed outside the active volume will be used to reduce the thickness 
of the detector. 

For the charged particle identification in the forward direction at momenta below 4- 
5 GeV/c the time-of-fiight technique can be applied. Time-of-fiight wall made of plastic 
scintillators, read out on both ends by fast photo tubes, will be placed at 7 m from the 
target. In addition to the main TOF wall, two side walls will be installed inside the 
dipole magnet gap, near the left and right wall of the magnet yoke, in order to register 
low momentum charged particles bent in the dipole field towards the yoke. The relative 
time-of-fiight between two charged tracks will be measured for particles reaching any of 
the time-of-fiight detectors including the barrel TOF. With the expected time resolution 
of 0" = 50ps Tr/i^ and -ft'/p separation on a 3 a level will be possible up to momenta of 
2.8GeV/c and 4.7GeV/c, respectively. 

2.3.2.4 Forward Electromagnetic Calorimeter 

Detection of photons and electrons in the forward direction will be performed with the 
Shashlyk-type calorimeter. The modules of Shashlyk calorimeter form a sandwich of in- 
tercalating plates of lead and plastic scintillator. The readout of the module is performed 
with wave-length shifting fibres passing through the block and coupled to photomultipliers. 
The detector will consists of 26 rows and 54 columns, 1404 modules in total, with the size 
of front face of the module 55 x 55 mm. With the Shashlyk modules constructed for the 
KOPIO experiment at BNL |16] the energy resolution of 4 %/ \/E has been achieved. 

2.3.2.5 Forward Muon Detectors 

The range tracking system for discrimination of pions from muons in the very forward 
direction similar to the muon system of the target spectrometer has been designed. 

2.3.3 Data Acquisition 

The high data rates of 2 ■ 10''events/s expected at PANDA and the complexity of the 
experiment required the data acquisition concept different from the traditional two layer 
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hierarchical approach. Traditionally a subset of specially instrumented detectors is used 
to evaluate a first level trigger condition. For the accepted events, the full information 
of all detectors is then transported to the next higher trigger level or to storage. The 
data acquisition (DAQ) concept of the PANDA assumes that each subdetector will be self- 
triggering, i.e. signals are detected autonomously by the sub-systems and are preprocessed. 
This requires hit-detection, noise-suppression and clustering at the readout level. The data 
related to a particle hit, with a substantially reduced rate in the preprocessing step, are 
marked with a precise time stamp and buffered for further processing. The trigger selection 
finally occurs in computing nodes which access the buffers via a high-bandwidth network 
fabric thus each of sub-detectors can contribute to the trigger decision on the same level. 

2.4 Electromagnetic calorimeter 

The detection of photons and electrons in a wide range of energy from several MeV up to 
15 GeV is required in the PANDA experiment within the full sohd angle. This task will 
be performed by an electromagnetic calorimeter (EMC). 

Calorimeters, in contrast to tracking detectors, degrade the initial energy of the inci- 
dent particle sharing it among the shower products, which are measured to determine the 
primary information on the particle. The important feature of calorimeters is that they are 
sensitive not only to charged particles but also to neutral particles. The electromagnetic 
shower propagation proceeds in the way that e.g. the initial electron creates a photon 
via bremsstrahlung process and in turn the photon provokes pair production, an electron 
and a positron. The shower develops resulting in increasing number of particles until the 
average energy of products decreases below a critical energy at which the main process of 
energy loss is ionisation. The longitudinal size of the shower is determined by the radiation 
length Xq of the material where the shower propagates. The transverse dimension of the 
shower is characterised by a Moliere radius, which, by definition, is the radius of a cylinder 
containing on average 90 % of the shower's energy deposition. 

Two types of electromagnetic calorimeters are used: homogeneous and heterogeneous 
(or sampling). In the homogeneous one the same medium is used to cause the shower 
development and to detect the produced particles, whereas sampling calorimeters have 
sandwich structure which consists of passive layers of converter with high Z alternated 
with active detection layers. 

There are several general requirements, which determine the design of the electromag- 
netic calorimeter in PANDA. The structure of EMC is determined by the design of the 
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PANDA Target Spectrometer, which assumes that EMC is placed inside the yoke of the 
solenoid. From the point of view of cost the size of the yoke is limited which in turns 
requires the calorimeter to be built from a material with small radiation length. On the 
other hand operation of EMC in strong magnetic field of the solenoid makes it impossible 
to use conventional PMTs for light readout. The fast response of EMC is necessary because 
of high event rates up to 2 ■ 10"^ events/s expected with PANDA. The requirements of good 
energy resolution over a wide energy range and high granularity are determined by the 
physics program of PANDA, which requires reconstruction of electrons and photons in the 
energy range from the low energy threshold at 10 MeV up to 15 GeV. Last but not least is 
the requirement of radiation hardness of the used components. More specific design issues 
are determined from the studies of possible registration in the particular physics channels. 

Combination of these requirements led to a design of the homogeneous electromagnetic 
calorimeter for the PANDA detector with PWO-II scintillator material and readout with 
Large Area Avalanche Photodiodes (LAAPD) and vacuum phototriodes. It will be shown 
below how the proposed design fulfils the requirements. 

The requirement of compactness determines to a significant extent the choice of PANDA 
EMC as a homogeneous calorimeter because a typical lateral size of the sampling calorime- 
ter modules needed to cover the maximum energy of photons is much larger than the lateral 
size of the crystal for a homogeneous calorimeter. Additional arguments in favor of the 
homogeneous calorimeter are: typically worse energy resolution of sampling calorimeters 
insufficient for the study of many physics channels and an elevated value of the low energy 
threshold for the detection of photons. 

2.4.1 Scintillator material (PWO) 

The choice of lead tungstate, PbWO^ (PWO), as a scintillator material for PANDA EMC 
is motivated by the requirements of compactness, fast response and radiation hardness. Its 
density p=8.28 g/cm? and radiation length, Xq = 0.89 cm, ensure that with this choice 
the most compact design will be achieved as compared with any other scintillator used 
in high energy physics. By its physical properties PhWO^^ is a birefringent, tetragonal, 
scheelite-type crystal belonging to the space group of monoclinic raspite, colourless and 
transparent in the visible spectrum. Its scintillation mechanism, optical properties and 
light yield are discussed below. 

The two most extensive recent applications of PWO are in the electromagnetic calorime- 
ters of the CMS (ECAL) |19] and ALICE (PHOS) detectors at the CERN LHC. High 
quality of scintillators in mass production for these detectors has been achieved, however 
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further optimisation for the apphcation in PANDA appeared mandatory, because of the 
required much lower energy range. In PANDA the detection of 10-20 MeV photons is 
required which asks for the increase of hght yield in comparison with 9-11 phe/MeV (mea- 
sured with a PMT with bi-alkali photocathode at room temperature) typical for crystals 
used in these two experiments. 

2.4.1.1 Production 

The world largest producer of PWO crystals is the Bogoroditsk Techno-Chemical Plant 
(BTCP) near Tula (Russia). It is considered as the main producer of crystals for PANDA. 
The Czochralski method is used for crystal growth from a platinum crucible. To grow high 
quahty crystals an additional precrystallization is required. The polycrystalline PbWO^ 
obtained during precrystallization is sintered and used as a starting material for the final 
crystal growth. Up to seven crystals can be grown from the same initial load of the cru- 
cible by adding some new material between the crystallizations. Ingots of 250 mm length 
and of slightly elliptical cross-section, up to 45 mm in diameter became available. Crystal 
growth along the direction of the a axis appears to be the best for producing long samples. 
An alternative producer of PWO crystals, the Shanghai Institute of Ceramics of the Chi- 
nese Academy of Sciences (SICCAS) uses the modified Bridgman-Stockbarger technique 
in closed platinum foils for the crystal growth. This technique produces simultaneously 
several ingots. 

2.4.1.2 Scintillation mechanism 

The scintillation mechanism of PWO is determined by the electronic structures of the con- 
duction and valence bands in this material, which are described in details in [T]. Lumines- 
cence appears in PWO due to charge-transfer transitions in anionic molecular complexes. 
Both regular W0\~ and irregular WO^, tungstate groups are the luminescence centers. Be- 
cause the W0\~ complex has a Td point crystalline field symmetry, the final configuration 
of the energy terms is found to be ^Ti, and ^Ti, with the ^Ai ground state. The 
^Ai — j-^ Ti ,^ T2 transitions are responsible for the blue luminescence (420nm) of the lead 
tungstate crystal. When an anion vacancy appears in a W0\~ anionic complex the local 
symmetry of the new WO^ complex is reduced to Csj^. Luminescence in green (490 nm) is 
ascribed to this center. 
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2.4.1.3 Thermal quenching in PWO 

Most of the crystals of the tungstate family have an intense but slow emission in the mil- 
lisecond range. PbWO^ has a rather weak but fast emission because of a strong quenching 
of the scintillation process. This quenching is rather complex in origin but to a large extent 
associated with high-temperature charge transfer processes and thermal decomposition of 
the excited states. As a result, the light yield of PWO is strongly temperature dependent. 
In the range between +20° and —25° the dependence of light yield on the temperature is 
almost linear with a gradient of about —2 %/°C. 

2.4.1.4 PWO-II 

Detailed studies have been performed to improve the light yield of PWO scintillators. 
During the initial development of PWO scintillators it was established that an additional 
doping with molybdenum (Mo) and lanthanum (La) at concentrations < 100 ppm increases 
the light output, however it introduces the slow scintillation components with r = 1 — 4 fis. 
The further research was in direction of increasing the structural perfection of the crystals 
and reduction of the concentration of point structure defects by doping with yttrium (Y), 
lanthanum (La) or lutetium (Lu) ions, which suppress oxygen and cation vacancies in 
the crystal matrix. Structural perfection of the crystal has been reached at BTCP by an 
improved control of the stoichiometric composition of the melt. These crystals contain at 
least two times less Frenkel type defects leading to an increase of the light yield by 80 % 
compared to the average CMS quality. The new generation of crystals is called PWO- 
II (second generation) according to the accepted terminology. The scintillators produced 
following this technology will be used for the construction of PANDA EMC. An additional 
property of PWO-II crystals is their larger average temperature gradient of the light yield 
— 3%/°C, which is 50% higher than for the standard PWO crystals. A sample of the 
machined and optically polished PWO-II crystals produced for the PANDA electromagnetic 
calorimeter is presented in Fig. |2.5[ 

2.4.1.5 Optical properties 

The scintillation emission spectrum of PWO results from the superposition of two broad 
and complex emission bands at 420 nm and 500 nm respectively. The optimization of 
crystals leads to a Gaussian-shaped spectrum (140 nm FWHM) peaking at about 440 nm 
with a range from 360 nm to 570 nm at 10 % of the maximum. The scintillation kinetics 
is fast at room temperature and is dominated by an exponential decay with a constant of 
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Figure 2.5: Machined and optically polished PWO-II crystals (200mm length) for the PANDA 
electromagnetic calorimeter |51j . 

r=6.5ns. The possible presence of slow component was verified for PWO-II crystals [51] at 
+25°C and —25°C. The measurements demonstrated that 95 % of the scintillation light is 
collected within a time window of 300-400 ns. 

Crystalline material of lead tungstate is transparent in the visible spectral range and 
fully colorless. An optical transmission of the crystals has been steadily improved, in 
particular in the region of scintillation, between 360 and 570 nm. The optical transmission 
of PhWOi crystals can be limited by the presence of macroscopic defects like inclusions, 
precipitates or veils which scatter light in all directions, or by the existence of traps which 
induce absorption bands. A better control of the raw material preparation and of the 
growth and annealing conditions, as well as the introduction of dopants, have led to a 
considerable improvement in the optical transmission of the full-size PbWO^ crystals. 

As a result of the improved technology of the PWO crystals production, the typical 
values of light yield of 17-20 phe/MeV are observed for full size (200mm length) PWO-II 
crystals measured at room temperature with a photomultiplier with bi-alkali photocathode 
(quantum efficiency ~ 20 %). Reducing the thermal quenching of the luminescence process 
by cooling the PWO-II crystals from +25°C to — 25°C leads to an increase of light yield 
by a factor ^4. 

2.4.1.6 Radiation hardness 

The important characteristics of scinitallators in application for calorimetry is their be- 
havior in radiation environment. The dose load to the crystals from the hadronic probes 
during the operation of the PANDA experiment was estimated [51] with a Dual Parton 
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Model (DPM) event generator. With the assumption of event rate of 10^/sec expected 
for PANDA the absorbed energy rates of 27 mGy/h, 1.5 mGy/h, 0.16 mGy/h and 0.03 
mGy/h were calculated for the polar angles of 5°, 25°, 90° and 135°, respectively. 

It has been established previously [84], that the radiation load does not damage the 
scintillation mechanism of PWO. The main effect is the loss of optical transmission due 
to creation of color centers with absorption bands in the visible spectral region. The 
deterioration of the optical transmission in the experiment is caused by the interplay of 
damage and recovery mechanisms. The latter are fast at room temperature and keep the 
loss of light transmission at a moderate level. The operation of the PANDA EMC foreseen 
at — 25°C required additional studies of the optical transmission upon irradiation with 
7-rays at this low temperature. The two irradiation facilities at IHEP (Protvino, Russia) 
and JLU (Giessen, Germany) have been adapted for that purpose. The measurements 
[T13j have shown that due to a slow relaxation of the color centers in the cooled crystal one 
deals, asymptotically after a prolonged irradiation depositing 30-50 Gy, with a reduction of 
20 %- 30 % in the scintillation response. This reduction is compensated with a significant 
margin by the factor of ~4 increase in light yield due to cooling (Sec. 2.4.1.5) relative to 
room temperature. 



2.4.2 Photosensors (APDs) 

The operation of the PANDA EMC inside the 2 T magnetic field of the solenoid makes it 
impossible to use conventional photomultipliers for light readout. Two types of photosen- 
sors insensitive to magnetic field have been chosen for the PANDA EMC, silicon avalanche 
photodiods (APD) for the barrel part and vacuum phototriodes (VPT) for the forward 
endcap. The current section is mainly concentrated on the description of the operation 
principle and properties of APDs since the experimental part of the thesis presented in 
Chapter |3] is devoted to a study of this type of readout. 

Besides insensitivity to magnetic field APDs have additional properties, which motivate 
their application for the PANDA EMC: internal signal amplification and weak response 
to ionizing radiation. An APD optimised to detect the scintillation light from the lead 
tungstate crystals was developed for the CMS detector in collaboration with the Hama- 
matsu Photonics [19] factory in Japan. These APDs have the following features: com- 
pactness, with the overall thickness of less than 2 mm; fast rise time of the order of 2 
ns; a quantum efficiency of 70 - 80 %; insensitivity to magnetic field and low cost in mass 
production. In contrast to most APDs used for detection of red and infrared light, the 
inverted structure has been applied for the CMS APDs to overcome such drawbacks as a 
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Figure 2.6: Schematic view of an APD with reverse structure 



reduced gain for wavelengths shorter than 500 nm and high sensitivity to ionizing particles. 



Schematic view of an APD with an inverted structure is presented in Fig. 2.6 In these 
APDs light enters via the p'^'^ layer and is absorbed in the p"*" layer lying below, where 
electron-hole pairs are generated. The electrons then drift in the electric field towards the 
pn junction. They are amplified there by an impact ionization and drift in the n- or p- 
material towards the n^~^ electrode where the charge is collected. In front of the p^'^ layer 
is a passivation layer which both protects the wafer and provides an anti-refiection coating. 



This passivation layer is made of either silicon dioxide or silicon nitride as in Fig. 2^ The 
amplification or gain of the APD is largest for wavelengths which are completely absorbed 
in the p^ layer, in the present case for the wavelengths shorter than 550 nm. In APDs 
with the reverse structure the response to ionizing radiation is much reduced and to first 
order it is proportional to the thickness of the p~^ layer in front of the amplification region. 

The main disadvantage of the S8664-55 APDs in application for the PANDA EMC is 
their small active area 5x5 mm^ in comparison with the surface of the crystal exit face 
which is 21.4 x 21.4 mm?. The PANDA collaboration invested into a collaborative effort 
with Hamamatsu Photonics, which resulted in development of large area APDs (LAAPDs) 
with an active area of 10 x 10 mm^ and the same reverse internal structure as in the CMS 
APDs. At present, there is an APD of the third generation in development whose active 
area 10 x 20 mm will permit to cover completely with two APDs of this type the exit 
face of a PWO-II scintillator. Having two APDs on each scintillator gives an additional 
advantage of having the possibility to cover different dynamic ranges of the registered 
photons. Moreover, one has a back-up solution in case one of them fails due to radiation 
damage. 
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2.4.2.1 APD gain 

One of the key parameters which define the operation of an APD is its gain. The APD gain 
depends on the electric field applied across the p-n junction. The higher the reverse voltage, 
the higher will be the gain. However, if the reverse voltage is increased above the limiting 
value, a break-down of the p-n junction occurs. The current fiowing through the device 
series resistance and the load resistance cause the voltage drop and protect from permanent 
damage of the device |103j . Moreover, the gain of an APD is temperature dependent. At 
a fixed reverse voltage the gain decreases with the increasing temperature. This happens 
because of the stronger excitation of the lattice vibrations - phonons, provoking more 
frequent collisions with the charge carriers and preventing them from reaching sufficient 
energy to trigger ionization. The temperature dependence of gain has a value of - 2.2 %/°C 
at an internal gain of M = 50 for a CMS- APD |109] . This dictates the requirement of 
stability of the operating temperature for the PANDA EMC at the level AT = iO.TC. 

2.4.2.2 APD noise and excess noise factor 

Since an APD is designed to operate under a reverse bias, its sensitivity at low light levels 
will be limited by the shot noise and the APD's leakage current. Shot noise derives from 
the random statistical Poissonian fiuctuations of the dark current. Id- The dark current 
Id of an APD can be divided into the bulk current, J;,, and the surface current, Ig- While 
the surface current is independent of the applied gain, M, the value of the bulk current 
increases with increasing gain, hence the total dark current can be defined as: 

h = h-M + h. (2.1) 

In addition, avalanche photodiodes generate excess noise due to the statistical nature 
of the avalanche process. These statistical fiuctuations of the APD gain are characterised 
by an excess noise factor, F, and have its origin in inhomogeneities in the avalanche 
region and in hole multiplication. a/F is the factor by which the statistical noise of APD 
current (equal to the sum of the multiphed photocurrent plus the multiplied APD bulk 
dark current) exceeds that which would be expected from a noiseless multiplier on the basis 
of Poissonian statistics (shot noise) alone. The total spectral noise current for an APD is 
given by: 

in = [2q{h + h-M^ ■F)B]"\ (2.2) 

where q is the electron charge and B is the bandwidth. 

The excess noise factor is a function of the carrier ionization ratio, k, where k is usually 
defined as the ratio of hole to electron ionization probabilities and is given by the following 
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equation: 

F = k-M + i2-k)il-j^), (2.3) 

and therefore the lower are M and k, the lower is excess noise factor. 

The excess noise factor measured for the APD S8664-1010 is F = 1.38 at the gain 
M = 50 and F = 1.57 at the gain M = 100 [SI]. 

2.4.2.3 Radiation hardness 

The effect of radiation on APD's performance is a critical issue for their operation in 
electromagnetic calorimeters. Two main damage mechanisms can take place. The first 
one is a bulk damage, due to displacement of Si atoms from their lattice sites, causing an 
increase in dark current. The second one is a surface damage, or the creation of defects 
in the front layer, which may lead to an increase of the surface current and a decrease 
of the quantum efficiency. The detailed studies of APDs radiation hardness have been 
performed with different sources of radiation at T = — 25°C, which corresponds to the 
foreseen operation conditions of the PANDA detector. The change of surface and bulk 
dark current induced by irradiation was studied with the proton beam at KVI, with an 
intense ^°Co 7-source at Giessen and with an "^^^Am — a — Be neutron source at Frankfurt. 
It was shown that the proton induced damages are mainly located inside the conversion 
layer and at the end of the avalanche region. The irradiation of APDs with an integrated 
fluence of protons 1.1 ■ 10^^ p which is comparable with the dose expected in 10 years of 
PANDA operation leads to a decrease of the APD gain by 50 - 60% [51]. Measurements 
with photons and neutrons give no indication for any kind of induced bulk damage and it 
is concluded that an increase of dark current is caused by the surface damage only. 

2.4.2.4 Vacuum phototriodes (VPT) 

The photon detection in the forward endcap EMC has to deal with rates up to 500 kHz 
and magnetic fields up to 1.2 T. Therefore vacuum phototriodes (VPT) were chosen with 
a diameter of the photocathode of 22 mm as photon sensors for this part of the EMC. 
The main reasons for this choice are rate capability, radiation hardness, absence of nuclear 
counter effect and absence of temperature dependence of the efficiency. In contrast to 
the barrel region, the magnetic field is oriented in the axial direction of the VPTs and 
thus makes it feasible to use VPTs for the endcap. Vacuum phototriode is essentially a 
photomultiplier tube with only one dynode and a weak dependence of the gain on the 
magnetic field. The VPTs with standard quantum efficiency of 15 % and gain of about 
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Figure 2.7: The barrel and forward endcap of the PANDA EMC [5T]. The antiproton beam is 
arriving from the left-hand side of the drawing along the symmetry axis of the barrel EMC. 



10 are available from Hamamatsu and the Research Institute Electron (RIE). However, 
Photonis is developing phototriodes with high quantum efficiency photocathodes and high 
gain >40 at 1000 V, which are more suitable for PANDA. 



2.4.3 Coverage of the EMC. Mechanical design 

The maximal coverage of the solid angle will be satisfied by dividing EMC into three parts: 
barrel and two endcaps. The barrel EMC, consisting of 11360 crystals, will cover polar 
angles from 22° till 140° and will have an inner radius of 57 cm and an outer radius of 94 
cm. The forward endcap, located at 2.05 m downstream from the target, will consist of 
3600 crystals and will cover polar angles from 23.6° down to 5° in the vertical and 10° in the 
horizontal direction. The backward endcap in the current design consists of 592 crystals, 
covers polar angles from 151.4° till 169.7°, and is located at 55 cm upstream from the 



target. The barrel and the forward endcap of the PANDA EMC are presented in Fig. 2.7 



following the Technical Design Report for the PANDA Electromagnetic Calorimeter [51] . 

The individual shape of a crystal in the barrel part is a tapered parallelepiped, with 
dimensions depending on the position along the beam axis. The front face of an individual 
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Figure 2.8: Geometrical arrangement of the crystals of the barrel in a cut along the beam axis 
[51]. Numbers indicate different crystal types. 



crystal is close to the square with the 20 mm side (21.18-22.02 mm), which corresponds 
to Moliere radius of PbWO^. The barrel part is divided into 16 identical sectors, each 
subtending 22.5° of the azimuthal angle, called slices. The crystals are grouped into packs 
of 4 X 10 along the azimuth. Along its length the barrel contains 71 crystals of 11 types 



depending on the position, as indicated in Fig. |2.8[ In one shce of 710 crystals, 3 or 5 alveole 
packs are grouped together into 6 modules: 4 modules of 120 crystals each, 1 module with 
70 and 1 with 160 crystals. The off-pointing geometry of EMC is foreseen, i.e. the crystals 
are not pointing towards the target position. This ensures, that tracks originating at the 
target never pass through gaps between the crystals, but always cross a significant part of 
a crystal. 

The dimensions of the crystals in the forward endcap are determined by the requirement 
to accommodate the VPT readout on the rear surface of the crystal. Currently available 
VPTs with the diameter of 22 mm can be accommodated on a crystal with a rear cross 
section of 26 x 26 mm^. The off-pointing geometry of the forward endcap EMC, i.e. an 
orientation of the crystals towards a point on the beam axis 950 mm upstream the target, 
determines that each crystal has a front-face of 24.4 x 24.4 mm^. A cluster of 4 crystals 
will form a mini-unit of trapezoidal cross section, mounted in a single carbon-fiber alveole. 
Four mini-units of crystals will be combined to form a 16-crystal subunit of ca. 19 kg that 
can be attached individually to the 30 mm thick aluminum mounting plate. 

To optimize light collection crystals will be wrapped with a reflective material consid- 
ered as Radiant Mirror Film ESR from 3M in the current design. 
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The cooling of crystals down to -25°C is foreseen with stabihty at the level of ±0.1°C, 
which will be satisfied by a cooling system in combination with a thermal shielding. The 
requirement of coohng power is estimated at Ri3300 W. 

2.4.4 Electronics 

According to the general scheme of the Data AcQuisition (DAQ) for the PANDA exper- 
iment, the readout of EMC assumes continuous digitization of the amplified and shaped 
signals from APD and VPT photosensors. The readout chain will include low-noise front- 
end electronics, digitizer modules and data multiplexers. 

The requirements to frontcnd electronics are determined by several factors. The detec- 
tion of photons with energies from the maximum 12 GeV per crystal down to 3 MeV, which 
is the single crystal threshold with the electronics noise level of 1 MeV requires the dynamic 
range of 12000. Expected event rates up to 500 kHz for the forward endcap and up to 100 
kHz for the barrel require short shaping times. Operation of the frontend electronics inside 
the cooled volume implies strict conditions on power consumption. Two complementary 
low-noise low-power charge sensitive preamphfier circuits have been developed. The first 
one, based on discrete components with very good noise performance using signal shaping 
with a peaking time of 650 ns, is foreseen for readout of the forward endcap EMC, where 
event rates up to 500 kHz are expected. A CMOS ASIC, based on 350 nm technology, 
will be used for the barrel EMC. It has a similar noise performance with a shorter peaking 
time of 250 ns and is designed for event rates up to 350 kHz. 

The digitizer modules will be located at a distance of 20-30 cm and 90-100 cm for 
the barrel EMC and the forward endcap EMC, respectively, away from the analogue cir- 
cuits and outside of the cold volume. The digitizers consist of high-frequency, low-power 
pipelined ADC chips, which continuously sample the amplified and shaped signals. The 
sampling is followed by the digital logic, which processes time-discrete digital values, de- 
tects hits and forwards hit-related information to the multiplexer module via optical fiber 
links. At this step the detection of clusters of energy deposition can be efficiently imple- 
mented. The multiplexer modules will be located in the DAQ hut and they will perform 
advanced signal processing to extract amplitude and signal- time information. 

2.4.5 Performance 

The main task of the electromagnetic calorimeter is reconstruction of photons and leptons 
with high efficiency. EMC measures the deposited energy of a particle, E, and provides 
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information on its direction of propagation by determining the point of impact [6, </>). 
Therefore, the energy and spatial resolution are the main parameters which define per- 
formance of the EMC. Besides the energy resolution, the energy coverage is an important 
characteristics of the EMC and especially the low energy threshold. In addition EMC 
can provide the time information and therefore the time resolution of the order of 1 ns 
achievable with the PWO scintillators can be used for event separation. EMC can also 
provide separation from the vr^ background by an analysis of the lateral shower shape. 
Such separation is possible because hadronic showers have different lateral shower shape 
defined by the difference in energy loss per interaction and the elementary statistics of the 
underlying processes. 

2.4.5.1 Energy threshold 

In study of radiative decay of the charmonium states produced as resonances in p — p anni- 
hilation the final states with vr^'s and r/'s produced non-resonantly are a significant source 
of background photons for the direct photons from the transitions in charmonium. If one 
of the photons from tt^ (or rj) decay has energy below the threshold, the 7r° (or rj) is not 
reconstructed. The second photon will be erroneously classified as coming from the reso- 
nance decay thus adding to the background. Therefore the probability of reconstruction 
depends on the the minimum photon energy Ethres accessible with the EMC; Ethres=^0 
MeV is a design goal for PANDA EMC. A Monte Carlo based study of this problem is 
presented in Chapter |4j 

2.4.5.2 Energy resolution 

The energy dependence of the Gaussian energy resolution as/E can be parameterized in 
the following way: 

^ = a®^ ^ , (2.4) 

where a is a constant term and 6 is a stochastic term. The stochastic term is mainly defined 
by photostatistics and can be expressed as 



N 



(2.5) 



where Npe is the number of electron-hole pairs released in the photodetector per GeV 
energy deposited in the scintillator and F is an excess noise factor, which parametrises 
fluctuations of the gain process in an APD. 
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The constant term has contributions from leakage of energy out of the scintilator, 
non-uniformity of light collection and intercalibration errors. The design requirements for 
PANDA EMC are a < 1 % and 6 < 2 % [H]. 

2.4.5.3 Position resolution 

Position resolution is mainly defined by granularity of the electromagnetic calorimeter. The 
reconstruction of an impact point of a particle is performed by defining the cluster, i.e. a 
group of neighboring crystals with the deposited energy above the single cluster threshold 
of Eciuster = 3 McV, and the following weighting of the position of the crystal with the 
deposited energy. Different algorithms of position reconstruction can be applied, which 
take into account also the lateral shower shape. The good space resolution is crucial for 
the splitting of clusters created by overlapping photons, for example from 7r° with a small 
opening angle. For such a splitting the local maxima produced by two photons should be 
separated by at least two crystal widths. 



Chapter 3 

Scintillator and APD test 
measurements 



3.1 Introductory remarks 

Interest of our equipe in development of the electromagnetic calorimeter was motivated by 
the proposed experiment related to 7-ray detection in the h^. — )■ r/^ transition in charmonium 
[52] and our former experience of work with cooled APDs 



Sec. |2.4| outlined the present status of the EMC correlated in time with writing of 
this Thesis, whereas when the work of its author just started (2004) many of the details 
considered now as self-evident, had to be still developed. We had realized that low-energy 
response of scintillators will be of interest for the future PANDA-EMC. The tools available 
at SINS for such a project were: a C-30 compact cyclotron and a Van-de-Graaff accelerator. 
The C-30 cyclotron accelerates protons up to about 26 MeV. Having a source of quasi- 
monoenergetic particles close to the desired energy, it seemed natural to start the project 
using the latter tool. It was obvious, however, that proximity of the target room to 
the cyclotron vault might adversely affect measurements, due to the cyclotron RF noise 
pick-up. It was decided, therefore, to readout a scintillator from both exit surfaces in 
coincidence, in expectation that in this way the potential source of pick-up noise will be 
suppressed. Moreover, this mode of operation makes it possible to study the time response 
of a scintillator and readout chain. 

The calorimeters existing at that stage suggested PWO [IS], [IS] and BGO [117] as 
the scintillators for possible application in the PANDA EMC, while LYSO was considered 
as a candidate for the EMC in the future detector at the ILC [45j. An APD S8664-55 
of Hamamatsu used for readout in the CMS and ALICE detectors at CERN had by far 
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3.2 Measurements with ^ 25 MeV protons at the C-30 cyclotron 



too small sensitive area of 5 x 5 mm to satisfy the low detection threshold required for 
PANDA. We have, therefore, started our measurements with APDs of Advanced Photonix 
Inc. (USA) which had sensitive area of ~ 2cm^; R&D work on an APD with the area of 
1 cm^ started almost simultaneously in a collaborative effort of PANDA with Hamamatsu 
Photonics Inc. We have tested several combinations scintillator-APD in the measurements 



with 25 MeV protons reported below in Sec. 3.2 



3.2 Measurements with ~ 25 MeV protons at the C- 
30 cyclotron 

3.2.1 Experimental setup 
3.2.1.1 The C-30 cyclotron 

The proton beam used in the measurements was produced with the C-30 cyclotron at SINS, 
Swierk. The C-30 cyclotron is a compact, fixed magnetic field, fixed frequency isochronous 
machine. Proton beams are obtained by accelerating H~ ions from the ion source followed 
by electron stripping in an aluminium foil, which leads to beam extraction. The 52 MHz 
radiofrequency electric field is used to accelerate protons. The extraction geometry and the 
magnetic field of 1.7 Tesla of the cyclotron magnet give the possibility to obtain a range 
of proton velocities extending from 0.2c up to 0.24c. This corresponds to proton beam 
energies from 19 MeV to 28 MeV, which is determined by a non intrusive time-of-flight 
system. An average external beam intensity up to 10 fiA can be achieved. 



3.2.1.2 Details of the experimental setup 

The experimental setup consists of a 070 cm diameter vacuum chamber with a 2 mg ■ cm~'^ 
Au target on a thin C backing in its center and a small chamber of 015 cm inner diameter 
fixed to the 30° exit port of the former. The small chamber houses the tested scintillator 
in the form of parallelepiped 20 x 20 x 40 mm^. The scintillator is readout from both 
20 X 20 mm^ end surfaces with a PMT XP2020 and an APD. The cathode pin of an APD 
is soldered with a thin covar wire to minimize heat leakage to the SHV male connector 
installed in the body of the small chamber. The connector provides signal path to and HV 
from a preamp. To enhance light transport from its origin in the scintillator volume to the 
APD sensitive area the scintillators were coated with a few layers of a teflon tape. The 



horizontal section of the setup is presented in Fig. 3.1 In this figure the small chamber is 
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Figure 3.1: Horizontal section of the experimental setup for measurements of energy resolution 
of scintillators using protons accelerated with the C-30 cyclotron. 



shown rotated by 90°, i.e. the PMT axis which is actually directed towards the floor (see 



Fig. 3.2) is drawn rotated into the plane of Fig. 3.1 to present the entrance slit and the 



scintillator location. 

The APD is cooled using two thermoelectric modules (Peltier elements DT12-4[^ (see 



Fig. 3.3 for details). The Peltier elements are fixed between the two side copper plates and 
the central plate, which is silver-soldered in between the arms of an U-shaped, water-cooled 
copper tube, presenting a heat-sink for the "hot" sides of the Peltier elements. The side 
plates are in thermal contact with the Peltiers' "cold" sides, which provide cooling of the 
LAAPD-holder via flexible copper joints (cold-fingers). In order to improve the thermal 
contact in vacuum, the contacting surfaces are polished and covered with a thin layer 
of thermoconducting (but electroisolating) white paste P12[^ The temperature control is 
provided by Pt-100 resistive thermometers installed in two points: the first on at the center 
of the right contact plate and the second one on the LAAPD holder. The small chamber is 
pumped down to 1. — 2. ■ 10"^ mbar together with the large one; vacuum effectively isolates 
all cooled elements from ambient and prevents from moisture deposition. 



iMarlow Industries, Inc., 10451 Vista Park Road, Dallas, TX 75238-1645, USA 
^Wacker Chemie GmbH, Germany 
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Figure 3.2: The experimental setup seen from the side containing the cooling elements (compare 



with Fig. 3.3). Seen are the two plastic tubes providing the cooling water inlet and outlet from 



the U-shaped copper tube (9); whereas in between them are seen the cables of power supply to 
the Peltier elements in a flexible metal screen. 
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Figure 3.3: Vertical section of the small chamber of the experimental setup, which houses 
the scintiUator. The constituent elements are: (1) APD, (2) scintihator, (3) PMT XP2020, (4) 
Teflon plate, (5) LAAPD holder, (6) Flexible copper joint (cold- finger), (7) Side copper plates, 
(8) Central plate, (9) Water-cooled U-shaped copper tube, (fO) Peltier elements DT12-4, (11) 
Vacuum feedthrough connecting to Pt-100 thermometers, (12) Vacuum feedthrough for power 
supply to DT12-4s, (13) Soft-iron housing of the PMT, (14) 7 or ^-ray source, (15) Pumping 
port to the chamber with a defining slit. 
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Material 


Density 


Radiation length 


LY 


Decay time 




{g ■ cm^^) 


(cm) 


% (LY Nal) 


(ns) 


PWO-I 


8.3 


0.89 


0.3 


40/10 


PWO-II 


8.3 


0.89 


0.6 


40/10 


BC-408 


1.03 


43 


24 


2.1 


BGO 


7.1 


1.12 


9 


300 


LYSO 


7.1 


1.2 


75 


40 



Table 3.1: Main parameters of the scintillation materials used in the measurements with 25 
MeV protons. 



3.2.1.3 Tested APDs and scintillators 

The main purpose of the study was to measure energy and time resolution for PWO scintil- 
lators with an APD readout. However, at the moment when these measurements have been 
performed BGO was considered as an alternative option for the PANDA electromagnetic 
calorimeter. It has an advantage of higher light yield, but has much slower fluorescence 
decay time in comparison with PWO. The plastic BC-408 is a very fast scintillator with 
large light output well matched to the maximum of spectral sensitivity of many APDs. It 
has been used for tuning fast timing in our experiment. In addition several measurements 
with a LYSO scintillator have been performed. LYSO is characterised by a larger radiation 
length than PWO, however it has an advantage of much higher light yield having at the 
same time a comparable speed. Its high price does not permit its application in calorimetry 
at present, however LYSO is considered as a perspective material for future applications, 
e.g. in the EMC at ILC [15]. The main parameters of the mentioned scintillators, such as 



density, radiation length, light output and decay time are summarized in Table |3.1[ For 
PWO both fast and slow components of the fluorescence decay time are listed. The two 
generations of PWO are introduced in the table; PWO-II is characterized by about twice 



higher light yield, it is discussed in more detail in Sec. 2.4.1 Only PWO-II has been used 



in our measurements. Light yield of a scintillator is the parameter which may vary from 



crystal to crystal, therefore the values listed in Table 3.1 are only indicative. 

Three different types of APDs were used during our measurements; two produced by 
Advanced Photonix Inc. (API): 630-70-74-510 and 630-70-73-510 and one produced by 
Hamamatsu Photonics: S8664-1010 . The Advanced Photonix APDs differ from each 
other by the quantum efficiency dependence on the wave-length. The 630-70-74-510 has 
maximum sensitivity in the blue part of light spectrum while 630-70-73-510 in ultraviolet. 



3 Scintillator and APD test measurements 



69 



100 
90 

ao 

70 
' 60 
I 5* 
I 40 
'' 30 
SO 
10 



EftBElive □ua.Hum Elfidency Gain 300 va. Wavelengt"! 



350 



: 












































































UV 

































































4&0 



5M 



6S0 7S0 



950 



10S0 



Figure 3.4: Dependence of the APD quantum efficiency on the wavelength for the two types of 
Advanced Photonix APDs 630-70-74-510 ("blue") and 630-70-73-510 ("UV"). 



see Fig. 3^ The APDs by Advanced Photonix are fabricated with three diameters of 
the sensitive area: 5, 10 and 16 mm. For our measurements the largest diameter type 
630-70-74 has been selected which has an area of ~2 cm^. This area is one half of the cross 
section area of the scintillator, thus about 50% of the light emerging from its exit face can 
be intercepted and processed. The sensitive area of Hamamatsu S8664-1010 is a square 
with 10 mm side. The studied products of these two companies also differ strongly in the 
required working bias voltage. A typical working bias for the Advanced Photonix APDs is 
1760 - 1870 V, with a gain of ~50 at the lower and ~300 at the upper limit of the interval, 
whereas Hamamatsu APDs work at a significantly lower bias of ~ 400 V. 



3.2.1.4 Readout electronics 

The readout electronics was a fast-slow coincidence system assembled using mostly the 
standard ORTEC NIM modules. The complete block-scheme is presented in Fig. |3.5 
Two signals from the PMT have been employed, the signal from the lO*'' dynode was 
amplified with an ORTEC 672 amplifier and after activating the Single Channel Analyser 
(SCA) was directed to the coincidence module together with a signal from the APD to 
produce a gate signal for the two-dimensional analyser. The second signal from PMT was 
collected from its anode and was used, after shaping with the aid of the constant-fraction 
discriminator ORTEC 935, as a start signal for time distribution measurements with the 
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Time-to-Amphtude Convertor (TAG) ORTEC 567. The stop signal for the TAG was 
derived from the fast timing branch (denoted "T") of the GSI GATS A preamplifier. The 
spectroscopic output from the preamplifier (denoted "E" in Fig. 3.5) was amplified with 
the standard ORTEG 452 amplifier, whose one output activated the SGA. If a coincidence 
occured with the "slow" PMT signal, the gate of the two-dimensional analyser is opened 
to accept the APD amplitude and time difference signals. The timing signal from the 
preamplifier after an additional amplification with the Timing-Filter-Amplifier was used 
as a stop signal for the TAG. The signal from TAG was fed into the second input of the 
analyser. To optimize time performance the timing signals from both APD and PMT 
were fed into the Gonstant Fraction Discriminators (GFD) with optimised settings of zero 
crossing and threshold level. As it is clear from the foregoing presentation, the registered 
events had two coordinates: the amplitude of an APD signal and the time difference 
between PMT and APD signal arrivals. This permitted to judge, neglecting the time 
dispersion of XP2020, on dependence of the APD signal time dispersion on its amplitude. 
Small signals, which might have been caused by background effects, were suppressed by 
setting appropriate thresholds on the signal amplitudes. A typical PMT spectrum from 
the 10th dynode of XP2020, with the level of threshold marked with a vertical broken line. 



is shown in Fig. |3.6[ The delay line amplifiers were used wherever it was necessary to 
synchronize the signals. 



3.2.2 Results 

The measurements were started with the plastic scintillator BG-408. Due to its high light 
yield it allowed to estimate the contribution of an APD to energy and time resolution of 
the combination scintillator plus APD from an APD only, i.e. in this case the number of 
photons should not be the limiting factor. Before the measurements with a proton beam the 
settings of the GFDs were adjusted using radioactive sources for optimum time resolution. 
The measurement with a plastic scintillator BG-408 has been performed with only one type 
of an APD, Advanced Photonix 630-70-74-510, which has maximum efficiency in the blue 
part of light spectrum. 

The time calibration of the time-to-amplitude converter has been performed with an 
ORTEG time calibrator, which produces start and stop signals separated by an integer 
multiple of a fixed time period, which has been selected 10 ns for our purposes. 

With the 26 MeV proton beam the coincidence 2-dimensional spectra have been col- 



lected. A sample 2-dimensional spectrum is shown in Fig. 3.7 (a) and in Fig. 3.7 (b) 



its projection on the APD amplitude axis. Because of its asymmetric shape the energy 
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Figure 3.5: Block-scheme of electronics used in the measurements with 25 McV protons. The 
hatched area indicates the region of ~ 4 mm width illuminated with protons (not to scale with 
the 40 mm length of the scintillator). 
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Figure 3.6: PMT threshold for coincidence measurements. 
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dispersion a is determined as the full width at half maximum (FWHM) divided by 2.35. 



The two-dimensional spectrum is presented in Fig. 3.7 (a). The time resolution has been 



estimated by projecting that part of the two-dimensional spectrum which corresponds to 
the signal from the proton beam, i.e. having APD amplitudes between the channels 280 
and 440. 



The projection of Fig. 3.7| (a) on the time axis, limited to APD amplitudes 280-440, 



containing clean proton hits, is presented in Fig. 3.7 (c). The time resolution, defined as 
FWHM is equal 0.37 ns. The two numbers extracted from the results of Fig. 3.7 (a) are 



representative for a fast scintillator with large light output. One may conclude that under 
these conditions it is feasible to reach subnanosecond timing using an APD. These results 
will be further compared with those obtained under less favorable conditions. 

The measurements with PWO-II scintillator which was the main purpose of this study 
have been performed with all the three types of APDs. An additional aspect of the mea- 



surements was a study of the effect of APD cooling on the APD gain. In Fig. 3^ the 
two spectra from PWO plus Advanced Photonix APD are presented. The upper one was 
measured with an APD at room temperature and the lower one with APD at — 20°C. We 
observe the APD gain increase due to cooling by a factor of 2.5. The dark current of the 
APD was controlled during the measurements and its decrease from 50-60 nA at room 
temperature down to 3-4 nA at — 20°C is observed. This indicates that for amplitude and 
time resolution an APD contribution to electronic noise must have also decreased. 

The APD amplitude spectra from the two types of Advanced Photonix Inc. (API) 
APDs (blue and UV) are presented in Fig. 3.9 a) and b). These two measurements have 
been performed with slightly different energies of the proton beam, 26 MeV with API 
"blue" and 22 MeV with API "UV". The obtained amplitude dispersions are determined 
as FWHM/2.35 and a/ E values are equal 17 % and 22 %, respectively. For the spectrum in 
Fig. |3.9 b) the fit assuming a gaussian peak plus a linear background has been performed. 
The difference in resolution seems reasonable since the worse energy resolution for APD 
"UV" is expected due to the lower quantum efficiency of the latter over PWO emission 
spectrum (see Fig. 3.4). Fig. 3.9 c) presents the time spectrum. The width at half maximum 
is equal 0.95 ns, a loss in resolution in comparison with Fig. 3.7| c) not unexpected in view 
of the decreased e-h statistics. 

The second series of measurement has been performed with an APD produced by 
Hamamatsu Photonics Inc. This APD possessed 10 x 10 mm sensitive area and was 
specially designed for the PANDA needs. However its area is twice smaller than that of 
API APDs which affects the measured energy resolution due to inferior light collection. In 



Fig. 3.10 a the amplitude spectrum with the extracted energy resolution is presented and 
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Figure 3.7: The two-dimensional spectrum from a BC-408 scintillator with an APD readout 
obtained with the 26 MeV proton beam (a) with its projections on the APD amplitude axis with 
a/E determined from the peak FWHM (b) and on the time axis limited to APD amplitudes 
280-440 (c). 
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Figure 3.8: The ampHtude spectra from PWO+APD at room temperature (upper spectrum) 
and cooled down to — 20°C (lower spectrum). 
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Figure 3.9: The response of PWO scintillator to the proton beam, a) energy spectrum taken with 
an API APD "blue", b) energy spectrum with an API APD "UV", c) time spectrum measured 
with an API APD "UV". 
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API("Blue") 


API("UV") 


Hamamatsu 


PWO-II 


17 ± 2 % 


22 ± 3% 


35 ± 3 % 


BC-408 


5.6 ±0.7% 






LYSO 






11.9 ±0.9% 


BGO 






9.8 ±0.8% 



Table 3.2: Summary of the measured energy resolutions for scintillators with (cooled) APD 
readout for the different combinations scintillator- APD. 



the time spectrum with its resolution is presented in Fig. 3.10 (b). 

The energy resolution worsened in comparison with Fig. 3.9tocr/i? = 35%. A comment 
on this observation will be made in the next subsection. On the other hand the time 
resolution of 1.0 ns does not differ very much from the one previously obtained with the 
API APD. 

BGO and LYSO were the two additional tested scintillators. In both cases the Hama- 



matsu S8664-1010 APD has been used. Fig. 3.11 a) presents the two-dimensional time- 
energy spectrum obtained in the measurements with BGO. Fig. 3.11 b) and c) represent 
energy and time spectra correspondingly which are the projections of this 2-dimensional 
spectra, the events are grouped in the three areas with the extreme top and bottom area 
equidistant from the central one. The time distance corresponds to the period of the 
cyclotron high frequency. Apparently, the TAG start and stop signals are derived from 
different cyclotron bunches in these extreme groups of events. The estimated amplitude 
resolution is based only on the events in the central group. This is equal to 5.6 % and the 
time resolution is 1.32 ns. 



Fig. 3.12 presents the spectra from the APD (a) and from the PMT (b) for the LYSO 
scintillator. The energy resolution for the APD is 11.9 % and for the PMT is 5.9%. 



The summary of the measured values for energy resolution is presented in Table |3.2 
and for time resolution in Table 13.31 



3.2.3 Interpretation of the results. Simulation of the light col- 
lection efficiency 

The measured value of the time resolution for PWO scintillator with an APD readout 
0.95-1.0 ns appears very promising for the purpose of the PANDA experiment where the 
time resolution below 1 ns is required. The quantity measured in our experiment is the 
dispersion of signal arrivals in time between APD and PMT. The PMT XP2020 has negli- 
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Figure 3.10: Response of the PWO scintillator readout with the Hamamatsu S8664-1010 APD 
to 26 MeV protons, reflected in a) energy spectrum and b) time spectrum. 





API("Blue") 


API("UV") 


Hamamatsu 


PWO-II 
BC-408 
LYSO 
BGO 


0.95 ± 0.01ns 
0.37 ± 0.01ns 




1.0 ± 0.01ns 
1.32 ± 0.01ns 



Table 3.3: Summary of the measured time resolutions for scintillators with (cooled) APD readout 
for the diff'erent combinations scintillators- APD. 
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Figure 3.11: The response of BGO scintillator to the proton beam, a) 2-dimensional time-energy 
spectrum, b) energy spectrum with Hamamatsu APD, c) time spectrum. 
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Figure 3.12: The response of LYSO scintillator to the proton beam, a) energy spectrum with 
Hamamatsu APD, b) energy spectrum with PMT. 
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gible signal dispersion in comparison with an APD. The time dispersion, o"time, of an APD 
is proportional to the ratio of two factors: 

C^time ~ CTnoisc/ {dV/ dt) , (3.1) 

where o"noise is the rms voltage noise and dV/ dt is the signal slope at threshold of the timing 
discriminator. In view of the rather low gain of ~300 of the applied APD in contrast to 
~10^ typical for PMTs the contribution of APD to the measured value of time resolution 
is dominating over the PMT. 

The measurements of the time resolution of full-size 200 x 20 x 20 mm^ PWO-II scin- 
tillator readout by Hamamatsu APD and with 10 bit 80 MHz sampling ADC used as a 
digitizer performed recently at the University of Groningen |51j gives the number of the 
order of 1 ns, which is comparable with our results. 

On the other hand the energy resolution for the PWO scintillator with Hamamatsu 
APD readout appears worse than we expected. The quenching effect was verified as a 
possible source of this discrepancy. Measurements of the energy resolution with low energy 
proton beam, which we have performed, can be compared to energy resolution for 7-rays 
on the basis of Birks' theory [31] . This theory empirically describes the effect of quenching 
of light yield for highly ionising particles. The light yield dL over the element dx of path- 
length is related to the energy loss dE by the following equation: 

dL/dx = S^^^, (3.2) 

dx 

where kB is a Birks' constant and S" is a scintillation efficiency. The total light yield 



produced by a particle is obtained by integration of Eq. 3.2| from its initial energy down 



to zero. The Birks' formula Eq. |3.2| is based on the assumption that the incident particle 
excites the number of molecular structures in scintillator, proportional to the stopping 
power, damaging a fraction BdE/dx of them. The latter are assumed to act as light- 
quenching centers with an efficiency k. The quenched fraction goes to zero in the limit 
of low ionization density, therefore the effect should be largest for incident ions with low 
energy. Although the high energy ions eventually slow down, the average dE/dx along the 
particle track is less for high energy incident particles. The factor kB is material dependent 
and should be extracted from measurements for each scintillator and particle type. For 
the PWO scintillator Ref. [75] presents the results of measurements with 85 MeV protons. 
A comparison with response to cosmic muons, which suffer a negligible quenching, yields 
the quenching factor Qf = 0.92. For the 25 MeV proton the quenching factor has been 
rescaled with the Birks' formula via its numerical integration over the energy range of 
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interest using the tabulated values of dE/dx. This estimate gives Qp = 0.9 for 25 MeV 
protons in PWO. This number tells that the quenching effect alone cannot be considered 
as the main source of poor energy resolution in our experiments. 



It has been argued in Sec. 'lA that the two main contributions to energy resolution are 
electronic noise and photon statistics. One of the factors which influence photon statistics 
is light collection efficiency. We have performed light transport simulations for our setup 
with the Litrani program to determine light collection efficiency. 

The setup implemented in simulations, consisted of a crystal, a photomultiplier and an 
APD. The 20 X 20 X 40 mm^ crystal was oriented with its long axis along the OZ axis. The 
PMT was placed under the crystal, as in Fig. 3^, and an optical contact with its entrance 
window was via an optical glue with a thickness of 0.03 mm. Its refractive index was not 
known, therefore for the simulations n = 1.704 was assumed as possessed by Meltmount, 
the optical couphng material implemented in Litrani. It should be stressed however, that 
changing the refractive index in the range 1.6-1.7 does not change the result of simulations 
significantly. The bottom side of the crystal was completely covered by the PMT entrance 
window which is 2 in diameter. The APD with square sensitive domain of 1 x 1 cm^ was 
located above the crystal not having an optical contact with it. Because of the limited 
cooling power of the used Peltier elements, such a contact would immediately deteriorate 
the minimum temperature reached by our APD. The size of the gap was of the paper card 
thickness and was assumed 0.3 mm in simulations. The teflon coating of the crystal was 
treated as a perfect diffusor in simulations. 

An area illuminated with protons extended over the band of 4 mm width along the 
length of the scintillator at its midplane and over its entire width, being deflned by the 
collimator at the entrance to the small chamber. The depth of proton penetration was 1.6 
mm, determined by the stopping power of PbWO^. The fluorescence light was generated 
with a wavelength of 420 nm, which is one of the two main components of PbWO^ scintil- 
lation. The absorption length was assumed the same as for the CMS PWO crystals, which 
is implemented within the program and is equal 18.5 cm at the 420 nm wavelength. The 
quantum efficiencies of PMT and APD were taken as 20 % and 70 %, respectively. 

The main output of the simulations with Litrani is the calculated ratio of the number 
of produced photoelectrons in PMT and e-h pairs in APD. This ratio is determined by 
the light collection efficiency for given geometry and quantum efficiencies of the given pho- 
todetector. It is equal 6:1 for PMT vs. APD, despite almost equal product of quantum 
efficiency and covered area. This is the effect of a lack of optical contact between the scin- 
tillator and the APD. To estimate the statistical term of the energy resolution the number 
of registered photoelectrons in APD should be calculated. The standard characteristics of 
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the hght yield of a scintillator is the number of photoelectrons per MeV measured with 
the PMT with bi-alkali photocathode with quantum efficiency of 20 %. For PWO-II crys- 
tals a typical value is 20 photoelectrons per MeV. This number of photoelectrons can be 
assumed to be equal to what is measured with PMT in our setup due to similar condi- 
tions of measurements. The number of photoelectrons registered in APD can be estimated 
as a number of photoelectrons registered in PMT multiplied by the ratio obtained from 
the Litrani simulation. One thus obtains 3.3 e-h pairs per MeV energy deposited in the 
scintillator and emitted as light. 

The statistical term in energy resolution is: 



E y Ne-h y/E{MeV) ' 



(3.3) 



With the excess noise factor F = 1.23 and the energy 22.5 MeV emitted as light, 
which corresponds to 25 MeV energy of the proton with 0.9 quenching factor, one obtains 
that the statistical term in energy resolution is equal a/E = 13.2%. This contribution 
is rather high, however it does not explain completely the measured resolution of 35 %. 
One of the sources of uncertainty is that the PWO-II light yield can differ from crystal to 
crystal and that our sample has inferior properties than assumed as the mean. The rest of 
the difference should be attributed to electronic noise, which could not be estimated with 
sufficient certainty. One should add on this occasion that we had to modify our setup shown 



in Fig. 3.1 (and in Fig. 3.3) when going from APDs with 2 cm^ sensitive area (Advanced 
Photonix Inc.) to APDs with 1 cm^ (Hamamatsu Photonics). The modifications, consisted 
in replacement of the large-sized preamplifier (GSI-CATSA), which was kept outside the 
small chamber and contacted with the APD via a short (~ 10 cm) cable, with a preamp of 
very small size which was installed within the extension shown on top of the small chamber 



(see Fig. 3.2 ). The length of the contact wire was 2 cm in this case. Moreover, this way the 
whole system was enclosed in a Faraday cage formed by the massive walls of the chamber 
and the extension. This reduced the noise pick-up significantly and made it possible to 
perform measurements with an APD of reduced area and a scintillator with low light yield 
(PWO). 

One should stress, however, that our team acquired a lot of practical experience in 
scintillator readout with APDs in the course of experiments at the cyclotron C-30. This 



experience was essential to perform the work described in Sec. |3.3[ The results of Sec. |3.2 
are partly summarized in a conference presentation pT]. 



3 Scintillator and APD test measurements 



83 



3.2.4 Outlook 



The conditions under which the measurements of Sec. |3.2| with 25 MeV protons have been 
performed were rather remote from those that are expected to be met with the PANDA 
EMC. This conclusion emerged from the measurements with a 3 x 3 matrix of full-length 
PWO-II scintillators performed at the MAMI microtron in the J. Gutenberg University 
in Mainz Our team was invited to participate in this experiment by the equipe led 
by Prof. Rainer Novotny from the University of Giessen, the spokesperson of PANDA 
EMC. Quasi-monochromatic gamma-rays were produced by bremsstrahlung of 855 MeV 
electrons in a thin C target with a simultaneous momentum analysis of the electron in a 
magnetic analyser. The latter has a position sensitive detector along its focal plane. By 
selecting a narrow band of electron momenta one obtains a coUimated beam of photons in 
the direction of the primary electron with well defined energy. The matrix was equipped 
with PMT readout and was kept at -25° C in a thermally isolated box whose interior 
was cooled with dry and cool nitrogen gas from a powerful cooler. The 9 scintillators were 
inserted one-by-one on the axis of the collimator to perform energy calibration of the whole 
matrix. This way it was feasible to determine eventwise the distribution of energies of the 
electromagnetic shower in the central crystal and its nearest neighbors when the central 
crystal was irradiated. By adding all these contributions it was possible to restore the line 
shape as a function of the photon energy. The line shape appeared close to Gaussian with 
the following dependence of the dispersion, <7{E), on photon energy: 

m^«J^^OM%, (3.4) 



where energies E are expressed in GeV, so that a{E)/E is dimensionless. 

This experiment was a milestone of the PANDA-EMC development and its results 
constituted a reference point for the development of our own instrumentation. A special 
small-size charge-sensitive preamplifier has been developed by the group of electronic en- 
gineers from the University of Basel [SI] for R&D work on the PANDA EMC. The above 
3x3 matrix was equipped with Hamamatsu's S8664-1010 APDs and the Basel preampli- 
fiers. The experiment was repeated with sixteen tagged photon energies in the range 40.9- 
674.5 MeV, which were used to determine the energy dependence of the response function 
[96j . The matrix was kept at 0°C. The following dependence of a{E)/E on E has been 
obtained in the quoted experiment: 

^ = 1^+0.65%. (3.5) 

E Je 
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3.3 Measurements with low energy 7-rays at the Van-de-Graaff accelerator 



One notes some increase in the statistical term in comparison with the PMT readout 



Eq. 3.4, however overall the resolution is excellent and a forecast for energy resolution 
close to 2 % at 1 GeV obtained. 

It is clear that the experiments with bremsstrahlung 7-rays performed at MAMI leave 
a wide gap of about 0-50 MeV in which measurements are necessary in order to verify 



whether the dependence of a{E)/E on the photon energy suggested by Eq. 3.4 or Eq. 3.5 
is also obeyed at these lower energies. It is necessary, however, that similar conditions are 
preserved: 

• use of monochromatic 7-rays, 

• use of full-size 20 x 20 x 200 mm PWO-II scintillators, 



use of the same reflector in tight contact with the scintillator, achieved with the aid 
of a thermoshr inking plastic tubing, 

use of an Hamamatsu S8664-1010 APD in contact with the scintillator via an optical 
oil of the same or a similar kind, 

use of the preamplifier of the same type, 

application of cooling of the scintillator and an APD, preferably down to — 25°C 



In the following section Sec. |3.3| an apparatus will be described which incorporates the 
above features. 



3.3 Measurements with low energy 7-rays at the Van- 
de-GraafF accelerator 

The measurements of PWO energy resolution with APD readout performed with the 25 
MeV proton beam, which were described above, had several drawbacks. The limited cooling 
capacity of the Peltier elements did not allow to cool scintillators but only APDs. In case 
of PWO its cooling from room temperature down to — 25°C should give fourfold increase 
in light output, that is the reason why a measurement at this temperature presents a lot of 
interest. The second drawback of the performed measurements with the proton beam was 
the distance between the APD and scintillator which had to be left to keep APD thermally 
isolated. Unfortunately, this worsens significantly the light collection which leads to high 
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statistical term in energy resolution, with the measured values worse than expected at the 
given energy. 

The experimental setup has been constructed for measurements of energy resolution 
of the cooled PWO with APD readout to detect low energy 7-rays produced in radiative 
proton capture reactions. The setup was designed to incorporate single crystals of the size 
20 X 20 X 200 mm, which is close to the final shape of the crystals foreseen in the PANDA 
experiment [51] . 



3.3.1 Production of photons in radiative capture reactions 

To determine the possible source of low energy 7-rays the following limitations as to the 
choice of the reaction had to be taken into account: 

1. the beam energy accessible with the Van-de-Graaff generator is limited to 2 MeV for 
protons and deuterons, 

2. the possibility to produce a target with reasonable thickness, 

3. internal spread of 7-ray spectra, 

4. the reaction cross-section, which defines the expected event rate. 

After considering the above mentioned limitations the following reaction was selected 
as the possible source of 7-rays: 

''B{p,^y'C. (3.6) 

This reaction was studied in detail in the 50-ties |119j and the 60-ties [3]. The study of 
excitation curve for this reaction has shown the existence of three resonances at 0.163 MeV, 
0.675 MeV and 1.388 MeV; the latter rather broad with F ^1.15 MeV. The 7-rays are 
emitted in a direct transition to the ground state of ^^C as well as in a cascade transition 
through the 4.44 MeV level of ^^C. For the ground state transition the angular distribution 
is peaked at 90° whereas for the excited state transition it is nearly isotropic for the energy 
range 4 MeV < -Ep < 14 MeV, where the data are available [3] . This favours measurements 



at 90° to the proton beam momentum. Table 3^ presents 7-ray energies for all the three 
resonances. 

For the proton energy Ep = 0.675 MeV the intensity of the transition to the ground 
state is ~ 20 times lower than that to the 4.44 MeV state which yields the correspondingly 
smaller intensity of the 16.6 MeV line in comparison with the 12.2 MeV one. For the 1.388 
MeV resonance the intensities of both lines are comparable at 90° to the beam direction. 
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Ep{MeV) 


E 


^ (MeV) 


0.163 


AAA 


11.6 


16.1 


0.675 


AAA 


12.2 


16.6 


1.388 


AAA 


12.8 


17.2 



Table 3.4: Gamma ray energies at the three resonances in ^^^(p, 7)12(7 reaction at 90° to the 
beam direction. 

It has been found in the performed measurements that our B targets are contaminated 
with ^^F, which produces 7-rays with the energy 6.13 MeV from the reaction ^^F{p, a'j)^^0. 
This reaction has a strong resonance at £'p= 1.374 MeV with cross-section much higher than 
the reactions with ^^B [71] . Subsequently, this reaction was used as the source of calibration 
7-rays with energy 6.13 MeV. 

A complementary reaction to those mentioned above is T{p,'j)'^He, producing a mono- 
chromatic 7-ray line with the energy 20 MeV. However, the low concentration of T in the 
available target did not allow to perform measurements in a reasonable time, therefore all 
the results presented below are related to the ^^B target. 



3.3.2 Concept of the measurements. Monte Carlo simulation of 
the experimental setup. 

The three main processes of 7-ray interaction with matter in the energy range 0-20 MeV 
are: photoelectric effect, Compton scattering and pair production in the nuclear electric 



field. The energy dependence of these processes for lead is presented in Fig. 3.13 The 



dominant process of interaction above some energy is electron-positron pair formation. 



For PbWOi material this energy is around 10 MeV, as follows from Fig. 3.13 and photons 
produced in radiative capture reactions in ^^B lie above this energy. It was assumed that for 
a significant part of events the first interaction of the detected 7 with the PWO scintillator 
is an electron-positron pair creation, however the Compton scattering is not negligible. The 
produced positron annihilates with emission of two 511 keV photons. Both of them can be 
absorbed in the scintillator but there is a high probability that one or even two of them 
escape. The coincidence setup has been proposed with the PWO scintillator surrounded by 
two plastic scintillators, which present two halves of a cylinder with the PWO scintillator 
located along its axis. The requirement of no signal in both plastic scintillators could be 
an indication that all the energy of the primary 7-ray is deposited in PWO, so that the 
recorded spectrum gives its proper resolution at the incident energy. 

The proposed setup might be termed "three crystal pair spectrometer". Such spectrom- 
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10 eV lliEV IMeV 1 GeV 100 GeV 

Photon Energy 

Figure 3.13: Photon total cross section as a function of energy in lead, showing contribution 
of different processes: fjp.e - atomic photoelectric effect, anayieigh - Rayleigh coherent scattering, 
o'Compton - incoherent scattering, Knuc - pair production (nuclear field), Kg - pair production 
(electron field) [4J. 

eters have been used since the sixties for the study of radiative capture nuclear reactions 

m- 

To prove its feasibiUty in the current appHcation detailed Geant 4 Monte-Carlo simu- 
lations have been performed. Simulations were done under the following assumptions: 

• Geometry. 3 active volumes were implemented - PWO crystal and 2 plastic scin- 
tillators. No support structure was taken into account. The PWO crystal was im- 
plemented as a G4Box with dimensions 20 x 2 x 2 cm^ and the 2 plastic scintillators 
were implemented as G4Tubes, i.e. sections of a tube with the internal radius 2.6 
cm, external radius 5.5 cm and 13.4 cm height, covering PWO over almost the full 
azimuthal angle since a gap between the 2 plastics of 2 mm is implemented. The 
PWO crystal is oriented along the OZ axis and its center is at 40 cm from the origin 
of coordinates, at which the source of 7-rays is assumed to be located. 

• Generated events. The 7's with energies 6.13 MeV, 12.8 MeV and 17.2 MeV were 
generated propagating in the OZ direction from a disk of 4 mm radius in the XOY 
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Figure 3.14: Energy distributions in PWO for the 17.2 MeV 7, resulting from the Monte-Carlo 
simulations described in Sec. |3.3.2[ The three panels correspond to: the spectrum of all events 
registered in PWO (upper panel), the spectrum in anti-coincidence with both plastics (central 
panel) and in coincidence with one of the plastic scintillators (lower panel). 



plane with its center at the origin of coordinates. 100 000 events for each energy were 
generated. 

• Physics list. The standard electromagnetic physics list was used, with the default 
cut of 1 mm for the range of e"*", e~ and 7 . 

The energy deposited in the PWO scintillator for the 17.2 MeV 7 is presented in Fig. 



3.14 The top spectrum represents the energy distribution of all recorded events. The 
middle and bottom spectra correspond to anti-coincidence and coincidence with one of the 
plastic scintillators. The three peaks can be seen in the spectra corresponding to the total 
energy, the one-7 escape and the double escape of the positron annihilation photons. In 
the upper spectrum the single escape peak is the dominant one. It can be also seen from 
the central spectrum that the requirement of no signal in any of the plastic scintillators 
increases the significance of the full energy peak but does not suppress the single- and 
double-escape peaks completely. The requirement of a coincidence (lower panel) removed 
the full-energy peak and enhanced the continuous background of Compton scattering in 
PWO. 



In Fig. 3.15 the energy deposited in plastic scintillators is presented. The peak which 
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I Energy distribuUon in plastic scyntillator No 2, E=17.2 MeV | 




Figure 3.15: Energy distributions in the two plastic scintillators for the 17.2 MeV 7-rays detected 



in PWO. The spectra result from the Monte-Carlo simulations described in Sec. 3.3.2 



is seen in both spectra corresponds to the Compton edge of the 0.511 MeV 7. According 
to the formula 

Emax = E ■ o , ^ 771 5 (3-7) 

the Compton edge of the 0.511 MeV 7 is at 0.341 MeV. The reader should be aware that 
the plastics, containing only H and C, respond almost exclusively via Compton scattering 
to low energy 7-rays. 



In Fig. 3.16 the energy distribution in PWO for all the 3 analysed energies (6.13, 12.8 



and 17.2 MeV) is presented folded with the energy resolution of a/E = 5%. The anti- 
coincidence with plastic scintillators is assumed. This plot gives an orientation as to the 
shape of energy spectrum expected in our measurements, however the actual values of 



energy resolution will appear to be worse than assumed in Fig. 3.16 (see Sec. 3.3.4). 



3.3.3 Experimental setup 
3.3.3.1 Van de Graaff accelerator 



The low energy 7-rays for the present study were produced in radiative proton capture 
reactions with protons accelerated in a Van de Graaff accelerator belonging to SINS and 
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Figure 3.16: Energy distribution in PWO for the 6.13, 12.8 and 17.2 MeV 7's having equal 



intensities, detected with the energy resolution a/E = 5% (see Sec. 3.3.2). 



the Warsaw University. This 3.5 MeV Van de Graff machine provides p, d, ^He and ^He 
beams with energy resolution of about 1 keV. The total accessible p and d beam currents 
on a target are up to 50 and for ^He and ^ife ions up to about 30 ^A. During 
our measurements a typical beam current on target was about 10 ^A. A photo of the 



accelerator located inside the closed pressure tank is presented in Fig. |3.17[ The beam 
is directed vertically into the target room lying underneath, protected with about 1 m of 
concrete. In the target room it is bent by 90° with a magnetic analyser and then it may 
be redirected in the horizontal plane with the aid of a distributing magnet into one of the 
three beam pipes terminated with an experimental set-up. One of them was used for the 
present scintillator studies, the other two are currently used in solid-state and chemical 
analysis studies, the latter with the PIXE method. Tools for steering, focusing and visual 
observation are provided along the path of the beam. 



3.3.3.2 Setup for the measurements 



The vertical section of the experimental setup is presented in Fig. 3.18 The common 
vacuum jacket houses two volumes separated from the ambient. The lower volume of 
the setup houses a PWO [PhWOi) scintillator, an APD (Hammamatsu S8664-1010), a 
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Figure 3.17: The SINS/UW Van-de-Graaff accelerator inside the pressure tank used to accel- 
erate protons for production of low energy 7-rays in radiative capture reactions. In the top part 
a crane is shown used to lift the top section for ion-source installation and/or repair. 
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preamplifier for APD signals and a Pt-100 sensor to measure temperature close to the 
APD. The PWO-II scintillator of the size 20 x 20 x 200 mm, produced by BTCP, was used 
in the measurements. 

The PWO crystal is surrounded by two cylindrical plastic scintillators (EJ-200) with 
PMT (XP4312) readout for coincidence measurements. The length of the plastic scintil- 
lators along the vertical axis is 134 mm, the internal and external radii of the cylinder 
are 26 and 55 mm, respectively. The EJ-200 material is characterised by high light out- 
put - 64 % of anthracene, long attenuation length and the emission spectra matching well 



the used photomultipliers. A drawing of the plastic scintillators is presented in Fig. 3.19 



Fig. 3.20 illustrates the way they are fastened to the aluminum tube, which is a part of 
the outer jacket of the cryostat. The lead collimator is placed between the target and the 
scintillator with a hole of 1.5 cm. The collimator prevents from direct photon hits of the 
plastic scintillators. On the other hand a rather precise relative positioning of the target, 
collimator and the axis of the PWO crystal is required. 



3.3.3.3 Cooling 

Through the upper volume of the cryostat the coolant is circulated from the JULABO F32- 
ME Refrigerated/Heating Circulator. The circulator is suitable for internal and external 
temperature application in the range of temperatures from — 35°C to +200°C. The nominal 
cooling capacity is 150 W at — 20°C and 30 W at — 35°C, i.e. decreases rapidly with the 
decreasing temperature. Two types of coolant have been tested: JULABO "Thermal" bath 
fluid H5S and Galden heat transfer fluid HT-70 of Solvay Solexis. The results presented 
below had been taken at -20° C which was achieved with the H5S coolant. It is slightly 
higher than — 25°C planned for the PANDA EMC environment. The HT-70 coolant was 
also tested in an attempt to reach lower temperatures. It is typifled by the lower viscosity 
of 1.09 mm^ ■ vs 4 mm^ ■ for the Julabo H5S coolant, so that it permits a higher 
velocity of the circulating liquid. However no signiflcant improvement has been achieved 
with the difference at the level of 0.5°C. We conclude that with the given thermal losses 
in our system we need a more powerful refrigerator to reach the desired -25°C. 

The cooled elements are isolated from ambient by vacuum common with the accelerator 
vacuum system, therefore the minimum achieved temperature depends on the vacuum level. 
The lowest temperature was achieved after a prolongued pumping when vacuum stabilised 



at 6- 10 Torr. The lower volume in Fig. 3.18 is fllled with He at slightly above atmospheric 



pressure which serves as a contact medium between the cooled bottom copper plate of the 
upper volume and the scintillator, APD and the preamplifler located in the lower volume. 



3 Scintillator and APD test measurements 



93 



COOLAMT IN 



COOLANT OUT 



He FILUNG 



POWER SUPPLY 
FEEDTHROUGH 




TO ACCELERATOR 
VACUUM 

PT100 

THERMOMETER 



Cu GASKET 
PREAMPLFIER 



POWER SUPPLY 
FEEDTHROUGH 



THIN WLLED SST 

PWO SCINTILUTOR 
20x20x200 mm 



LEAD 

SHIELDING 



WATER COOLING 



Figure 3.18: Vertical section of the setup used for the measurements of energy resolution of 
cooled PWO scintillators with low energy 7-rays from the radiative capture reactions. 
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3.3 Measurements with low energy 7-rays at the Van-de-Graaff accelerator 




Figure 3.19: A drawing of the plastic scintiUators, which are used for coincidence measurements. 
EJ-200 is machined to form a cylinder cut into two halves along the vertical axis. A cylindrical 
plexiglass light-guide is glued at the center of each of the two half-cylinders to maximize light 
transmission to the photocathode of a 3 diameter XP-4312. The PMs are protected with /i-metal 
shields from stray electromagnetic fields. Transistorized voltage dividers are seen at the two PM 
extremes. 

Temperature in the lower volume was monitored with a Pt-100 sensor, whose resistance 
was measured with a Keithley digital multimeter, which provides the possibility to store 
the temperature information in a PC via an USB port. 

The working temperature of -20° C was achieved on the average after 6 hours of cooling. 



3.3.3.4 Targets 

The high melting point of boron 2300°C) makes it difficult to prepare thin boron 
targets of uniform thickness by an evaporation method. The targets were prepared by 
evaporating natural boron on Mo backings 0.5 mm thick, using electron bombardment 
of a boron pellet in a Ta boat, as described in [56]. Several targets have been produced 
with thicknesses ranging from 2 fim, corresponding to ~100 keV energy loss for 1 MeV 
protons, to 10 fim. The target holder is constructed so, that the back surface of the Mo 



backing is water cooled (see the bottom part of Fig. 3.18). It is tilted at 45°, while the 
axis of the Pb collimator is at 90° to the beam direction. The quality of the target has 
been studied with a large Nal(Tl) scintillator 04" x 4" coupled to a 4" diameter PMT. The 
measurements at the beam energy Ep = 0.67MeV have been performed. The two high 
energy lines originating from the transitions from the capturing state to the ground state 



and the 4.44 MeV ffist excited state in ^^C are clearly seen in the spectrum of Fig. 3.21 
However, the peaks corresponding to the transition 4A4MeV — )■ g.s. [ground state) , which 
would be also useful for the PWO response studies, are obscured by the ^^F contaminant. 
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Figure 3.20: A photograph of the setup demonstrating the fastening of plastic scintihators. The 
cryostat is shifted out of the heavy Pb shielding. A beam-line with a gate-valve and a small- 
diameter stainless-steel tube, which terminates with a water-cooled target are seen below the 
cryostat. 
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Figure 3.21: Energy spectrum of 7-rays produced in a 2 fim thick boron target at Ep = 0.67 



MeV, detected with the (/>4 x 4 NaI{Tl) detector mentioned in Sec. 3.3.3.4 



which has high chemical affinity to B. It should be stressed that a)^^0*(6.13 MeV) 

proceeds with a "normal nuclear" cross-section of the order of m6, while 7)^^C is 

mediated by an electromagnetic interaction, therefore its intensity is relatively suppressed. 



3.3.3.5 Readout electronics 

The readout chain of PWO starts with a preamplifier, which has been designed and pro- 
duced at the Basel University especially for PANDA EMC design studies [U]. It is the 
Low Noise / Low Power (LNP) charge preamplifier with charge sensitivity 0.5 V/pC. The 
LNP preamplifier fits detector capacitances in the range from 50 pF to 300 pF, therefore 
it is well applicable to work with (10 x 10 mm?) Large Area APDs S8664-1010 produced 
by Hamamatsu possessing the capacitance of about 270 pF. The measured noise in combi- 
nation with an APD and the postamplifier with the shaping time of 750 ns is 1900 crms 
at — 25°C [SI]. The signal from the preamplifier goes into the standard spectroscopic am- 
plifier ORTEC 452 with the shaping time of 1 /is. The amplified signal is sent in parallel 
to a Linear Gate and a Single Channel Analiser (SCA) with a threshold set to suppress 
the electronic noise, whose output triggers the linear gate. For both plastic scintillators 
the signals from PMTs are amplified with ORTEC amplifiers with the same shaping time 
of 1 ^s. The three output signals are fed into a 4-channel Flash ADC CAEN N1728B 
(100 MHz, 14 bits). CAEN N1728B Flash ADC has a 10 ms wide buffer which is filled 
continuously with digitized amplitudes from each of the three channels. In the acquisition 
software for the N1728B FADC the signal coming from PWO is assigned as a trigger, i.e. 
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Figure 3.22: A scheme of readout electronics used for the measurements of energy resolution of 
PWO scintillators with low energy 7-rays. 

when it exceeds some threshold the inputs in all the three channels are stored from the 
FADC buffer 10 fis before and 20 fis after the trigger into hard drive of a PC. A slightly 



simplified scheme of the readout is presented in the Fig. 3.22 



3.3.4 Results 

The measurements of energy resolution of PWO scintillators with Large Area APD readout 
have been performed at the proton beam energy of 1.4 MeV. It is close to the 1.390 MeV 



resonance in 7) C reaction at which (see Sec. 3.3.1) 7-rays with energies 12.8 MeV 



and 17.2 MeV should be produced with comparable cross-sections. 

The cryostat was cooled during 20 hours before the measurements were started to 
stabilize the temperature at — 20°C. The bias voltage on the APD was set at +377 V and 
its inverse dark current was continuously monitored at an average level of 25 nA. The data 
have been collected for 10 hours. The collected energy spectrum is plotted in Fig. |3.23| in 



the logarithmic scale. Three energy peaks are seen: the intense peak is assigned to 6.13 
MeV from the 0:7)^^0 reaction and the two weaker maxima to 12.8 MeV and 17.2 

MeV lines from 7)^^(7, following our previous experiment with the Nal(Tl) shown 



in Fig. 3.21 



Fig. 3.24 presents the energy spectra measured with the plastic scintillators. The edge 
of the Compton continuum, corresponding to the 0.511 MeV annihilation 7-ray, is clearly 
seen in both spectra. A single ORTEC power supply provided HV for both PMTs, so 
that the differences in amplitude of the Compton edge reflect the asymmetry of PMT 
gains at -1.1 kV of the common voltage. Vertical lines indicate the range of energies of 
annihilation quanta used in coincidence or anticoincidence requirements. A signal within 
the shown range was considered as the detection of a 0.511 MeV annihilation 7-ray, whereas 
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Figure 3.23: Spectrum of B{p,^y^C and ^^F{p,ay^O (target contaminant) gamma-rays 
emitted at Ep=1.39 MeV, detected with the 20 x 20 x 200 mm PWO scintillator cooled down to 
-20° C (logarithmic scale). 



its absence in both plastic scintillators served as an indication that the entire energy was 
deposited in the PWO scintillator. It was found that the anticoinidence requirement tends 
to improve the widths of the observed peaks and shifts them to higher energy since without 
this condition peak positions correspond to large extent to a single escape peak (see also 
Sec.[3X2|. 



Fig. 3.25 presents for convenience the energy spectrum of 7-rays in the linear scale. The 
left and the right parts are scaled differently. The presented spectrum corresponds to the 
anticoincidence requirement with the plastic scintillators. A fit with gaussians has been 
performed to extract energy resolutions. The presence of transitions to the higher states 
of ^^O with energies 6.92 MeV and 7.12 MeV was taken into account in the fit. According 
to [H] the intensities of these two lines are at the level 3-8% of the 6.13 MeV transition 
depending on the proton beam energy. For the fit the intensities of each of them were fixed 
at the level of 5% of the transition 6.13MeV — )■ g.s.. Moreover, due to the asymmetry of 
the 6.13 MeV peak the fit was performed only over its right wing. The fitted Gaussian 
energy resolution is 19 ± 3% for the 6.13 MeV, 15 ± 2% for the 12.8 MeV and 9.6 ± 1.4% 
for the 17.2 MeV peak. 

A comparison with the results of PWO resolution measurements using bremsstrahlung 
photons in the 40.9 - 675 MeV range [96] at MAMI is of interest. The measurements at 
MAMI were performed with a 3 x 3 matrix of 20 x 20 x 200 mm^ PWO scintillators with 
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Figure 3.24: Energy spectra in the two plastic scintillators. The range of amplitudes between 
the vertical lines is used for the anticoincidence requirement. 
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Figure 3.25: Spectrum of B{p,j) C and F{p,a) O (target contaminant) gamma-rays 
emitted at Ep=1.39 MeV, detected with the 20 x 20 x 200 mm PWO scintihator cooled to -20°C 
(linear scale, note that the two parts of the plot have different scales) in anticoincidence with the 
plastic scintillators. The solid line is a fit to the measured spectrum using three Gaussians with 
the indicated relative dispersions. 



100 



3.3 Measurements with low energy 7-rays at the Van-de-Graaff accelerator 




Figure 3.26: Measured gaussian resolutions (solid points with error bars) compared with the 
solid curve parameterizing the results obtained in the 40.9-675 MeV range at MAMI [96]. 



APD readout cooled down to 0°C APD and preamplifiers of the same type as in our work 
have been used in [96]. The measured energy resolution is well described by the following 
dependence on the photon energy (in GeV): a/E = l.^QVo/ E{GeV) + 0.65%. It is 
estimated by the authors of [HE] that lowering the temperature from 0°C to -25°C should 
increase light output by the factor of 2 thus the statistical term in energy resolution will 
be reduced by the factor v^, so that the resolution to be compared with our data reads: 

a/E = 1.31%/ ^/E (GeV) + 0.65%. (3.8) 



Fig. 3.26 compares our experimental points and the extrapolation of the results from MAMI 



corrected for the difference in working temperatures (solid curve). One may conclude that 



our results prove that Eq. 3.8 is a valid extrapolation of the high energy data into the 



lower energy range, where resolution has not yet been measured. The surplus of counts 



above the solid line seen in Fig. 3.25 between the 6.13 and 12.8 MeV peaks suggests that 
an assumption of pure gaussian line shapes is an oversimplification of the reality. However, 
the measured spectrum does not leave enough space to test more complex hypotheses (e.g. 
Gaussians with tails). Assuming such an eventuality would make our present estimates 
lower limits of the true resolution. The proper testing ground might be offered by the 
T{p, '^YHe reaction, which is expected to demonstrate the true line shape in the space 
between the 20 MeV and the 6.13 MeV peaks, the latter ever present in the measured 
spectrum. The tritium target needed for such an experiment is, however, an expensive 
investment. 



Chapter 4 



Monte-Carlo simulation of PANDA 
physics performance 



This part of the PhD Thesis is devoted to demonstration of the PANDA detector physical 
performance in detection of the he state in charmonium. The importance of a study of 
he, i.e. of the measurement of its mass, width and precise determination of its quantum 



numbers has been justified in Sec. 1.3, The first part of this chapter (Sec. 4.1) is devoted to 



the description of software used in the analysis and the second part (Sec. 4.2 4.8) provides 
the results of simulations. 



4.1 Panda offline software. 

The offline software developed for the PANDA physics studies has been used in this analysis 
and the following description is based on the Physics Performance Report for PANDA [52] . 
The software consists of several tools adapted from other HEP experiments which include 
event generators (EvtGen, DPM), particle tracking with GEANT4, digitization of the 
signal, reconstruction and identification of charged and neutral particles and high level 
analysis tools providing vertex and kinematical fits. 

4.1.1 Event generators. 

Two event generators have been used in this analysis. EvtGen is used to produce the 
reaction of interest with the possibility to define decay modes of the particle and angular 
distributions. Dual Parton Model (DPM) can be used to produce generic background from 
pp annihilations. 
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4.1.1.1 EvtGen 

EvtGen event generator [83j was originally developed for the Babar experiment at SLAG 
to study the physics related to B-mesons. The modular design of EvtGen allows to extend 
it by incorporating arbitrary physics channels and decay models, for these reasons it was 
adapted by the ATLAS and PANDA experiments. Decay amplitudes, instead of probabil- 
ities, are used for the simulation of decays. The amplitude for each node in a decay tree 
is used to simulate the entire decay chain, including all angular and time-dependent cor- 
relations. The implementation of each decay amplitude is independent of how the mother 
particle was generated or how the daughter particles are going to decay. EvtGen computes 
the spin density matrices and the decay probability and decay models, which implement 
a single node in a decay tree, must only specify the decay amplitude for each combination 
of the mother and daughter spin states. 

4.1.1.2 DPM 

In high energy proton-antiproton collisions soft processes dominate in multiparticle produc- 
tion and perturbative QGD is not applicable for their description. However, suitable large 
N expansions of QGD provide a topological classification of diagrams and a potentially 
useful non-perturbative approach. This topological expansion, when supplemented with 
generally accepted theoretical principles like duality, unitarity, Regge behavior and the 
parton structure of hadrons, provides the basis underlying the dual parton model (DPM) 
[39] . It has been shown that DPM provides a complete phenomenological description of 
all facets of soft processes. 

The main objects of the model are constituent quarks having masses ~ 300-350 MeV, 
strings, and string junctions for baryons. It is assumed that mesons consist of a quark and 
an anti-quark which are coupled by color forces. The vortex lines of the field are concen- 
trated in a small space region forming a string-like configuration. So, mesons are considered 
as strings with small masses. Various processes are considered in baryon-antibaryon inter- 
actions, such as string junction annihilation and creation of three strings, quark-antiquark 
annihilation and string creation between diquark and anti-diquark, quark-antiquark and 
string junctions annihilation and so on. The energy dependence of cross sections for these 
processes is parametrized to reach an agreement with experimental data. The string frag- 
mentation is performed by the mechanism similar to that included in the LUND model [5]; 
after the string fragmentation all unstable hadrons decay with the aid of the code DEGAY 
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4.1.2 Particle Tracking and Detector Simulation. 

The detector simulation starts with propagation of the generated particles through the 
PANDA detector using GEANT4 transport code. It takes into account all the possible 
processes of particle interactions and decays. As an output GEANT produces the collec- 
tion of hits, which contains information on the intersection points and energy lost by all the 
particles in corresponding parts of the detector. For the description of the detector geome- 
try the approach developed by the CMS experiment is followed. The Detector Description 
Database (DDD) is used in which the detector is represented as a multigraph structure with 
a compactified description. The Detector Description Language (DDL) based on XML is 
used to encode the geometrical properties of the components, their relative positions and 
the materials. The simulations have been performed with the complete detector setup and 
the Straw Tube Tracker (STT) chosen as the central tracker. The pellet target scenario 
was considered in order to take into account the material budget presented by the target 
pipe. The Gaussian spread of interaction points typified by a— 0.275 mm was taken into 
account in simulations. 

4.1.2.1 Digitization 

The output of Geant simulation is used as an input for the digitization step, which models 
the processing of the signal from individual subdetectors by their front-end electronics. The 
output of the digitization should correspond to the detector response in the real experiment. 

Readout of MVD. MVD will use two different sihcon detector types, sihcon strip and 
pixel detectors. The readout of the sihcon devices is for both types different and is treated 
differently in the digitization scheme. The hit position on the sensor surface defines the 
channel number and the deposited energy defines the charge collected by the electronics. 
Strip sensors will be sensitive on both sides and the formation of digitized channels is done 
independently on both sides of the sensor. In the case of pixel detectors, the trajectory is 
projected to the surface and depending on its relative orientation, all excited pixel cells are 
calculated and the charge signal is shared among all pixel cells depending on the fraction 
of the local track. 

Straw Tube Tracker and Drift Chambers. Digitization for the Straw Tube Tracker 
(STT) and the Drift Chambers (DCH) has been treated in similar way. Both devices consist 
of wires inside an Ar — CO2 gas mixture volume. If a charged particle traverses this gas 
volume, the local helix trajectory is derived from the corresponding GEANT4 intersection 
points. The drift time of the ionisation electrons is estimated from the smallest distance of 
this helix trajectory to the wire dpoca ipoca - point of closest approach). The uncertainty of 
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tfie drift time is taken into account by smearing dpoca witfi a Gaussian distribution witfi a 
standard deviation of a— 150 jim for the STT, and a — 200 jim for the DCH devices. The 
energy signal of a straw tube is finally calculated by taking into account Poisson statistics. 

GEM readout. Each GEM station consists of two detection planes. The distance 
between the detection planes is 1 cm. It has been assumed that each detection plane has 
two strip detection layers with strips having perpendicular orientation to each other. The 
gas amplification process and the response of the strip detector has not been simulated in 
detail. Instead, the entry point of a charged track into the detector plane has been taken 
directly from GEANT4 and smeared with a Gaussian distribution of 70 /xm width in each 
strip orientation direction. 

DIRC detectors. The light propagation in the Cherenkov radiators, the signal pro- 
cessing in the front-end-electronics, and the reconstruction of the Cherenkov angle have 
been modeled in a single effective step. The resolution of the reconstructed Cherenkov an- 
gle (7(7 is mainly driven by the uncertainty of the single photon angle a^-y and the statistics 
of the relatively small number of the detected Cherenkov photons, Nph. A single photon 
resolution of crc,7 = 10 mrad was used, corresponding to the experience with the existing 
DIRC detectors. The number of detected photons was calculated from the velocity, /3, of 
charged particles passing through the quartz radiators and the path length, L, within the 
radiator. The sensitive wavelength interval, [A„(j„; \rnax\: was chosen as [280nm; 350 nm], 
and a total efficiency oi e = 7.5 % was used to take into account the transmission and 
reflectivity losses as well as the quantum efficiency of the photo detectors. The precision of 
the measured Cherenkov angle obtained with the digitization and reconstruction procedure 
has a gaussian shape with (T=2.33 mrad. 

EMC readout. The parameters of the EMC which were taken into account for digi- 
tization in the target spectrometer are electronics noise and photon statistics. According 
to the prototypes of PANDA EMC front-end electronics the noise level at — 25°C is at 
the level of 0.5-1 MeV. The Gaussian distribution of the noise with 1 MeV width was 
used in simulations as a conservative estimate. For the PWO-II scintillator around 80 e-h 
pairs/McV arc expected with the LAAPD readout. Additional contribution to photon 
statistics is given by the excess noise factor of an APD, which was measured for Hama- 
matsu APDs at 1.38 for the gain M = 50. For the shashlyk calorimeter in the Forward 
Spectrometer, consisting of Icad-scintillator sandwiches only a fraction of roughly 30 % of 
the incident energy is deposited in the scintillator material. Based on this energy deposit 
the electronics noise with a= 3MeV has been considered for digitization. 
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4.1.3 Reconstruction 

4.1.3.1 Charged particle track reconstruction 

MVD cluster reconstruction. The MVD provides very precise space point measure- 
ments as a basis for the track and vertex reconstruction. The hit resolution of individual 
MVD measurements has a = 21.0 fim for hits where only one pixel contributes to the hit 
cluster and a = 3.8 fim for multihit clusters, where the charge weighting between pixel 
cells in the cluster can be used to calculate the mean position of the hit. 

Global track reconstruction. The track object provides information about a charged 
particle path through space. It contains a collection of hits in the individual tracking 
subdetectors. An idealized pattern recognition has been used for track building based on 
Monte Carlo information to assign reconstructed hits to their original tracks. The tracks 
in the target spectrometer are fitted with the Kalman Filter algorithm, which considers 
not only the measurements and their corresponding resolutions but also the effect of the 
interaction with the detector material, i.e. multiple scattering and energy loss. A typical 
choice for the parametrisation of the track in a solenoidal field is a five-parameter helix 
along the principle field direction z. The task of the track-fitting algorithm is to determine 
the optimal parameter vector and its covariance matrix as a function of the flight length / 
in order to create a representation of the track as a piecewise helix. In the physics analysis 
the particle position and momentum can then be accessed through this piecewise helix 
representation. 

Tracking performance. With the used reconstruction software more than 90 % track 
reconstruction efficiency has been achieved for the transverse momenta pt > 0.2GeV/c 
dropping down to 70% for pt = 0.1 GeV/c. The momentum uncertainty can be charac- 
terised with ap/p = 1 % for the pions of 1 GeV/c momentum at a polar angle 20° and the 
achieved vertex space resolution is cr^o = 47.5 fj,m and a^o = 49.2 fj,m for the 3 GeV/c pions. 

4.1.3.2 Photon reconstruction 

Most of the EMC reconstruction code is based on algorithms developed and applied by 
the Babar experiment p,l6]. The first step in EMC reconstruction procedure is a cluster 
reconstruction. Electromagnetic shower during its development deposes energy in a set of 
neighboring crystals which is called cluster. The energy deposits in the crystals, which form 
the cluster, and their positions allow to determine the energy and direction of the initial 
photon. Cluster finding algorithm starts from the crystal with largest energy deposit. Its 
neighbors are then added to the cluster if the energy deposit is above a certain threshold 
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Exti- Tfiis tfiresfiold, on tfie one liand, sfiould be as low as possible, to detect low energy 
photons, on the other hand, it must be sufficiently high to separate the reconstructed 
photons from the readout noise. Ej-u was set to 3 MeV, which corresponds to the energy 
equivalent to the 3a of the expected noise of the electronics. The procedure of adding new 
elements to the cluster is continued for the neighbors until no further module fulfils the 
threshold criteria. A cluster is accepted if the total energy deposit in all of its elements is 
above a second threshold E^, which was set to 10 MeV. 

The next step in EMC reconstruction procedure is to search for bumps within each 
reconstructed cluster, which can appear if the cluster is formed by more than one particle 
detected at small relative angular distances. In this case the cluster has to be subdivided 
into regions which can be associated with the individual particles. This procedure is called 
bump splitting. The bump is defined by two criteria. First is that the energy deposit 
in one scintillator module, Eiocai, must be above the threshold, Emax, while all neighbor 
modules have smaller energies. Emax = 20 MeV was chosen for the Target Spectrometer. 
In addition the highest energy E^max of any of the N neighboring modules must fulfill the 
requirement 0.5(A^ — 2.5) > Ej^max/Eiocai- The total cluster energy is then shared between 
the bumps, taking into account the shower shape of the cluster. For this step an iterative 
algorithm is used, which assigns a weight Wi to each scintillator module, so that the bump 
energy is defined as E^, = WiEi. E^ represents the energy deposit in the i*'* module and 
the sum extends over all modules within the cluster. The crystal weight for each bump is 
calculated with: 

^ Eiexp{-2.bri/rm) , ^ 

E,-^.exp(-2.5r,/r^)' ^ ' ^ 

where - Moliere radius of the scintillator, r^, Vj - distances of the i*^ and module to 
the center of the bump. The procedure is iterated until convergence. 

The spatial position of a bump is calculated via a center-of-gravity method. The radial 
energy distribution, originating from a photon, decreases approximately exponentially. 
Therefore, a logarithmic weighting with 

Wi = max(0, A{Ek) + ln{Ei/Eb)) (4.2) 

was chosen, where only modules with positive weights are used. The energy dependent 
factor A{Ei,) varies between 2.1 for the lowest and 3.6 for the highest photon energies. 

In case of the shashlyk calorimeter in the Forward Spectrometer the same reconstruction 
algorithms are used however with different thresholds: E^ti — 8 MeV, Ed — 15 MeV and 
= lOMeV. 
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The sum of energies deposited in the scintillator material of the calorimeters is in general 
less than the energy of the incident photon. While only a few percent is lost in the TS 
EMC, which mainly originates from energy losses in the material between the individual 
crystals, a fraction of roughly 70 % of the energy is deposited in the absorber material 
for the shashlyk calorimeter. The reconstructed energy of the photon in the TS EMC is 
expressed as a product of the measured total energy deposit and a correction function which 
depends logarithmically on the energy and - due to the layout - also on the polar angle. 
The correction function f(\nE,6) has been obtained from the Monte Carlo simulations 
with single photons. In case of shashlyk calorimeter in FS a correction has been considered 
depending only on the collected energy. 

4.1.4 Charged particle identification (PID) 

In the PANDA experiment several subdctcctors will provide useful PID information for 
specific particle species and momenta. Energy loss measurements within the trackers pro- 
vide good criteria for the distinction between the different particle types below 1 GeV/c; 
the DIRC detector is the most suitable device for the identification of particles with mo- 
menta above the Cherenkov threshold; EMC, in combination with the tracking detectors, 
is the most powerful detector for an efficient and clean electron identification; the Muon 
detector is designed for the separation of muons from other particle species. The best 
PID performance, however, can be obtained by taking into consideration all the available 
information from all subdetectors. 

The PID software is divided in two parts. In the first stage the recognition is done for 
each detector individually, which results in probabilities for all five particle hypotheses (e, 
fj,, TT, K and p). In the second stage the global PID combines this information by applying 
a standard hkehhood method. 

4.1.4.1 Subdetector PID 

dE/dx meeisurements. The energy loss of particles in thin layers of material directly 
provides an access to the dE/dx. As can be seen directly from the Bethe-Bloch formula, 
a given momentum particles of different types have different specific energy losses, dE/dx. 
This property can be used for particle identification. The method however suffers from 
two limitations. First of all, at the crossing points, there is no possibility to disentangle 
particles. Secondly, the distribution of the specific energy loss displays a long tail which 
constitutes a limitation to the separation, especially when large differences exist between 
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different particle yields. In PANDA, two detectors will give access to a dE/dx measurement, 
the MVD detector setup, and the central spectrometer tracking system. 

PID with DIRC. Charged tracks are considered if they can be associated with the 
production of Cherenkov light in the DIRC detector. Based on the reconstructed momen- 
tum, the reconstructed path length of the particle in the quartz radiator, and the particle 
type hypothesis the expected Cherenkov angle and its error are estimated. Compared with 
the measured Cherenkov angle the likelihood and significance level for each particle species 
are calculated. 

Electron identification with the EMC. While muons and hadrons in general loose 
only a certain fraction of their kinetic energy by ionisation processes, electrons deposit their 
complete energy in an electromagnetic shower. The ratio of the measured energy deposit in 
the calorimeter to the reconstructed track momentum (E/p) will be approximately unity. 
Furthermore, the shower shape of a cluster is helpful to distinguish between electrons, 
muons and hadrons. Since the chosen size of the scintillator modules corresponds to the 
Moliere radius of the material, the largest fraction of an electromagnetic shower originating 
from an electron is contained in just a few modules. Instead, a hadronic shower with a 
similar energy deposit is less concentrated. These differences arc reflected in the shower 
shape of the cluster. A neural network based algorithm was implemented for electron 
recognition based on different parameters characterizing cluster shape. 

PID with the muon detector. The particle ID for the muon detector is based on 
propagation of the charged particles from the tracking volume outward through the neigh- 
boring detectors hke DIRC and EMC and finally through the iron of the solenoid yoke, 
where the detection layers of the muon device are located. Then the extrapolated inter- 
section points with the detection cells are compared with the detected muon hit positions. 
In case the distance between the expected and the detected hit is smaller than 12 cm, 
according to Aa of the corresponding distribution, the muon hit will be associated with 
the corresponding charged track. The procedure results in a good muon identification for 
momenta above approximately 1 GeV/c. While a complete muon-electron separation is to 
be expected, a contamination of a few percent by pions can at best be achieved. 

4.1.4.2 Global PID 

The global PID, which combines the relevant information of all subdetectors associated 
with one track, has been accomplished with the aid of a standard likelihood method. Based 
on the likelihoods obtained from each of the individual subdetectors, the probability for 
a track originating from a specific particle type p{k) is evaluated from the likelihoods as 
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follows: 



p{k) 



(4.3) 



where the product with index i runs over all the considered subdetectors and the sum with 
index j over the five particle types e, /i, vr, K and p. 

Due to the variety of requirements imposed by the different characteristics of the bench- 
mark channels various kinds of particle candidate lists depending on different selection 
criteria on the global likelihood are provided for the analysis. The usage of the so-called 
VeryLoose and Loose candidate lists allows to achieve good efficiencies, and the Tight and 
VeryTight lists are optimized to obtain a good purity with efficient background rejection. 

4.1.5 Physics analysis tools 

The analysis user has the choice to reconstruct decay trees, perform geometrical and kine- 
matical fits, and to refine the event selection by using Beta, BetaTools and the fitters 
provided by the analysis software [70] directly in an application framework module, or by 
defining the analysis in a more abstract way using SimpleComposition tools. This TCL 
based high level analysis tool package provides an easy-to-learn user interface for the defi- 
nition of an analysis task and the production of n-tuples, and it allows to set up analysis 
jobs without the need to compile any code. SimpleComposition as well as Beta, Beta- 
Tools and the geometric and kinematic fitters were taken over as well tested packages from 
BaBar, and adapted and slightly extended for the usage in PANDA. The last n-tuple based 
analysis steps were carried out utilizing the ROOT toolbox |36], which provides powerful, 
interactively usable instruments among others for cutting, histogramming, and fitting. 



For the exclusive benchmark channels it turned out that especially a 4C-fit of the re- 
constructed decay tree is a powerful tool to improve the data quality and to suppress 
background. Kinematic fits allow to test if the measured final state is compatible with 
the hypothesis of the particular assumed reaction kinematics. 4C-fit uses the fact that the 
four-momentum of the initial pp state is precisely known and should be equal to the sum 
of 4-momenta of the final state particles. Least-square fit is performed for the kinematic 
variables of the measured final state with 4 constraints: 



4.1.5.1 4C-fit 



(4.4) 



k 
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Decay mode 


BR 


K+K-TT + TT- 

i^*(892)F*(892) 

77 


(1.9 ±0.4) ■ 10-2 
(1.9 ±0.4) ■ 10-2 
(1.5 ±0.6) ■ 10-2 
(1.2 ±0.3) ■ 10-2 
(9.2 ±3.4) ■ 10-3 
(2.7 ±0.9) ■ 10-3 
(2.4 ± 1.1) ■ 10-^ 



Table 4.1: Decay modes of tJc with the corresponding branching ratios [1] 



Et = J2E>' = H \l\tk? + ml (4.5) 

k k 

The constraints are taken into account based on the method of Lagrange multiphers. The 
fit gives for each event the probabihty that the considered final state is compatible with 
the beam energy-momentum. 4C-fits and cuts on the results can also be defined with the 
SimpleComposition tool. 



4.2 he detection. General remarks 

According to theoretical predictions and previous experimental observations |Ul llllj . one of 
the most promising decay modes for observation of the h^. is its electromagnetic transition 
to the ground charmonium state: 

hc^Vc + 7, (4.6) 

where the energy of the photon is Ej = 503 MeV. The rjc can be detected through many 
exclusive decay channels, neutral {rjc — )■ 77) or hadronic. Several selected decay modes of 



the rjc with the corresponding branching ratios, according to [1], are listed in Table 4.1 

The previous observation of he was done in neutral decay mode of rjc in pp annihilation 
[6] or in hadronic decay modes of rjc formed along the chain e~^e~ — t- ^(25*) — ?■ ir^hc, h^ — ?■ 
Vcl |lllj - In the case of PANDA experiment, where h^ is produced as a resonance in pp 
annihilation, the detector is capable to study hadronic final states, hence, the he can be 
observed exclusively whether in one of the hadronic or neutral {rjc — )■ 77) decay modes of 



T]c (see Table 4.1 ). 

The expected signal to background ratio is one of the criteria of feasibility to detect 
he experimentally. For hadronic final states the evaluation of signal to background ratio 
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is even more important due to the fact that off-resonance production of the hadronic final 
states is much more significant in pp annihilation in comparison with e^e~. 

In order to estimate the signal cross-section we calculate the value of the Breit-Wigner 
formula at the resonance energy, E^: 

CTp = —BppBrj^^, (4.7) 

where fc^ = {E]^ — Am^) and the -B's represent the branching ratios of he into the initial 
and final states, respectively. The antiproton momentum at resonance pr = 5.609 GeV/c 
fits both low- and high-resolution ranges of the HESR synchrotron operation (see Sec. |2.2 ). 

Using the value measured by E835 [6J TppBrj^^=10 eV and assuming a value of 0.5 MeV 
for the he width jl^ we obtain dp =33 nb. 

Physics performance of the PANDA detector was studied in the present Thesis for 
several selected decay modes of r]c (77, K^K^tt^ , K*K* and 00). 



4.3 he — > 37 decay mode 

This decay mode was observed at Fermilab by the experiment E835 ^ . It is characterized 
by a fairly clean final state, but the low value of the r^c ~^ 11 branching ratio (4.3 ■ 10~^) 



see Table 4.1|) results in a relatively low event rate in comparison with the hadronic decay 

as a 



modes of rjc- The energies of 7's produced in this decay mode are plotted in Fig. |4.1 
function of the polar angle in the laboratory system. The lower band corresponds to 7's 
from the radiative transition he — )■ ?7c7- The upper band of highly energetic photons is due 
to the rjc — )■ 77 decay. One may note, that observation of both decays, requires registration 
of 7's in the energy range from 150 MeV up to 5.5 GeV. 

4.3.1 Background considerations 

The observation of h^ in 37 final state in proton-antiproton collisions by E835 |S] has 
demonstrated that it is feasible to achieve sufficient background suppression for this final 
state, however differences with PANDA in the detector setup require estimation of the 
degree of background suppression also for the planned experiments. The main contribu- 
tors to the background for the 37 final state are 7's from the 7r°, 77 and rj' decay in 77 
decay modes: the loss of one or more 7's outside the detector acceptance or below the 
energy threshold of the electromagnetic calorimeter (EMC), will result in a 37 final state. 



The background channels considered in this analysis are listed in Table 4.2 with the corre- 
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Figure 4.1: Distribution of events on the 7-ray energy - emission angle plane for the he — )■ f/cT 
37 decay. 



spending cross-sections measured by E760 and E835 [H] integrated over the angular range 
\cos{9)\ < 0.6. 

The angular dependence for all the studied background channels is strongly peaked in 
the forward and backward direction, which is typical for two and three meson production 
in antiproton-proton annihilations at energies of interest. For the Monte-Carlo study the 
angular dependence of the cross-sections was parameterized with 6th or 7th order polyno- 



mials in cos(^). As an example we show in Fig. 4.2 the angular distribution for the 
channel and in Fig. 4.3 for the channel vr'^vr^, respectively. 



Channel 


Cross-section, nb 


PP — !■ TT^Tl^ 


31.4 


PP — 7- 7l^'~f 


1.4 


PP — 7- TT^r] 


33.6 


PP — > rji] 


34.0 


PP — 7- TT^r]' 


50.0 



Table 4.2: The main background contributors to he — )• 87 with the corresponding cross-sections 
[7t I14j integrated over the range of angles | cos(^)| < 0.6. 
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cross-section 
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X' I ndf 


19.6/18 


pO 


0.5184 + 0.1053 


P1 


-0.3077 + 0.9066 


P2 


-4.677 + 3.575 


p3 


25.57 + 18.37 


p4 


40.69 + 28.37 


p5 


-190.5 + 106.4 


p6 


-5.561 ± 57.68 


P7 


354.6 ±178.4 
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Figure 4.2: Angular dependence of the tt^^ cross-sections parametrized (solid line) with a sum 
of powers in cos(^) used in Monte-Carlo simulation. The coefficients of the fit together with their 
errors are indicated in the inset. 



The distribution of events in 7-ray energy emission angle plane for the background 
channel pp — )■ tt^tt^ is shown in Fig. 4.4, In contrast to the signal (see Fig. 4.1 ) events cover 



the plane starting from zero energy. The forward cut-off in Fig. |4.4| is caused by event 
generation over limited range of angles dictated by the experimental data available only in 
the range | cos(6')| < 0.6. 



4.3.2 Event selection and background suppression 

The number of Monte-Carlo events used for this analysis for the signal and all the back- 



ground channels is collected in Table 4.3 



The event selection is done in the following steps: 

1. An rjc candidate is formed by pairing two 7's with an invariant mass in the window 
[2.6; 3.2] GeV. The third 7 is added to this pair to form the he candidate. 

2. A 4C-fit to beam energy-momentum is applied to the he candidate and the informa- 
tion on the he and the updated information on the daughter 7's is stored into the 
root ntuple. 
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4.3 he — !■ 37 decay mode 




Figure 4.3: Angular dependence of the vr'^vr*^ cross-sections parametrized (solid line) with a sum 
of powers in cos{6) used in Monte-Carlo simulation. The coefficients of the fit together with their 
errors are indicated in the inset. 
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Figure 4.4: Distribution of events on the energy-emission angle plane for the pp — >• ir^n^ back- 
ground channel. 



4 Monte-Carlo simulation of PANDA physics performance 



115 



Channel 


Number of events 


pp ^ he ^ 37 


20 k 


PP — )■ 71^71^ 


1.3M 


pp — )■ 7r°7 


100 k 


PP — )■ 71*^77 


1.3M 


pp — )■ 7777 


1.3M 


pp — )■ n^T]' 


100 k 



Table 4.3: Number of generated Monte-Carlo events used for different background channels in 
the he —7- 37 analysis. 



3. The following cuts are applied at the ntuple level to suppress background: 

(a) Events with 37's were selected. This cut keeps 47% of the initial events, 

(b) Cut on the confidence level of the 4C-fit: CL > 10"^, 

(c) Cut on the CM energy of the 7 from the he — )■ ?7c7 radiative transition: 0.4 GeV 
<E^< 0.6 GeV. 

(d) Angular cut | cos(^^)| < 0.6, to reject the background which is strongly peaked 
in the forward and backward directions. The cos(^) distributions for the back- 



ground channel (tt^tt^) and for the signal are shown in Fig. 4.14 and Fig. 4.15 
respectively, 

(e) The cut on invariant mass for combinations M(7i, 73) >1.0 GeV and M(72, 73) > 
1.0 GeV (the value of the cut is determined by the rj' mass). 

An impact of the applied selection criteria on the signal and background events is discussed 
below. 



Fig. 4.5 demonstrates the distribution of 7's in multiplicity for the signal events; 53% 
of the events have exactly 3 reconstructed neutral particle candidates. A small fraction 
of events has less than 3 neutral candidates due to detector acceptance. The events with 

> 3 are caused by the electromagnetic split-offs. For pp vr^vr" the distribution of 7's 



in multiplicity is shown in Fig. 4.6 8% of these events have three neutral candidates. One 
may conclude, taking into account the ratio of signal to background cross-sections, that 
the number of background events will exceed the number of expected 37 events for the 
signal. A similar situation is encountered for all the remaining background channels listed 



in Table 4.2 , i.e. they have a significant fraction of 37 events. 
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Figure 4.5: Multiplicity of the reconstructed 7's in the he ??c7 ^1 decay. 
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Figure 4.6: Multiplicity of the reconstructed 7's (73) for the pp — >■ -k^tt^ background reaction. 
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The distribution of confidence level for the 4C-fit to beam energy-momentum for the 



he Vcl 37 events is presented in Fig. 4.7 on the linear and in Fig. 4.8 on the logarithmic 
scale of the abscissa in order to expand the region with CL < 0.1. The cut is applied at 
the level CL > 10~'^, as indicated above in the step 3b. Lower values of CL correspond 
to higher of the fit. The CL distribution is almost fiat in the range [0.1; 1.0] for the 
signal events, whereas for the background channels the probability of higher CL values is 
lower in general as it can be seen from Fig. 4.8 for the pp vr^vr^ background channel. For 
many background events the 4C-fit does not converge. Such events get assigned negative 
CL values and are removed by the CL > 10~^ cut from subsequent analysis. 

The 7's which are combined to rjc are labeled 1 and 2, whereas number 3 is assigned 
to the 7 emitted in the radiative transition he — )■ ?7c7- The distribution of reconstructed 
energies of the 3^^^ 7 in the pp CM system is presented in Fig. 4^ The expected energy of 
the 7 emitted in the he — )■ ?7c7 transition should be 503 MeV according to the compilation 
[1]. The reconstructed distribution is peaked close to this value. The range [0.4; 0.6] GeV 
between the two vertical lines defines the limits of the effect (cut 3c). The corresponding 
distribution for 7's emitted in pp — )■ tt^h^ is presented in Fig. 4.10 40 % of the events for 
this particular background channel pass the latter selection cut. 

The distribution of events in the i?^3-invariant mass m{rie) plane is presented in Fig. 



4.11 



for the he — )■ ?7c7 channel. All the events are concentrated along a line which is defined by 
the applied 4C-fit, consequently, after the fit the E^s and m{rie) variables are no longer inde- 
pendent. The size of boxes is proportional to the number of events in a certain E^^, m{rie) 
range. One may note that events are mostly concentrated around the point £^^3=0.503 
GeV, m{rie)=2.98 GeV. For the case of pp — )■ tt^tt'^ background reaction (see Fig. 4.12) 
events are more uniformly distributed along the same line. 

Fig. 4.13 presents the reconstructed distribution of the r]c invariant mass (projection 
of Fig. 4.11 on the Mirje) axis). A Breit-Wigner formula, fitted to the simulated data, is 
superimposed on the figure. The fitted width parameter F = 28 MeV is slightly larger than 
the PDG value Tpj:,G=yj MeV, which illustrates the influence of instrumental resolution, 
Ae^p ^ 22 MeV, expected with the PANDA-EMC. 



It was mentioned in Sec. 4.3.1 , that the background cross-sections and the corresponding 
intensities of 7's are peaked in the forward and backward directions. The distribution of 
the reconstructed 71^2 in cos(^) in the CM system is shown in Fig. 4.14 for the pp — 
Ti^Ti^ channel. An abrupt drop in intensity above | cos( 



> 0.6 reflects the limits of the 
generated events, because of the missing experimental information (see the remarks related 
to Fig. 4.3[ ). Vertical lines at | cos(^)| = 0.6 indicate the width of the window. In case of 
the signal (Fig. 4.15) events have rather flat distribution in cos(^), therefore applying a 
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Figure 4.7: Distribution of confidence level of 4C-fit for he — )• t/c7 ~^ 37. 



cut eliminating forward and backward angles is a priori expected to improve the signal to 



background ratio. The two dips seen in Fig. 4.15 are related to the EMC installation: the 
one at ^ 0.35 reflects irregularity between the barrel and the forward endcap, the other 
one at ~ 0.95 is at the angle of transition between the forward endcap and the shashlyk 
calorimeter in the Forward Spectrometer. 



Fig. |4.16| presents the Dalitz plot for 37 events, which pass all the previously mentioned 
cuts. The signal events are concentrated along the diagonal extending from the top left to 
the bottom right. Background events from different reactions are marked as red dots. The 
blue lines mark the selection cuts, i.e. M(7i,73) > l.OGeV^ and M(72,73) > l.OGeV. The 
value of 1 GeV is chosen to eliminate events originating from the pp — > TT^rj' background 
channel, i.e. by the mass of 77'. None of the background events passes this selection cut and 
only 37% of the signal events is able to survive, which results in the total signal efficiency 
of about 8%. 



4.3.3 Results 



The selection efficiencies for different cuts are collected in Table l44l Efficiencies are cumula- 
tive, i.e. applied one after another. Taking into account the signal cross-section crpp_i./i^=33 
nb at resonance, branching ratio BR{r]c — ?■ 77) = 4.3 ■ 10~^ and background cross-sections 



from Table 4.2 one gets the expected signal to background ratios as listed in Table 4.5 
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Figure 4.8: Distribution of confidence level of the 4C-fit for he 
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background events on a logarithmic scale of the abscissa. The vertical line indicates the assumed 
position of a cut on CL > 10~^ (see Sec. 4.3.2). 
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Figure 4.9: Distribution of reconstructed energy of the 7 (73) from the he — )• r/c7 radiative 
transition. Vertical lines indicates the window imposed to select this transition. 
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Figure 4.10: Distribution of reconstructed energy of the 73 emitted in the background channel 
pp — )• TT^TT^. Events between the vertical lines would be erroneously assigned to he — )• ?7c7 if not 



eliminated by other selection criteria (see Sec. 4.3.2). 
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Figure 4.11: Distribution of reconstructed energy of 7 (73) versus r]c invariant mass for he — )• ?/c7- 
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Figure 4.12: Distribution of reconstructed energy of 7 versus ijc invariant mass for the pp 
background reaction. 
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Figure 4.13: Reconstructed invariant mass of rjc- Solid line is the Breit-Wigner fit with the 
indicated parameters. 
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Figure 4.14: Distribution of the reconstructed cos 6 of the 7's assigned to rjc — s- 77 decay in CM 
system for the pp — t- vr^vr'^ background reaction. 
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Figure 4.15: Distribution of the reconstructed cos 9 of the 7's assigned to 7]c — )• 77 in CM system 
for the he — )• r/c7 decay. 
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Figure 4.16: Dalitz plot for the pp — t- 777 events. Vertical lines denote the positions of cuts to 
eliminate tt'^, tt't] and tt't]' events originating from the background reactions listed in Table 4.4 



The expected event rate for running in high luminosity mode, L = 2 ■ 10"^^ cm^^s^^, is 
20 events/day, and for high resolution mode with L = 2 ■ 10'^^cm~^s~^, 2.0 events/day, 
respectively. 
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Table 4.4: Selection efficiencies for he — ?• 87 and the indicated background reactions. 
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4.4 he — )■ 007 decay mode 



Channel 


S/B ratio 


pp — 7- TT^TT^ 


> 94 


PP — i- TT^'-f 


> 164 


PP — i- TT^f] 


> 88 


PP — r7?7 


> 87 


PP Tc^r]' 


> 250 



Table 4.5: Signal to background ratio for he — t- 87 and different background cliannels. 

4.4 he ^ (fxf)^ decay mode 

As a benchmark channel with a hadronic decay mode of the t]c the (jxp final state was 
studied with BR = 2.6 • 10^'^. The (p was detected through the decay (j) — )• K^K^ , with 
BR =0.49 H]. 

4.4.1 Background considerations 

For the exclusive decay mode considered in this study: 

pp ^ he ^ rje'^l —7- 007 — 7- K^K^K^K^'j, 
the following three reactions are considered as the main contributors to the background: 

1. pp -> K+K-K+K-TT^, 

2. pp^ (pK+K^Ti^, 

3. pp — 7- (f)(j)7[^. 

With one photon from the 7r° decay left undetected, these reactions have the same final 
state particles as the studied he decay. 

Additional possible sources of background are: 

4. pp A++(1232)A""(1232)7r° pn+pn-n^, 

5. pp — 7- K^K^n^Ti^TC^, 

6. pp — )■ TT + TT^ 71 + 71^71*^. 
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Figure 4.17: Energy distribution of the 7 for the — )• /ic — )• ??c7 and pp — )• (fxpTT^ reactions. 



The reaction 7^4 has been considered as one of the main sources of background in the 
Fermilab experiment E835 |6], which studied the same decay mode of he. This reaction has 
similar kinematics as the reaction of interest, i.e. 4 charged tracks in the forward direc- 
tion. PANDA has significant advantages for discrimination of this source of background in 
comparison with E835, being equipped in magnetic analysis and PID tools. Nevertheless, 
an attempt was made to estimate explicitly the suppression factor for the reaction #4. 

The last two reactions could contribute to background because of pion misidentification 
as kaons. In addition to PID such events can be effectively suppressed with the 4C-fit to 
beam energy-momentum. However, because of large cross-sections for these reactions, it 
seemed necessary to estimate the signal to background ratio for these channels in a detailed 
simulation. 



Fig. |4.17| presents the distributions of 7's in the laboratory energy for the signal and 
the background channel #3. 

The energy range of 7's from the he decay is [0.15; 2.0] GeV. The distribution for the 
background has different features. First, it is wider distributed extending beyond this 
energy range. Second, it has an increasing tendency towards zero energy. If we want to 
recover vr^'s to separate the signal from background, we should succeed to lower photon 
detection threshold as much as possible. Moreover, when one 7 from the 7r° decay has low 
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energy, the other 7 together with the charged hadrons has total momentum closer to the 
total momentum of the initial pp system. Such events pass cuts on the fit probability of 
the 4C-fit; this situation increases the relative importance of having low energy threshold 
in 7-ray detection in order to apply a veto on vr^'s in an event and correspondingly to 
suppress background due to vr^'s. 

There are no experimental measurements, to our best knowledge, of the cross-sections 
for the first three background reactions, which are supposed to be the main contributors 
to background. The only way to estimate their cross-sections was found to use the DPM 
(Dual Parton Model) event generator [39]; 2 ■ 10^ events were generated with DPM at the 
p beam momentum pz = 5.609 GeV/c, which corresponds to the studied he resonance. The 
corresponding numbers of events are 115 and 12 for the first two background channels. No 
events for the pp — )■ (jxpir^ reaction were observed in the generated sample. With the total 
pp cross-section at this beam momentum of 60 mb, the cross-sections for the corresponding 
background channels are estimated at 360 nb, 36 nb and below 6 nb, respectively. For 
the pp — )■ A+"'"(1232)A (1232)7r° and pp — > 7r+7r^7r+7r^7r° background channels the cross- 
section are known from the measurements at the p beam momentum p = 5.7GeV/c [17] 
and [H] . The values are equal a = 530 fib and a = 750 fib, respectively. For the purpose of 
estimation of the signal to background ratio it was assumed that the phase space engaged 
in the final state is uniformly populated. An estimate of the pp — > K^K^ 11^71^71^ cross- 
section was done by extrapolating from a lower energy according to the total inelastic 
cross-section. The result is a = 30 fib. 



4.4.2 Analysed events and selection criteria 



The numbers of analysed events are listed in Table 4.6 For the pp — > K^K^ii^n^'n-^ 
channel 15 million out of the 20 million of the events were simulated with a filter applied 
on the invariant mass of the pair of kaons. The events with m{K^ K~) in the range [0.95; 
1.2] GeV were selected. The efficiency of the filter is 29.9 %, i.e. only 29.9 % of the events 
produced by the event generator are simulated with GEANT4 transport code with the 
following event reconstruction. This gives the effective number of simulated events ~ 55 
million. 

The following selection criteria were applied: 



1. candidates were defined as , K pairs with the invariant mass in the window 
[0.8; 1.2] GeV. Two candidates in one event with invariant mass in the window 



4 Monte-Carlo simulation of PANDA physics performance 



127 



Channel 


N of events 


PP ^ he 4>4>'y 


20 k 


PP K+K-K+K'Ti^ 


6.2M 


PP — )■ (pK^K^Tx^ 


200 k 


PP — i- (j)(j)TX^ 


4.2 M 


PP A++(1232)A""(1232)7r° 


100 k 


PP — )■ K^K^n^Ti^TT^ 


5 M + 15 M 


PP — )• Tr^TT^TT^Tr^Tr'^ 


1 M 



Table 4.6: The number of analysed events for the he decay and the selected background reactions 



[2.6;3.2] GeV defined an rjc candidate which, combined with a neutral candidate, 
formed an h^ candidate. 



2. A 4C-fit to beam energy-momentum was applied to the he candidate, which was 
stored to a root ntuple together with the updated information on its decay products. 

3. The following additional cuts are performed at the ntuple level for background sup- 
pression: 

(a) cut on the confidence level of the 4C-fit to beam energy-momentum, CL > 0.05, 

(b) ?7c invariant mass within [2.9; 3.06] GeV, 

(c) within [0.4; 0.6] GeV, 

(d) invariant mass within [0.99; 1.05] GeV, 

(e) no TT^ candidate in an event, i.e. no 27 invariant mass in the range [0.115; 0.15] 
GeV with two different low energy photon thresholds: 30 MeV and 10 MeV. 



Fig. 4.18 presents the multiplicity distribution of reconstructed EMC clusters for the 
signal and one of the background channels. One may note that the mean number of neu- 
tral candidates exceeds one, the value expected for the signal, or two expected for the 
background from 7r° decay. This is the result of hadronic split-offs, caused by nuclear 
interactions of charged hadrons with the material of the detector. Secondary particles 
produced in such interactions can create additional clusters in EMC, which cannot be as- 
sociated with primary charged hadrons, and thus are interpreted as neutral candidates. 
This effect makes it impossible to select as the signal the events with only one cluster. 
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Figure 4.18: The number of reconstructed EMC clusters for the — )• /ic — )• ?7c7 and pp 
reactions. 



because it leads to a significant drop in efficiency. This observation emphasizes the impor- 
tance of other selection criteria, in particular of the veto on 7r° in an event. The effect of 
the latter requirement strongly depends on the assumed low energy photon threshold. 



The distribution of - K pair invariant mass is presented in Fig. 4.19 for the signal 



and in Fig. 4.20 for the pp — K^K K^K vr'' background. Windows imposed on the mass 
of are marked with vertical lines (in red) . 



4.4.3 Signal to background ratio and low energy 7-ray threshold 



The efficiencies of various cut selection criteria are listed in Table 4/7 for the signal and 
the first three of the considered background channels. 

For pp —7- A++(1232)A (1232)7r'^ and pp — t- 7r+7r~7r+7r~7r'' no event has been found to 
pass pre-selection. 

Assuming the he production cross-section of 33nfe at resonance, one obtains the signal 
to background ratios given in Table |4.8[ 

For the pp A+''"(1232)A (1232)7r° and the pp — )■ 7r"'"7r~7r"'"7r~7r° channels, because 
of their large cross-sections, one needs around 10* events to make valid conclusions on 
background suppression at S/B levels above 1:1, which is not easy to fulfil because of the 
required excessive CPU time. The pp — )■ K~ tx^ tt~ ir^ channel is more favorable from this 
point of view, having an order of magnitude smaller cross-section. Fortunately, it permits 
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Figure 4.19: Invariant mass of K^K^ with mass window on (j) for he — )• (jxpl- 
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Figure 4.20: Invariant mass of K^K witii mass window set on (j) for the pp — )• K^K K^K vr*^ 
background. 
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Table 4.7: Efficiency of different event selection criteria. 



channel 


Signal/background ratio 


PP — 7- (fxp-K^ 

PP — )■ K'^K^TT^'K^TT^ 


8 
8 

> 10 

> 12 



Table 4.8: Signal to background ratio for different he background channels. 



to study the influence of same factors, i.e. PID misidentification and power of the 4C-fit 
on background suppression. For the latter channel 2 ■ 10'' events have been analysed and 
the achieved conclusion tells that the expected S/B ratio will exceed 12:1 (see Table 4.8). 

For the pp — ?■ (fxpn^ background channel the reduction of low energy 7-ray threshold 
from 30 MeV to 10 MeV gives 20 % improvement in the signal to background ratio, for the 
PP — )■ (f)K~^K~TT^ the corresponding improvement is 40%. 

It is worthwhile to note that finite energy resolution is taken into account in the digiti- 
zation of EMC signals which results in around 14 % at 10 MeV in Monte Carlo simulations 



using the adopted software. For comparison the measurements presented in Sec. 3.3.4 give 
around 17 % at this energy. However the readout with two APDs considered for PANDA 
will increase the collected light and improve resolution especially at these low energies 
where the statistical term is of higher importance. Taking this remark into account one 
may conclude that the resolution used in Monte Carlo simulations is consistent with the 
above presented measurements. 

With the final signal selection efficiency of 24% (see Table 4.7) and the assumed lumi- 
nosity in high luminosity mode of L = 2 ■ lO^^cm-^s-^, the expected signal event rate is 92 
events/day. For the high resolution mode with L = 2 ■ lO^^cm-^s-^, the expected signal 
event rate is 9 events/day, respectively. 
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4.5 he ^ KgK^TT 7 decay mode 

The decay mode rjc — ?• KsK^n~ has the highest branching ratio among the decay modes 
of rjc (see Table 4.1). The presence of Ks in the final state allows a clear selection of this 
channel using the geometrical constraint on the secondary decay vertex. The lifetime of 
Kg corresponds to cr = 2.68 cm. The dominant decay mode of Kg is tt^tt^ with 

branching ratio BR=68.95 %. The complete reaction and decay path is as follows: 

PP ^ he ^ TjcPj — > KsK^Tl^^ —7- 71^71^ K^TT~'J. 



4.5.1 Background considerations 

The main source of background considered in this study is pp — )■ KsK^tt^tt^ , i.e. it 
differs from the studied signal by 7 replaced with tt^. Such a source of background is 
very typical for the study of electromagnetic transitions, where non-resonant background 
contains the same hadrons in the final state and one 7 from vr^ decay mimics the 7 from 
the considered electromagnetic transition. The background cross section is unknown at the 
beam momentum p^ = 5.609 GeV/c, however it was measured at lower beam momentum 
Pz = 3.66 GeV/c in 1966 in Brookhaven [TT]. It is equal 74 and a rough estimate of the 
cross section at pz = 5.609 GeV/c was obtained by scaling it down with the total inelastic 
cross section. The obtained value is 60 fib. 

Additional sources of background are K^K^tt^ti^ti^ and Tr+Tr^Tr+Tr^Tr" assuming PID 
misidentification of tt and K. However it will be shown below that the background channel 
pp — )■ KsK^TT'Tc^ cannot be sufficiently suppressed and these two additional channels have 
not been studied in detail. 



4.5.2 Event selection and background suppression 

The analysis has been performed with 80 k of the signal events pp ^ he KsK^tt"'-/ and 
100 k of the background events pp — )■ KsK^h'tt^ . 

The event selection was done in the following steps: 

1. Ks candidates were selected combining two charged vr with the requirement of the 
common vertex and mass window [0.3; 0.8] GeV. Ks were combined with K^ and tc" 
and provided their invariant mass fallen into the window [2.6; 3.2] GeV they defined 
an rjc candidate, r/c combined with a neutral candidate, formed an he candidate. 
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4.5 he KgK^Ti 7 decay mode 



Selection criteria 


signal 


KsK^ 


^Tr- 




pre-selection 


0.41 


2.82 


io 


~2 


CL > 0.001 


0.27 


4.29 


10 


-3 




0.26 


1.17 


10 


-3 


no 7rO(30Me1/) 


0.22 


7.90 


10 


~4 


no 7r°(10Me1/) 


0.20 


5.00 


10 


~4 



Table 4.9: Efficiency of different event selection criteria for KsK^tt 7 decay mode of he- 



2. A 4C-fit to beam energy-momentum was applied to the he candidate, with the result 
stored to a root ntuple together with the updated information on its decay products. 

3. The following additional cuts are performed at the ntuple level for background sup- 
pression: 

(a) cut on the confidence level of the 4C-fit to beam energy-momentum, CL > 0.001, 

(b) r]c invariant mass within [2.9; 3.1] GeV, 

(c) within [0.4; 0.6] GeV, 

(d) no 77° candidates in the event, i.e. no 27 invariant mass in the range [0.115; 0.15] 
GeV with two different low energy photon thresholds: 30 MeV and 10 MeV. 



4.5.3 Results 

The efficiency of different selection criteria for the signal and background are presented in 
the Table SH 

The signal to background ratio 1:240 has been obtained with the signal efficiency of 
20 %, which is prohibitively small, in our point of view, to perform this experiment. The 
signal to background ratio is improved by 20 % by reducing the 7-ray registration threshold 
from 30 MeV to 10 MeV, however this does not reverse the conclusion on difficulties with 
he registration in this final state. This estimate is based on the assumption that the 
background cross section decreases with energy as the inelastic pp cross section. If it drops 
faster the signal to background ratio can be more favorable for performing this experiment. 

Additional explanation of the difficulty with signal and background separation follows 



from the inspection of Fig. |4.21[ This Fig. demonstrates the 2-dimensional distributions of 
the signal and background events in the m{ric) vs plane. The panels a) and b) contain 
the events coming directly from the event generator, i.e. with no detector effects included. 
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Figure 4.21: Signal and background event distributions in the m{ric) vs plane for the direct 
output of the event generator [a) and b)] and after full Monte Carlo analysis [c) and d)]. 



The panel b) is an expanded part of a) containing the region of signal events. One is prone 
to conclude from Fig. 4.21 b) that pure kinematics might suffice to separate the signal 
and background events. However, the finite detector resolution smears the region that 



corresponds to the signal so that it overlaps strongly with the background [see Fig. 4.21 
d)]. 



4.6 he K^*K^*j decay mode 

The analysis of this decay mode is very similar to the two hadronic final states presented 
above. The branching ratio for the rjc — ?■ K'^*K'^* decay mode is BR = 8.5 ■ 10~^. The main 
decay mode of the K^* resonance is K^* K^n^ with the branching ratio BR=0.67. 
The complete reaction and decay path is: 

4.6.1 Background considerations 

The two background channels were considered for this decay mode of he- 
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4.6 he i^°*ii'°*7 decay mode 



Channel 


N of events 


pp K^*K-n+TT^ 
pp — )■ K~^K^n^7i^7i^ 


80 k 
1.1 M 
200 k 



Table 4.10: The numbers of analysed events for he K^*K^*^ decay and background reactions 



pp — K^K TT + TT 7r°, 

and 

i.e. one intermediate K'^* resonance is engaged in the latter background channel. The 
background cross sections were estimated in the same way as for the KsK'^iT~'~f decay 
mode. The measurements at the beam momentum p = 3.66 GeV/c give o"pp_^^+^-^+^-^o = 
Sifib [TT]. Scaling down according to the inelastic cross-section gives a = 30 fib at p = 
5.6GeV/c. At the 3.66 GeV/c antiproton beam momentum 30 % of reactions goes via one 
resonance. Assuming the same fraction at 5.6 GeV/c momentum one gets an estimate 
(Tpp — )• K^*K~K^Ti^ = 8.4 fib. No evidence for K^*K^* emission was observed in [TT]. 

4.6.2 Event selection and background suppression 

The number of generated events for the signal and two background channels is collected in 
Table lilOl 

The event selection was done in the following steps: 

1. K*^ candidates were selected combining and vr" pairs with the requirement for a 
common vertex and invariant mass in the window [0.7; 1.1] GeV. K*^ is produced in 
the charge conjugated combination within the same mass window. The pair of K*'^ 
and K*^ with invariant mass in the window [2.6; 3.2] GeV defines an rjc candidate. 
7]c combined with a neutral candidate, formed an h^ candidate. 

2. A 4C-fit to beam energy-momentum was applied to the he candidate, which was 
stored to a root ntuple together with the updated information on its decay products. 

3. The following additional cuts are performed at the ntuple level for background sup- 
pression: 
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Selection criteria 


signal 








K+K- 


'7T + 7T 


pre-selection 


0.43 


1.4- 


10- 


^2 


6.6 


■ 10-3 


CL > 0.1 


0.31 


2.3 ■ 


10- 


'3 


1.0 


■ 10-3 




0.27 


3.5 ■ 


10- 


-4 


1.9 


■ 10-4 




0.13 


4.5 ■ 


10- 


-5 







no 7r°(30 MeV) 


0.11 


2.4- 


10- 


-5 







no 7r°(30 MeV) 


0.10 


1.3 ■ 


10- 


-5 








Table 4.11: Efficiency of different event selection criteria for the K°*K^*^ decay mode of he- 

(a) cut on the confidence level of the 4C-fit to beam energy-momentum, CL > 0.1, 

(b) r]c invariant mass within [2.93; 3.03] GeV, 

(c) within [0.4; 0.6] GeV, 

(d) invariant mass of and within [0.86; 0.94] GeV mass window. 

(e) no vr'' candidates in the event, i.e. no 27 invariant mass in the range [0.115; 0.15] 
GeV with two different low energy photon thresholds: 30 MeV and 10 MeV. 

4.6.3 Results 

The efficiency of different selection criteria for the signal and background are presented in 
the Table KTT\ 

The obtained signal to background ratio is 1:10 for the K'^*K-'n-+7T° and for the 
K+K-n+n-TT^ channel the limit > 1:15 was established. So the conclusion is that this 
decay mode presents difficulties for an experimental observation. 

4.7 Determination of the he width 

In this section we investigate the performance of PANDA in determination of the width. 
For this purpose we performed Monte Carlo simulations of energy scans around the res- 
onance. Events were generated at 10 different energies around the he mass, each point 
corresponding to 5 days of measurements in high resolution mode. 

The expected shape of measured resonance in pp ^ ^ ?7c7 is the convolution of the 
Breit-Wigner resonance curve with the normalised beam energy distribution and an added 
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4.7 Determination of the he width 



Tji^MC, MeV 


Tu^reco, MeV 


ATr, MeV 


1 


0.90 


0.22 


0.75 


0.71 


0.17 


0.5 


0.52 


0.14 



Table 4.12: Reconstructed he width (Tji^reco) by fitting to the "experimental points" generated 



in Monte-Carlo simulations described in Sec. 4.7, assuming that the width is Tr^mc- 



background term. The expected number of events at the z-th data point is: 

= X y Ldt], X + X j ^^-^^^-^^-^^iEl (4.8) 

where cxj is the beam energy resolution at the i-ih. data point, and Mr the resonance 



width and mass, ap incorporates branching ratios for the formation and decay (see Eq. 4.7), 



the factor in square brackets in front of the r.h.s. of Eq. 4.8 is the product of e, an 
overall efficiency and acceptance factor and the integrated luminosity at the z-th point 
of measurements. To extract the resonance parameters the likelihood function, -ln£, is 
minimized assuming Poisson statistics, where: 



For our simulation we assumed a signal to background ratio of 8:1 and we used the signal 
reconstruction efficiency of the he — )■ ?7c7 — ^ 007 channel (see Sec. 4.4). The simulated 



data were fitted to the expected signal shape with four free parameters: Er, F/j, Ohkgd-, 
CTp. The background was assumed energy independent. The study has been repeated for 
three different Yr = 0.5, 0.75 and 1.0 MeV. The results of the fit for 0.5 MeV and 1.0 
MeV are presented in Fig. 4.22 and Fig. 4.23, respectively. The extracted Trs (Tji^reco) 
with errors (AF/j) are summarized in Table 4.12 As a conclusion one may state that with 
the designed parameters of the beam energy resolution the precise determination of the 
width of the he resonance is feasible with PANDA. However the scan should be performed 
in High Resolution mode of HESR which provides lower luminosity and therefore is rather 
time consuming. 
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E, GeV 



Figure 4.22: Fit of the Breit-Wigner distribution convoluted with gaussian resokition of the 
beam momentum (soUd hne) to the "experimental points" generated assuming that the he reso- 
nance width is r = 0.5 MeV. 




Figure 4.23: Fit of the Breit-Wigner distribution convoluted with gaussian resolution of the 
beam momentum (solid line) to the "experimental points" generated assuming that the he reso- 
nance width is r = 1.0 MeV. 
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4.8 Reconstruction of the angular distribution of photons in the he decay 



4.8 Reconstruction of the angular distribution of pho- 
tons in the he decay 

From the measurements of he in the CLEO experiment the decay of he is consistent 
with the hypothesis of spin-parity J^'" = l'^" , however with rather high uncertainty. The 
present Monte Carlo study has been performed to determine the minimum duration of 
an experiment with the PANDA detector which would determine angular distribution of 
photons with sufficient precision that could serve a proof of J^^ = 1+^. 

To study the feasibility of reconstruction of the angular distribution in he — )■ ?7c+7 decay, 
the proper angular distribution was implemented into Monte Carlo simulations. EvtGen 
event generator provides the possibility to define angular distribution of the decay. The 
W{9) ~ sin^ 9 dependence was implemented for pp ^ he ^ rje + •y reaction which follows 



from Sec. 1.4 , where 9 is the polar angle of t]c emission with respect to the proton momentum 



in the pp center-of-mass system. The rje — )■ 00 channel was selected for reconstruction of 



the angular distribution as the most promising for he detection (see Sec. 4.4) from the 
point of view of signal/background ratio. It should be taken into account that efficiency of 
the EMC is not uniform in cos 9 to reconstruct properly the angular distribution. Events 
with phase space uniform distribution for he ^ rjc + j decay have been generated and 
reconstructed to determine correction factors taking into account deviations from uniform 
detection as a function of the polar angle. Efficiency of he reconstruction as a function of 



cos^ is presented in Fig. 4.24 The characteristic drop in efficiency can be seen around 
cos^ = 0.4. This region corresponds to the transition between the barrel EMC and the 
forward endcap. 

The reconstructed angular distribution of 7-rays from the pp he ^ rje + '~f reaction, 



corrected for acceptance of the detector, is presented in Fig. |4.25| as a distribution in cos 6*. 
The presented data points correspond to 10 days of running with the estimated production 
rate of 92 events/day in the High Luminosity mode. The entire angular range in cos 6* 
was divided into 20 intervals. The fit with W{9) = sin^ 9 hypothesis gives x^/^d.oj = l-l- 
We conclude that the accumulated statistics is sufficient to identify J^*" = l+~ for the he 
resonance, uniquely. 
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cos(e) 

Figure 4.24: Efficiency of he reconstruction in the 20 intervals of cos 6. 




cos(e) 



Figure 4.25: The distribution of photons emitted in — )• /ic — )• rye + 7 in the 20 intervals of 
cos 6, assuming J^*^ = 1"* for he, which leads to El emission in the transition he — ?• tJc- 



Chapter 5 
Conclusions 



The PANDA experiment at the FAIR facihty is in the stage of finahzing its design with 
the beginning of construction expected in the nearest one or two years. With the start 
of operation expected in 2015 it will contribute to many topics in hadron physics and the 
high precision charmonium spectroscopy is a one of the major tasks of PANDA. Electro- 
magnetic transitions to lower lying states are among the main decay modes with the aid of 
which the states of charmonium can be identified. This emphasizes the importance of an 
electromagnetic calorimeter (EMC) as one of the main components of the PANDA detec- 
tor. The electromagnetic calorimeter of PANDA will be based on lead tungstate (PWO) 
scintillators readout with avalanche photodiodes (APDs). 

In this thesis the results of measurements of energy resolution of lead tungstate scin- 
tillators readout with APDs are presented. Energy resolution was measured in the energy 
range from 6 to 17 MeV and around 25 MeV, which is close to the low energy threshold 
assumed for the PANDA EMC. Two series of measurement have been performed: using 
protons and low energy 7-rays produced in radiative capture reactions. 

Protons with energy around 25 MeV accelerated with the C-30 cyclotron, were used to 
measure energy and time resolution of small size PWO-II scintillators 20 x 20 x 40 mm^ 
readout with several types of avalanche photodiodes. One of the purposes of these mea- 
surements was to gain a practical experience of work with this scintillator having a rather 
low light yield and to draw conclusions on how the readout performance is improved when 
the APD is cooled to -20°C. The energy resolution measured with two types of API APDs 
17 % and 22 % appeared better then 35 % measured with the Hamamatsu APD, which is 
to large extent an effect of twice smaller active area of the latter. However the Hamamatsu 
S8664-1010 APD was selected for the PANDA experiment because of its mechanical design 
and lower operating voltage are better fitting large scale applications. Simulation of light 
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propagation in a scintillator and its APD detection performed with the Litrani program 
showed that the setup was not optimum from the point of view of light collection, which 
resulted in high contribution of the statistical term to energy resolution. At the same 
time the time resolution below 1 ns was demonstrated for PWO-II scintillator with APD 
readout. The results of these measurements were taken into account during construction 
of the setup for measurements of energy resolution with low energy photons produced in 
radiative capture reactions. 

The setup was constructed to measure energy resolution of cooled PWO scintillator 
with Hamamatsu APD readout at the three discrete 7-ray energies of 6.13 MeV, 12.8 
MeV and 17.2 MeV, the latter two produced in 7)^^C reaction and the former one 

in ^^F{p,a)^^0* on the target contaminant. The crystals of the size 20 x 20 x 200 mm^ 
were used, which is close to the final design of crystals for the PANDA EMC. The goal 
temperature — 25°C has not been achieved due to limitation of the cooling capacity of the 
used JULABO F32-ME Refrigerated/Heating Circulator and a measurement at — 20°C 
has been performed. The energy resolution 19 ± 3%, 15 ± 2% and 9.6 ± 1.4% has been 
obtained for the three energies, respectively. These confirm nicely that the smooth energy 
dependence of resolution, established in the energy range 40.9-674.5 MeV by the experiment 
performed at the MAMI facility in Mainz (Germany) [Eq. 3.8 ] gives a valid extrapolation 
of resolution into this low energy range. 

The Monte Carlo simulations aimed to study the performance of the PANDA detector 
in detection of the charmonium state demonstrated feasibility of the study of this state 
with high precision with PANDA. With the main decay channel /ic — )■ ?7c + 7 it can be 
registered exclusively in either neutral or hadronic final states of rjc- Several possible decay 
modes of t]c have been considered {KsK~^tt~ , K*K*, cfxp and 77). It was demonstrated that 
the most promising decay mode for hf. observation is 007, where signal to background ratio 
around 8:1 is expected. With the signal selection efficiency of 24% around 92 events per 
day are expected in the high-luminosity mode. An alternative channel for he registration 
is a neutral final state /ic — )■ + 7 — ^ 777? however with lower predicted event rate 
due to small efficiency and low branching ratio of the r/c 11 decay. A study of the 
PANDA sensitivity in determination of the h^. width has been performed and estimates 
of precision and needed running time for the three hypotheses concerning the h^. width: 
r=0.5 MeV, 0.75 MeV and 1 MeV were obtained. The feasibility of reconstruction of the 
angular distribution of the decay photons with the required precision for confirmation 
of the spin-parity assignment J^^ = 1^^ was also demonstrated. 
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